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CHAPTER  I 
INTRODUCTION 


The  concept  of  shearing  maps  as  generalizations  of  Dehn  twists  and  transforma¬ 
tions  on  points  in  Teichmiiller  space  was  first  developed  by  Thurston  [49,  51],  and 
has  been  used  by  Kerckhoff  in  proving  the  Nielsen  Realization  Conjecture  [32]  and 
Bonahon  in  his  investigations  into  Thurston’s  boundary  of  Teichmiiller  space  [6]. 
These  earthquakes  have  been  extended  by  McMullen  [35]  to  include  another  geologic 
action,  that  of  grafting,  on  points  in  Teichmiiller  space.  Taken  together,  through 
composition,  we  have  complex  earthquakes. 

Essentially,  Thurston’s  shearing  earthquake  opens  a  hyperbolic  surface  (point  in 
Teichmiiller  space)  along  a  measured  geodesic  lamination,  shears  the  surface  along 
the  seam  in  a  manner  determined  by  the  measure,  and  re-attaches  the  surface.  Mc¬ 
Mullen’s  grafting  maps,  rather  than  shearing  along  the  lamination,  insert  or  remove 
a  cylinder  and  re-attach  the  surface.  In  Chapter  II,  we  describe  these  earthquake 
actions  in  detail,  giving  several  examples  paying  particular  attention  to  the  geometric 
nature  of  the  transformations. 

A  circle  packing  is  a  collection  of  circles  with  a  prescribed  pattern  of  tangencies. 
Of  course,  given  a  collection  of  circles,  there  is  no  guarantee  that  they  can  be  “fit 
together”  in  a  particular  pattern;  in  fact,  some  patterns  may  not  be  possible.  For 
example,  given  that  all  the  circles  have  the  same  radius,  one  circle  cannot  be  tangent 
to  more  than  six  other  circles.  The  “prescribed  pattern,”  is  a  strictly  combinatorial 
structure  with  no  inherent  geometry.  As  the  circles  adjust  their  radii,  trying  to 
meet  the  constraints  of  the  pattern  prescribed,  a  rigid  geometry  is  realized.  The 
interplay  between  the  combinatorial  structure  and  the  rigid  constraints  inherent  in 
the  geometry  of  the  circles  provides  a  deep  link  to  geometric  function  theory  and 
the  structure  of  Riemann  surfaces  which  we  will  exploit  in  describing  the  relationship 
between  earthquakes  and  circle  packing. 

Imposing  the  geometry  of  a  circle  implies  the  existence  of  a  metric;  thus  we  can 
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speak  about  circle  packings  on  any  surface  with  a  metric.  Consequently,  we  may 
discuss  circle  packing  on  any  Riemann  surface.  In  fact,  it  is  known  that  given  any 
reasonable  pattern  of  tangencies  there  exists  a  unique  Riemann  surface  which  sup¬ 
ports  a  circle  packing  having  that  pattern  of  tangencies  [4],  Not  all  surfaces,  though, 
support  a  circle  packing.  It  has  been  shown  by  Brooks  [12]  and  by  Bowers  and 
Stephenson  [10,  8],  however,  that  these  packable  Riemann  surfaces  are  dense  in  Te- 
ichmiiller  space.  To  go  beyond  this  notion  of  density  and  enter  into  a  deterministic 
discussion  of  how  the  pattern  of  tangencies  affects  the  geometry  of  the  surfaces  and 
their  circle  packings  can  be  very  difficult  [4,  8,  10,  42], 

In  Chapter  III,  we  outline  the  essential  background  material  in  the  theory  of  circle 
packing.  We  then  prove  two  results  regarding  circle  packings  and  surfaces.  First,  we 
prove  a  proposition  of  independent  interest  regarding  the  interaction  between  the 
pattern  of  tangencies  and  the  geometry  of  a  surface  and  the  unique  packing  on  that 
surface.  We  then  prove  a  lemma  regarding  the  density  of  a  subclass  of  packable 
surfaces;  specifically,  we  show  that  we  can  approximate  any  Riemann  surface  with 
a  sequence  of  packable  surfaces  such  that  the  radii  of  the  circles  in  the  packings 
admitted  on  the  surfaces  tend  toward  zero. 

In  Chapter  IV,  we  turn  to  our  main  result,  using  techniques  developed  by  Williams 
[53,  54]  to  describe  a  discrete,  conformal  (combinatorial  or  circle  packing)  version 
of  the  earthquakes  developed  by  Thurston  and  McMullen.  This  process  involves 
opening  a  packing  along  the  pattern  of  tangencies,  shearing  or  grafting  within  the 
combinatorics  and  repacking  to  obtain  a  new  surface.  We  then  prove  that  if  we  have 
a  compact  hyperbolic  Riemann  surface  R  and  a  sequence  of  packable  surfaces  Rk 
converging  to  R  with  mesh  going  to  zero,  then  we  may  approximate  the  action  of 
a  complex  earthquake  on  R  through  combinatorial  earthquakes  on  the  sequence  Rk- 
Further,  we  show  that  this  sequence  of  approximations  is  convergent. 

Having  constructed  approximation  methods  for  earthquakes  on  hyperbolic  Rie¬ 
mann  surfaces,  we  note  that  the  literature  contains  no  reference  to  earthquake  actions 
on  Euclidean  surfaces,  such  as  tori.  In  Chapter  V,  we  extend  the  work  of  Thurston 
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and  McMullen  by  developing  a  theory  of  earthquakes  on  compact  tori.  Once  we  have 
these  earthquake  actions  defined  and  we  demonstrate  that  the  action  thus  described 
is  a  transformation  on  the  Teichmiiller  space  of  tori,  we  develop  discrete  conformal 
(combinatorial  or  circle  packing)  approximations  for  these  Euclidean  earthquakes  sim¬ 
ilar  to  those  described  in  Chapter  IV,  and  show  that  these  approximations  are  again 
convergent.  Finally,  we  discuss  earthquakes  on  the  one-point  Teichmiiller  spaces  of 
the  remaining  compact  surfaces,  giving  a  complete  characterization  of  earthquakes 
on  compact  Riemann  surfaces. 

The  views  expressed  in  this  article  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  United  States  Air  Force,  Department  of  Defense,  or 
the  U.S.  Government. 
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CHAPTER  II 

REAL  AND  COMPLEX  EARTHQUAKES 


2.1  Preliminaries  and  Definitions 

We  begin  by  describing  the  action  of  finite  left  earthquakes  on  the  Poincare  unit 
disk.  We  use  a  definition  for  finite  laminations  of  the  unit  disk,  D  =  {2GC:|2;|<1}, 
from  the  excellent  exposition  of  finite  left  earthquakes  by  Gardiner  and  Lakic  [23,  24] 
rather  than  the  original  definition  for  a  more  general  lamination  used  by  Thurston 
and  others  [21,  35,  49,  51,  54].  This  definition  is  used  here  as  it  is  more  accessible  to 
geometric  and  graphical  examples,  and,  as  we  shall  see,  this  more  “accessible”  finite 
definition  suffices  to  demonstrate  the  results  germane  to  this  investigation. 

Definition  2.1.1.  Let  S  be  a  subset  of  the  unit  circle  consisting  of  n  >  4  points.  A 
finite  geodesic  lamination  C  of  the  unit  disk  D  associated  with  S  is  a  collection 
of  n  —  3  disjoint  hyperbolic  lines  (geodesics)  joining  the  points  of  S  so  that  no  line  in 
C  joins  adjacent  points  of  S.  More  generally,  a  geodesic  lamination  of  D  is  a,  not 
necessarily  finite,  collection  of  pairwise  disjoint  (inU>)  hyperbolic  geodesics. 

The  collection  of  disjoint  geodesics  in  a  finite  geodesic  lamination  is  maximal  in  the 
sense  that  we  may  place  at  most  n  —  3  such  geodesics.  Note  that  a  finite  geodesic 
lamination  divides  D  into  n  — 2  disjoint  hyperbolic  triangles.  Given  such  a  lamination, 
we  may  place  a  weight  (measure)  on  each  geodesic. 

Definition  2.1.2.  For  a  finite  set  Sc®  with  n  >  4  points  and  a  lamination 
C  =  {Tj}”!"]3  associated  with  S  consisting  of  n  —  3  disjoint  hyperbolic  geodesics,  Li,  let 
a  be  a  non-negative  measure  on  C.  That  is,  o  :  Li  i-»  G  [0,  oo),  i  =  1,  2, . . . ,  n  —  3. 
We  call  the  2-tuple  (£,  a)  a  finite  measured  geodesic  lamination  on  D. 

In  Figure  2.1  we  see  three  examples  of  finite  laminations  of  D.  These  are  clearly 
not  the  only  finite  geodesic  laminations  on  the  given  number  of  points  in  each  case.  So, 
given  that  we  can  construct  these  laminations,  we  ask  the  natural  question.  “How 
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(a)  Lamination  on  4  Points 


(b)  Lamination  on  6  Points 


(c)  Lamination  on  16  Points 


Figure  2.1:  Examples  of  Finite  Geodesic  Laminations 


many  distinct  laminations  are  there  on  a  given  set  of  points?”  This  question  is 
essentially  a  restatement  of  Euler’s  polygon  division  problem,  Proposition  2.1.1.  In 
placing  the  n  points  of  S  on  the  unit  circle,  we  divide  it  into  n  subarcs;  that  is,  we 
induce  an  n-sided  polygon.  Then,  in  constructing  a  lamination  through  the  inclusion 
of  n  —  3  geodesics  connecting  points  of  S ,  we  divide  this  polygon  into  three-sided 
figures  or  generalized  triangles. 


Proposition  2.1.1  (Euler’s  Polygon  Division).  A  plane  convex  polygon  ofn>  4 
sides  can  be  divided  into  triangles  by  diagonals  in 

1 


Cn-2  = 


n 


2n  —  4 
n  —  2 


(2.1.1) 


distinct  ways,  where  Cn- 2  is  the  Catalan  number  on  n  —  2. 


This  answer  to  the  question  of  how  many  distinct  finite  laminations  are  possible 
on  n  points  has  been  noted  by  Gardiner,  Hu,  and  Lakic  [23,  24].  See  Dorrie  [19]  for 
the  classic  elementary  proof  of  this  result.  The  laminations  thus  generated  will  fall 
into  classes  of  laminations  that  are  unique  up  to  rotations.  In  Figure  2.2  we  see  the 
topological  classes  of  laminations  possible  on  n  =  7  vertices.  Note  that  there  are  six 
distinct  classes  of  lamination  and  seven  laminations  possible  in  each  class,  obtained 
by  rotations  of  a  representative  of  the  class.  Thus,  we  have  C5  =  42  total  possible 
finite  laminations  on  those  seven  vertices. 


Figure  2.2:  Topological  Types  of  Finite  Laminations  on  n  =  7  Vertices 

Definition  2.1.3.  Given  a  set  S  C  dD  with  n  >  4  points  and  a  lamination  C,  the 
dual  graph  for  the  lamination  is  called  the  tree  of  triangles. 


Figure  2.3:  Finite  Lamination  with  Dual  Graph/Tree  of  Triangles 

In  Figure  2.3  we  see  an  example  of  a  finite  geodesic  lamination  and  the  correspond¬ 
ing  dual  graph  or  tree  of  triangles.  To  call  this  dual  graph  a  tree  makes  combinatorial 
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sense,  since  the  standard  definition  of  a  combinatorial  tree  from  graph  theory  is  a  con¬ 
nected  graph  that  contains  no  subgraph  isomorphic  to  a  cycle  [26] .  The  constraints 
placed  upon  the  definition  of  a  finite  lamination  guarantee  that  the  dual  graph  will 
satisfy  this  definition. 

Definition  2.1.4.  Let  C  be  a  finite  geodesic  lamination  of  D  on  a  set  S  containing 
n  >  4  points.  If  a  triangle  T  induced  by  this  lamination  has  two  sides  that  are  arcs 
of  <90,  then  T  is  called  a  border  triangle. 

Again,  this  definition  of  a  border  triangle  comes  from  Gardiner  and  Lakic  [23,  24]. 
We  have  a  useful,  and  necessary,  result  regarding  these  border  triangles. 

Proposition  2.1.2.  If  L  is  a  finite  geodesic  lamination  of  D  on  a  set  S  containing 
n  >  4  points,  then  C,  contains  at  least  two  border  triangles. 


2.2  Shearing  Maps  and  Real  Earthquakes 
The  following  proposition  is  frequently  used  or  implied  in  work  on  finite  earth¬ 
quakes  since  it  allows  an  explicit  description  of  a  disk  automorphism  used  to  construct 
those  earthquakes  [23,  24].  Due  to  the  relatively  elementary  nature  of  the  proposition, 
however,  no  proofs  appear  in  the  literature.  In  the  interest  of  a  complete  exploration 
of  the  properties  and  construction  of  finite  left  earthquakes,  we  prove  it  here. 

Proposition  2.2.1.  Let  L  C  ID)  be  a  hyperbolic  geodesic  with  endpoints  a,b  G  dD,  a  ^ 
b;  call  a  and  b  the  left  and  right  endpoints,  respectively,  of  L.  Let  p  >  0.  The  map 

(aeM  —  b)z  +  ab(  1  —  eM) 


C(z) 


(2.2.1) 


(e^  —  1  )z  +  (a  —  be11) 

is  a  disk  automorphism  with  exactly  two  fixed  points  at  a,  b  G  <9D  such  that  C  (L)  =  L 
and  points  on  L  are  translated  a  fixed  hyperbolic  distance  p  to  the  left  along  the 
geodesic  (away  from  b  and  toward  a). 

Proof.  To  show  that  C(z )  has  exactly  two  fixed  points  at  a,  b  G  <9ID)  is  a  simple  exercise 
in  algebra;  we  simply  solve  the  equation  C(z )  =  z  for  z.  This  immediately  gives  that 
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C(z )  has  exactly  two  fixed  points  at  a,  &  £  c®.  To  prove  that  the  transformation 
C(z )  is  a  disk  automorphism,  we  recall  [38]  that  a  Mobius  transformation 


M(z) 


u>i  z  T  U2 

UJ3Z  +  CO4 


(2.2.2) 


where  U1U4  —  cj2uj3  /  0,  is  a  disk  automorphism  if  and  only  if  there  exist  constants 
a  £  D  and  9  €  E  such  that 

M(z)  =  (2.2.3) 

Since  ( ae P  —  b)(a  —  be M)  —  (e*1  —  l)(a&(l  —  eM))  =  e^(a  —  b)2  =  0  if  and  only  if  a  =  b, 
C(z)  is  indeed  a  Mobius  transformation.  We  can  express  C(z)  as 


(ae^  —  b)z  +  ab(l  —  eM)  _  /  aeM  —  b\  iz~  \ 

(ef1  -  1  )z  +  (a  -  be11)  \a-bef1)  J  ' 


Therefore,  in  order  to  show  that  C(z)  is  a  disk  automorphism  we  need  only  verify 
that  the  right-hand  side  of  Equation  (2.2.4)  satisfies  the  necessary  criteria. 

ae 9  —  b_  I  (aei*  —  b)  (aeJ1  —  b) 
a  — be11  y  (a  —  bet1)  (a  —  be t1) 


e 2t*  —  2e^R  e(ab)  +  1 
e2!1  —  2e>lRe(ab)  +  1 


ablet1  —  1) 
aet1  —  b 


e*  -  1 


\aeti  —  b\  yj  e2^  —  2e^Re(a6)  +  1  s/e2*1  —  2e^  +  1 


=  1  (2.2.6) 


/  ab^et1  —  1)  \  _  ab^t 1  —  1)  _  aa(l  —  eM)  _  1  —  _  1  —  .  . 

\  aet1  —  b  J  aeM  _  b  a  —  a  —  -e*1  a  —  bet1 

Thus,  Clz )  is  a  disk  automorphism  preserving  the  unit  circle.  Now,  since  C(z)  is 

a  Mobius  transformation,  the  image  of  a  circle  is  a  circle.  Further,  since  Clz)  is 

analytic  except  for  a  single  simple  pole  in  C\B,  Clz)  is  conformal  in  a  neighborhood 

of  each  fixed  point.  Thus,  any  circle  T  passing  through  a  and  b  must  have  as  its 

image  under  Clz)  a  circle  through  a  and  b  meeting  the  unit  circle  in  the  same  angle 

as  T.  This  is  sufficient  to  ensure  that  C'(T)  =  T  for  any  circle  T  through  a  and  b ; 

in  particular,  C'(L)  =  L.  Finally,  to  show  that  points  on  L  are  translated  a  fixed 

hyperbolic  distance  n  to  the  left  along  the  geodesic  (away  from  b  and  toward  a)  recall 


[2]  that  the  translation  length  p(z,M(z ))  of  a  hyperbolic  Mobius  transformation  of 
the  form  given  in  Equation  2.2.2  is  given  by  the  relation 


\  \tr{M)\  =  cosh  C^p(z,M(z)) 


(2.2.8) 


where 


[tr{M)f  = 


2  _  (W1  +  w4)2 


—  UJ2U3 


Applying  this  relation  to  C(z )  we  have 

c°sh  (^p{z,M(z))^j  = 


efl  + 1 
V2e^' 


and  this  reduces  to 


p(z,  M(z))  =  2  cosh  1  (^==f  )  =  ln  (e#1)  =  P- 


This  completes  the  proof. 


(2.2.9) 


(2.2.10) 


(2.2.11) 

□ 


Now,  given  a  finite  measured  geodesic  lamination  (£,  a)  on  D,  we  use  the  hyper¬ 
bolic  Mobius  transformation  described  in  Proposition  2.2.1  to  construct  a  discontin¬ 
uous  map  ha  :  B  — >  D  (which  extends  to  a  continuous  map  on  dH>)  with  the  following 
properties. 


1.  ha( D)  =  O. 

2.  h  :  S  — >  (9D  defined  by  h  =  ha\s  is  a  cyclic  and  order-preserving. 

3.  On  the  interior  of  any  triangle  induced  by  C,  ha  is  a  hyperbolic  isometry. 

4.  If  L  €  C  is  the  common  hyperbolic  geodesic  between  neighboring  triangles  Tx 
and  T2,  then  the  points  of  T2  move  to  the  left  relative  to  the  points  of  T\. 


We  will  call  the  map  ha  a  finite  left  earthquake  on  D.  Note  that  if  we  fix  the 
orientation  of  the  geodesics  as  given  in  the  finite  lamination  C  (i.e.,  we  continue  to 
consider  left  and  right  endpoints  as  above)  and  globally  replace  the  positive  weights 
in  the  measured  geodesic  lamination  with  negative  weights,  we  obtain  a  finite  right 
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earthquake.  All  of  the  results  and  characteristics  pertaining  to  left  earthquakes 
also  hold  for  these  right  earthquakes.  In  addition,  if  we  relax  the  restriction  that  all 
weights  be  either  positive  (to  induce  a  finite  left  earthquake)  or  negative  (to  induce 
a  finite  right  earthquake)  and  allow  both  positive  and  negative  weights,  we  have 
legitimate  finite  earthquake  actions  at  the  expense  of  uniqueness  results  to  be  given 
in  Proposition  2.2.2.  This  more  general  class  of  finite  “mixed”  earthquakes  will  be  the 
class  we  consider  in  constructing  discrete  conformal  approximations  (circle  packing 
approximations)  for  the  earthquake  actions  in  Chapter  IV. 

We  now  explicitly  construct  the  map  ha.  Let  T\  be  any  triangle  induced  by  the 
finite  measured  geodesic  lamination  (C,a).  At  least  one  side  of  T\  must  be  a  geodesic 
L\  in  C  with  left  and  right  endpoints  a\  and  ,  respectively,  as  seen  from  the  interior  of 
Ti,  and  having  weight  cr(Li)  —  /q.  Using  these  parameters,  construct  the  hyperbolic 
Mobius  transformation  6\,  called  the  comparison  map,  defined  in  (2.2.1),  and  apply 
Ci  “down  the  tree  of  triangles.”  That  is,  apply  C\  to  all  edges  (geodesics)  and  nodes 
(triangles)  below  in  the  tree.  In  a  sense,  we  apply  the  discontinuous,  piecewise 
function 

{C\  (z)  if  z is  “below”  Ti  on  the  tree  of  triangles 

(2.2.12) 

z  otherwise 
to  the  closed  unit  disk. 

Now,  choose  any  triangle  adjacent  to  T\.  (If  T  has  more  than  one  neighbor  in  the 
tree  of  triangles,  choose  one,  follow  the  construction  down  that  branch  of  the  tree, 
return  to  Ti,  and  continue  down  any  other  branches  associated  with  those  additional 
neighbors.)  Call  this  neighbor  T2.  If  T2  is  a  border  triangle,  then  the  earthquake 
action  is  complete  on  this  branch  of  the  tree  of  triangles.  If  T2  is  not  a  border 
triangle,  then  one  side  of  T2  is  a  geodesic  L2  not  equal  to  Lx.  Note  that  this  is  not 
necessarily  a  geodesic  in  the  lamination  £;  it  may,  rather,  be  the  image  of  one  such 
geodesic  under  the  action  of  C\.  Since  Ci  is  a  Mobius  transformation,  though,  we  are 
assured  that  L2  is,  in  fact,  a  geodesic,  and  we  may  repeat  the  construction  above.  As 
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with  Lx,  L2  has  left  and  right  endpoints  a2  and  b2,  as  seen  from  the  interior  of  T2,  and 
weight  <t(L2)  =  p,2;  the  weights  on  our  geodesics  are  not  changed  by  the  application 
of  prior  comparison  maps.  Note  that  exactly  one  of  a2  and  b2  is  equal  to  exactly  one 
of  a\  and  bi,  while  the  other  of  a2  and  b2  is  equal  to  either  Ci(ai)  or  Cj(&i).  As 
before,  construct  the  hyperbolic  disk  automorphism  C2  defined  in  (2.2.1)  from  these 
parameters,  and  apply  C2  down  the  tree  of  triangles.  Continue  in  this  manner  down 
every  branch  of  the  tree  of  triangles  until  reaching  a  border  triangle,  at  each  step 
constructing  comparison  maps  from  the  images  points  in  S  under  applications  of  the 
preceding  comparison  maps.  Thus,  if  there  are  n  points  in  S,  there  will  be  n  —  3 
comparison  maps,  one  for  each  geodesic  in  the  lamination  £. 

When  the  tree  of  triangles  has  been  exhausted,  we  have  a  map  hff  from  the  closed 
unit  disk  to  the  closed  unit  disk  such  that  the  restriction  of  hc  to  <90  is  a  cyclic, 
order-preserving  homeomorphism  from  the  unit  circle  to  itself.  Clearly  the  map  h„  is 
discontinuous  on  the  interior  of  O,  but  our  main  concern  will  be  the  action  of  h„  on 
the  90. 

In  Example  2.2.1  we  construct  an  explicit  finite  left  earthquake  from  a  given  finite 
measured  geodesic  lamination. 


Example  2.2.1.  Let  S  =  {si,  s2,  s3,  s4,  s5}  =  |l,el?,z, e*T,  — 1  j.  Let  £  =  {Li,L2}, 
and  let  Ti,  T2,  and  T3  be  as  shown  in  Figure  2.4a.  Define  a  by  cr(Li)  =  ln(2)  and 
ct(L2)  =  ln(3).  These  values  for  the  weights  on  Li  and  L2  may  appear  arbitrary,  but 
they  were  chosen  for  reasons  of  computational  convenience.  Notice  that  in  Figure 
2.4a  an  additional  dashed  circle  is  drawn,  and  we  have  included  two  additional  points 
on  the  geodesics  Li  and  L2.  These  objects  have  no  direct  bearing  on  the  action  of 
the  finite  left  earthquake;  they  do  serve,  however,  as  a  reference  and  to  illustrate  that 
action.  Using  the  triangle  labeled  7\  as  the  “top”  of  the  tree  of  triangles,  we  construct 
the  map 


Ci(z)  = 


(2  —  i)z  —  i 


(2.2.13) 


z  +  (1  —  2  i) 

from  the  left  and  right  endpoints  and  the  weight  associated  with  L\.  We  apply  this 
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disk  automorphism  in  the  piecewise  manner  described  above;  the  result  is  shown  in 
Figure  2.4b.  Note  that  the  shaded  portion  of  the  closed  disk  is  unchanged,  while 
points  in  the  unshaded  portion  have  been  translated  to  the  left  as  seen  from  T\. 


(a)  Lamination  on  ©  (b)  Action  of  C\  ( z ) 


Figure  2.4:  Finite  Left  Earthquake 


Now,  we  construct  a  second  map  using  the  new  endpoints  of  the  geodesic  L2.  The 
explicit  map  thus  constructed  is  given  by 

cm*)  =  (2-2'14) 

We  apply  this  transformation  down  the  tree  of  triangles  to  obtain  the  result  shown 
in  Figure  2.4c.  Again,  the  shaded  portion  of  the  closed  disk  is  unaffected  by  the 
piecewise  application  of  the  comparison  map  C2,  and  points  in  the  unshaded  region 
have  been  translated  to  the  left  relative  to  T2. 

In  Example  2.2.1  we  illustrated  the  action  of  a  simple  finite  left  earthquake.  In 
doing  so,  however,  we  chose  a  triangle  from  the  lamination  to  act  as  the  “first” 
in  the  tree  of  triangles  induced  by  the  lamination.  Thereafter,  the  action  of  the 
earthquake  was  determined,  but  the  requirement  that  we  choose  a  starting  point  for 
the  execution  of  a  finite  left  earthquake  raises  the  question  of  uniqueness.  Would  a 
different  choice  for  the  initial  triangle  result  in  a  different  earthquake?  The  answer  is 
that  the  fundamental  action  of  a  finite  left  earthquake  is  independent  of  the  choice 
of  initial  triangle.  This  result  and  an  explanation  of  what  is  meant  by  “fundamental 
action”  is  given  in  Proposition  2.2.2. 
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Proposition  2.2.2.  The  effect  of  a  finite  left  earthquake  on  S  C  c®  induced  by 
a  finite  measured  lamination  (£,  a)  associated  with  S  is  uniquely  determined,  up  to 
post-composition  by  a  disk  automorphism,  by  the  measured  lamination. 


For  a  proof  of  this  result  see  the  work  of  Gardiner,  Hu,  and  Lakic  [23,  24].  We 
do  not  include  their  proofs  here  since  a  concrete  example  of  this  property  is  perfectly 
illustrative  and  demonstrates  beautifully  the  essential  mechanism  of  their  proofs,  the 
explicit  formulation  of  an  appropriate  Mobius  transformation.  As  such,  we  rather 
include  Example  2.2.2  as  a  illustration  of  the  essentially  unique  determination  of  the 
action  of  a  finite  left  earthquake  on  the  points  of  S. 


Example  2.2.2.  Suppose  that  S  =  je  (ie1)’  j  ,  and  let  C  be  the  finite  geodesic 
lamination  associated  with  S  as  shown  in  Figure  2.5a.  We  have  numbered  the  fourteen 


triangles  induced  by  the  lamination  as  Ti,T2,  . . .  ,Ti4  for  reference,  without  necessarily 


implying  any  a  priori  decisions  regarding  the  manner  in  which  the  finite  left  earth¬ 
quake^)  are  executed.  Place  a  uniform  weight  of  p  =  In  (|)  on  each  geodesic  of  the 
lamination. 


(a)  Lamination  on  B  (b)  Start  with  Ti  (c)  Start  with  T4 


Figure  2.5:  “Uniqueness”  of  Finite  Left  Earthquakes 


In  Figures  2.5b  and  2.5c  we  show  the  action  of  the  finite  left  earthquakes  induced 
by  placing  Ti  and  T4  at  the  top  of  the  tree  of  triangles,  respectively.  The  action 
of  the  earthquake  appears  very  different  in  each  case.  As  Gardiner,  Hu,  and  Lakic 
[23,  24]  point  out,  however,  if  we  construct  a  Mobius  transformation  by  composing 
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the  comparison  maps  associated  with  the  unique  sequence  of  geodesics  separating  the 
two  triangles  Ti  and  T4  as  simple  disk  automorphisms  rather  than  as  piecewise  func¬ 
tions,  we  obtain  immediately  the  desired  Mobius  transformation.  In  this  example, the 
desired  maps  are 


Ci(z) 


(.79  -  .ni)z  -  (.35  +  .35*) 
.5z  —  (.06  +  1.06?) 


C,(z)  ~  (1-46  —  i)z  —  (.02  +  .5*) 
°2{Z)~  .5z  +  (.94  -1.5*) 

(2.49  +  .16*)z  +  (-49  +  .08*) 

^3(z)  ~  " 


(2.2.15) 

(2.2.16) 
(2.2.17) 


,5z  +  (2.48  +  .24*) 

Note  that  the  Mobius  transformations  shown  in  (2.2.15),  (2.2.16),  and  (2.2.17)  are 
each  created  based  on  the  images  of  their  endpoints  under  the  previous  comparison 
maps  applied  in  the  application  of  the  finite  left  earthquake  induced  by  placing  7\  at 
the  top  of  the  tree  of  triangles.  We  do  not  use  the  endpoints  of  those  geodesics  in  the 
original  lamination,  though  any  geodesic  that  forms  a  side  of  the  initial  triangle  in 
the  tree  is  always  fixed  by  the  action  of  a  finite  left  earthquake,  and  we  will  thus  use 
the  initial  values  for  the  endpoints  of  those  geodesics.  Composing  these  three  Mobius 
transformations  we  have 


M(z)  =  (C3oC2oC1)(z)  « 


(1.95  -  5.57 i)z  -  (4.28  +  1.19*) 
(2.64  -  3.57 i)z  -  (4.97  +  3.19*) ' 


(2.2.18) 


A  few  routine  calculations  show  that  the  image  points  of  S  under  the  two  earthquake 
maps  shown  in  Figure  2.5  differ  only  by  post-composing  the  disk  automorphism  M(z) 
given  in  (2.2.18)  with  the  finite  left  earthquake  action  resulting  from  placing  T\  at 
the  top  of  the  tree  of  triangles.  Note  that  should  we  wish  to  alter  the  order  here, 
that  is  begin  with  the  earthquake  induced  by  placing  T4  at  the  top  of  the  tree  and 
post-compose  by  a  Mobius  transformation,  then  we  do  not  simply  reorder  the  maps 
Ci,  C<2  and  C3;  rather,  we  must  calculate  a  completely  new  set  of  maps  based  on  the 
images  of  geodesics  in  that  finite  left  earthquake. 


Thus  we  have  seen  that  given  a  finite  measured  geodesic  lamination  on  D,  we 
can  explicitly  construct  and  describe  the  action  of  the  (essentially)  unique  finite  left 
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earthquake  induced  by  that  lamination.  The  next  natural  step  is  to  ask  whether  or 
not  given  the  action  of  a  finite  left  earthquake  we  can  determine  the  lamination  which 
induced  it.  That  is,  given  a  finite  set  of  points  S  C  83  and  their  images  under  a 
finite  left  earthquake,  can  we  find  a  finite  measured  geodesic  lamination  (£,  a)  that 
induces  the  given  map  on  S ?  The  answer  was  given  by  Thurston  [49]. 

Theorem  2.2.1  (Thurston’s  Finite  Earthquake  Theorem).  Suppose  that  h  is  a 
given  order-preserving  map  from  a  finite  subset  S  of  the  unit  circle  into  the  unit  circle. 
Then  there  exists  a  unique  finite  measured  geodesic  lamination  (£,  a)  associated  with 
S  such  that,  up  to  post- composition  by  a  Mobius  transformation,  h  is  the  restriction 
to  S  of  the  finite  left  earthquake  induced  by  (C,a). 

Gardiner  and  Lakic  [24]  and  Gardiner,  Hu,  and  Lakic  [23]  give  excellent  proofs  of 
this  theorem.  Their  proofs  are  of  particular  value  and  interest  in  that  they  are  espe¬ 
cially  constructive;  their  proofs  suggest,  inductively,  a  method  to  explicitly  calculate 
the  desired  (unique)  finite  measured  geodesic  lamination  given  the  sets  S  and  h(S). 
We  do  not  describe  the  algorithm  here,  but  we  have  implemented  the  algorithm  in 
Mathematica,  and  the  code  used  is  provided  in  Appendix  C. 

Example  2.2.3.  We  take  a  “random”  set  of  five  points  on  the  unit  circle  and  another 
set  of  “random”  points  to  act  as  their  images  under  h.  To  say  that  the  points  used 
are  random  is  somewhat  misleading.  In  fact,  the  points  in  this  example  were  chosen 
randomly  within  the  constraints  that  we  require  that  h  is  cyclic  and  order-preserving. 
The  points  used  in  this  example  are  shown  (to  six  significant  figures)  in  Table  2.1.  Us¬ 
ing  the  Mathematica  code  given  in  Appendix  C,  we  find  the  two  appropriate  geodesics 
and  weights  on  each  as  shown  in  Figure  2.6. 

Observe  that  any  finite  left  earthquake  on  this  finite  measured  geodesic  lamination 
fixes  the  three  points  in  the  triangle  placed  at  the  top  of  the  tree  of  triangles.  Thus, 
the  earthquake  induced  by  the  calculated  lamination  does  not  give  the  correct  images 
in  h(s );  in  fact,  no  finite  left  earthquake  can  do  so  alone.  The  three  points  fixed  by 
the  earthquake  point  out  a  choice  for  the  Mobius  transformation  which  will  map  the 
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Table  2.1:  The  Points  S  and  h(S)  Used  to  Compute  a  Lamination 


s 

m 

Si 

0.259205  +  0.965822* 

h(si) 

-0.485661  +  0.874147* 

S2 

-0.294304  +  0.955712* 

h(s2) 

-0.932807  +  0.360375* 

s3 

-0.410866  -  0.911696* 

h(s3) 

-0.171589  -  0.985169* 

s4 

0.733058  -  0.680167* 

/*(s4) 

0.8436120  -  0.536941* 

s5 

0.887984  -  0.459873* 

h{ss) 

0.987219  -  0.159371* 

Figure  2.6:  Measured  Geodesic  Lamination  on  S  Given  h(s) 

points  as  desired.  We  calculate  the  unique  Mobius  transformation  sending  the  three 
points  in  the  first  triangle  to  their  images;  in  fact,  the  algorithm  used  to  compute  the 
lamination  makes  this  the  necessary  choice  of  Mobius  transformation.  Thus,  in  this 
case  we  have  three  possible  finite  left  earthquakes  and  three  corresponding  Mobius 
transformations  given  in  (2.2.19),  (2.2.20),  and  (2.2.21),  corresponding  to  the  finite 
left  earthquakes  induced  by  placing  the  triangles  Ti,  T2,  and  T3,  respectively,  at  the 
top  of  the  tree  of  triangles,  as  shown  in  Figure  2.6. 


Mx{z) 


-(1.470  +  1.478*)*  +  (0.0489  +  0.496*) 
(0.472  -  1.580*)*  -  (1-950  +  0.738i) 


(2.2.19) 
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(2.2.20) 


M2(z)  = 


(2.287  +  3.370 i)z  -  (0.211  +  0.611*) 
(-0.628  +  0.152*)^  +  (3.571  +  1.957*) 


M3(z) 


(-0.807  +  0.640*)z  +  (0.107  -  0.189*) 
(0.173  -  0.132 i)z  -  (0.524  -  0.886*) 


(2.2.21) 


2.3  Grafting  Maps  and  Complex  Earthquakes 
In  geologic  terms,  earthquakes  do  not  only  involve  lateral  “shearing”  along  a  fault 
line.  Earthquakes  may  also  involve  a  separation  in  which  two  tectonic  plates  move 
apart  and  new  surface  rises  to  fill  the  gap  or  a  subduction  in  which  one  tectonic  plate 
slides  beneath  another,  and  surface  is  lost.  Just  as  we  have  maps  which  mimic  a 
geologic  shearing  action,  so  we  have  maps  which  mimic  these  separation  and  subduc¬ 
tion  actions.  We  will  consider  these  together  as  “grafting”  actions  which  are  positive 
or  negative,  respectively,  in  their  behavior.  Shearing  and  grafting  maps,  combined 
through  composition,  comprise  what  we  will  call  complex  earthquakes. 

In  this  section  we  first  revisit  the  shearing  actions  developed  in  Section  2.2,  giving 
an  alternate  description  through  their  local  projective  behavior.  We  then  use  this 
local  projective  structure  to  describe  a  new  form  of  earthquake,  a  grafting  action 
rather  than  a  shearing  action,  introduced  by  McMullen  [35],  from  whom  much  of  the 
discussion  below  originates,  and  investigated  by  others  [40,  41]. 


Figure  2.7:  Local  Projective  Description  of  a  Shearing  Map 


First,  we  revisit  the  action  of  a  shearing  map.  Suppose  a,b  G  <90  such  that 
a  /  5,  and  L  C  D  is  the  unique  hyperbolic  geodesic  connecting  a  to  b.  Apply  the 
Mobius  transformation  M(z)  which  takes  the  unit  disk  D  to  the  upper  half  plane  H 
so  that  the  image  of  L  is  the  positive  imaginary  axis,  the  image  of  50  is  the  real  axis, 
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M(a)  =  co,  and  M(b)  =  0.  Now,  we  apply  a  simple  scaling,  rz,  where  r  >  0,  in  one 
quadrant  of  H  and  the  identity  map  on  the  other.  Finally,  we  apply  the  inverse  of  the 
Mobius  transformation  M.  As  we  can  see  in  Figure  2.7,  the  effect  of  thus  composing 
these  maps  is  that  of  a  hyperbolic  shearing.  In  fact,  should  we  perform  this  general 
composition  and  simplify  the  resulting  function,  we  will  find  that  the  composition 
map  is  the  comparison  map  given  in  (2.2.1). 

Before  we  turn  ourselves  to  the  description  of  grafting  maps  and  earthquakes,  we 
need  to  recall  some  properties  of  the  upper  half  plane  H  as  a  model  for  the  hyperbolic 
plane  [2].  First,  recall  that  the  density  we  use  to  define  a  hyperbolic  metric  on 
H  =  {z  G  C  :  Im(;z)  >  0}  is  given  by 


1 


Im(z)  ‘ 


(2.3.1) 


As  a  consequence,  what  we  think  of  as  a  strip  in  a  Euclidean  sense  may  not  be  a  strip 
in  HL  Similarly,  a  strip  in  the  upper  half  plane  may  not  appear,  to  our  Euclidean  eye, 
as  a  strip  at  all.  Specifically,  we  note  that  a  Euclidean  wedge  in  H  with  its  corner 
touching  the  real  line,  as  shown  in  Figure  2.8,  is  a  strip  in  the  hyperbolic  metric  on 
Et,  since  the  distance  from  edge  to  edge  of  the  wedge  is  constant  in  the  hyperbolic 
metric.  The  circular  arcs  shown  in  Figure  2.8  are  all  of  equal  length  in  the  hyperbolic 
metric  on  EL 


Im(z) 


Re(z) 


Figure  2.8:  Hyperbolic  Strip  in  H 


To  define  the  grafting  maps  and  imaginary  earthquakes  described  by  McMullen 
[35],  we  use  the  approach  thus  described.  The  difference  is  that  rather  than  a  piecewise 
scaling  action  given  as  multiplication  by  a  real  number,  we  apply  a  piecewise  rotation 
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by  some  angle  a  given  as  multiplication  by  the  complex  number  eia.  This  piecewise 
rotation  leaves  a  gap;  we  fill  this  gap  by  grafting  in  a  hyperbolic  cylinder,  as  shown 
in  Figure  2.9c.  Now,  we  can  map  the  image  of  the  upper  half  plane  under  the  maps 
so  far  applied  back  to  the  upper  half  plane  with  a  power  map  z&,  where 


7 r 


tt  +  a 


(2.3.2) 


This  map  restores  the  real  axis,  and  allows  us  to  directly  map  the  image  region  back 
to  the  unit  disk  and  see  the  action  of  the  earthquake  on  the  disk.  This  is  illustrated 
in  Figure  2.9. 


Figure  2.9:  Local  Projective  Description  of  a  Grafting  Map 


Note  that  in  the  application  of  the  power  map  z&,  the  arcs  of  circles  and  geodesics 
in  Figure  2.9c  are  taken  to  curves  which  are  no  longer  either  geodesics  or  circles  in 
Figure  2.9d.  This  is  especially  clear  as  we  return  to  the  unit  disk  B  in  Figure  2.9e. 
The  two  curves  connecting  the  points  a,  b  G  <90  are  clearly  not  geodesics,  and  the 
two  disjoint  halves  of  the  dashed  reference  circle  shown  in  Figure  2.9a  are  no  longer 
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circles.  This  failure  to  preserve  geodesics  and  circles  is  key.  Given  a  finite  lamination 
of  the  disk,  we  may  execute  shearing  maps  on  those  geodesics  sequentially  since  the 
maps  involved  are  hyperbolic  isometries  and  thus  preserve  the  geodesics.  In  the  case 
of  grafting  maps,  however,  we  cannot  perform  the  grafting  actions  in  such  a  sequence, 
since  the  first  application  transforms  the  other  geodesics  of  the  lamination  into  curves 
which  are  no  longer  geodesics.  To  make  the  action  of  these  grafting  maps  computable 
is  the  major  motivation  of  this  research. 

Now  we  simply  define  a  complex  earthquake  as  the  composition  of  shearing  and 
grafting  maps,  applied  in  that  order.  Consider  the  shearing  and  grafting  maps  devel¬ 
oped  in  Figure  2.7  and  Figure  2.9.  This  composition  is  illustrated  in  Figure2.10. 


(a)  Geodesic  (b)  Project  to  H  (c)  “Scale”  in  H 


(d)  “Rotate”  in  H  (e)  Power  Map  (f)  Return  to  H3> 

Figure  2.10:  Local  Projective  Description  of  a  Complex  Earthquake  Map 


To  call  such  a  map  a  “complex”  (as  opposed  to  real  or  imaginary)  earthquakes  is 
quite  natural.  A  real  weight  indicated  a  scaling,  multiplication  by  a  real  number,  and 
an  imaginary  weight  indicated  a  rotation,  multiplication  by  an  imaginary  number 
with  unit  modulus.  Now,  a  complex  weight  x  +  iy  on  the  geodesic  indicates  that  the 
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transformation  applied  to  points  in  one  quadrant  of  H  amounts  to  multiplication  by  a 
complex  number  ex+w  —  exeiy .  From  this  point,  when  we  refer  to  measured  geodesic 
laminations,  we  allow  the  measure  to  take  complex  values.  We  will,  in  the  case  of 
earthquakes  on  compact  hyperbolic  surfaces,  restrict  the  imaginary  part  of  complex 
measures  to  non-negative  values.  Otherwise,  we  place  no  a  priori  restriction  on  the 
weights  in  a  measured  geodesic  lamination. 

2.4  Riemann  Surfaces  and  Teichmiiller  Theory 
Our  eventual  goal  is  to  define,  discuss  the  properties  of,  and  approximate  the 
action  of  earthquakes,  both  grafting  and  shearing  maps,  on  the  Teichmiiller  spaces  of 
Riemann  surfaces.  Here  we  review  some  of  the  important  definitions  and  properties 
of  Riemann  surfaces  and  their  Teichmiiller  spaces  [14,  15,  24,  29,  30,  33,  48]. 

Definition  2.4.1.  A  Riemann  surface  is  a  one  complex- dimensional  manifold 
with  charts  whose  overlap  maps  are  analytic.  The  maximal  collection  of  charts  on  a 
Riemann  surface  is  the  conformal  structure  for  that  surface. 

The  conformal  structure  for  Riemann  surface  R  (that  is  not  the  sphere)  gives  a 
convenient  way  of  representing  the  surface  in  the  complex  plane  C,  if  R  is  Euclidean, 
or  in  the  hyperbolic  plane  D,  if  R  is  hyperbolic. 

Definition  2.4.2.  Let  R  be  a  Riemann  surface  with  fundamental  group  tti(R).  A 
closed,  connected  subset  P  —  of  C  (orB)  is  called  a  fundamental  region  for  R  if 

1.  for  each  z  E  C  (or  z  E  B)  there  exists  at  least  one  u  E  P  such  that  z  =  g(co), 
where  g  E  7Ti  (R)  (we  say  z  ~  u);  and 

2.  for  every  Wi,u>2  E  P,  there  does  not  exist  an  element  g  E  it i(R)  such  that 
g(u>i)  =w2. 

A  torus  T  with  fundamental  group  given  by  <  z-\-l,z+{l+i)  >  has  as  a  fundamental 
region  the  one  shown  in  Figure  2.11.  Here  we  have  Zi  ~  w2  for  each  i  =  1,2, 3, 4, 5, 
and  oji  7^  g{w 2)  for  any  g  E<  z  +  1,  z  +  (1  +  i)  >. 
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Im(z) 


Figure  2.11:  Fundamental  Region  for  a  Torus 

Two  Riemann  surfaces  Ri  and  R2  are  said  to  be  conformally  equivalent  if  and 
only  if  there  exists  a  conformal  homeomorphism  /  :  Ri  — »  R2.  The  set  of  equivalence 
classes  of  surfaces  of  the  same  topological  type  as  Ri  under  this  equivalence  relation  is 
the  moduli  space  of  R\.  The  equivalence  relation  defined  by  conformal  equivalence 
which  determines  the  moduli  space  of  Ri  is  not  adequate  for  our  purposes.  We  will 
require  that  for  surfaces  i?i  and  R2  to  be  considered  equivalent  they  must  first  be 
conformally  equivalent  (the  same  point  in  moduli  space)  and  we  will  further  require 
that  the  generators  of  their  fundamental  groups  correspond. 

Definition  2.4.3.  Let  R  be  a  Riemann  surface,  and  let  E  be  a  collection  of  canonical 
generators  for  tti(R).  The  collection  E  is  called  a  marking  for  R.  Two  markings  on 
R  are  equivalent  if  and  only  if  they  differ  by  the  choice  of  their  base  point. 

This  now  gives  us  a  new  equivalence  relation  which  we  use  to  describe  another 
space  of  surfaces. 

Definition  2.4.4.  Two  marked  Riemann  surfaces  ( R ,  E)  and  (R,  E)  are  said  to  be 
equivalent  if  and  only  if  there  exists  a  conformal  map  f  :  R  — >■  R  for  which  the 
marking  /(E)  is  equivalent  to  E.  The  Teichmuller  Space  of  R  is  the  set  of  these 
equivalence  classes. 

An  illustration  of  how  two  surfaces  may  be  conformally  equivalent,  and  thus  equiv¬ 
alent  points  in  moduli  space,  while  differing  in  their  markings,  making  them  different 
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points  in  Teichmiiller  space,  is  given  in  Figure  2.12.  The  tori  (or  fundamental  regions 
and  markings  for  the  tori)  T  and  T  shown  in  Figure  2.12a  and  Figure  2.12b,  respec¬ 
tively,  are  conformally  equivalent  [30,  45].  As  points  in  Teichmiiller  space,  however, 
they  differ  in  the  choice  of  markings.  The  generator  z  +  1  corresponds  to  “going 
through”  the  hole  in  the  torus  once,  and  the  generators  z  +  u>  and  z  +  Cj  each  corre¬ 
spond  to  “going  around”  the  hole  in  the  torus  once.  Each  of  z+uj  and  z+ui,  however, 
correspond  to  closed  loops  on  the  torus  in  different  homotopy  classes,  since  z+co  does 
not  “go  through”  the  hole  and  z  +  u  “goes  through”  the  hole  once,  as  illustrated  in 
Figure  2.12c. 


(c)  Marking  for  z  +  &  on  T 


Figure  2.12:  Moduli  Space  vs.  Teichmiiller  Space 


Another  useful  description  of  the  Teichmiiller  space  of  a  surface  involves  equiva¬ 
lence  classes  of  maps  from  a  reference  or  base  surface.  First  we  require  a  generalization 
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of  the  concept  of  a  conformal  map. 

Definition  2.4.5.  Let  f  be  a  continuous,  orientation-preserving  map  from  a  domain 
f2  C  C  into  the  complex  plane,  and  fix  z0  E  D.  Define  Le  and  te  by 

Le  =  max  {|  f(z)  -  f(zo) |  :  \z  -  z0 \  =  e}  , 

and 

4  =  min{\f(z)  -  f(zo) |  :  \z  -  z0\  =  e} , 
as  shown  in  Figure  2.13.  The  dilatation  of  f  at  z0  is 

Df(z 0)  =  limsup 

£ — >-0 

If  supz€fl  {Df(z)}  <  K,  then  we  say  that  f  :  Q  — >■  C  is  K -quasiconformal. 

f 


Figure  2.13:  Action  of  a  Map  /  with  Dilatation  Df(zo)  >  1 

Before  moving  on,  we  state  as  Lemma  2.4.1  a  well-known  and  useful  property  of 
quasiconformal  maps. 

Lemma  2.4.1.  If  f  :  LI  — >  A  is  K\- quasiconformal  and g  :  A  — >  C  is  K2- quasiconformal, 
then  their  composition  g  o  f  :  ft  — »  C  is  [Kf)  {Kf) -quasiconformal. 

Note  that  the  class  of  quasiconformal  maps  contains  those  maps  which  are  conformal. 
Since  a  conformal  map  sends  infinitesimal  circles  to  infinitesimal  circles  [18],  the 
limit  of  the  ratios  described  in  Definition  2.4.5  is  1.  That  is,  a  conformal  map  is 
1-quasiconformal.  Now,  the  class  of  quasiconformal  maps  allows  us  to  define  a  new 
equivalence  relation  of  the  set  of  Riemann  surfaces. 

Definition  2.4.6.  Quasiconformal  maps  fi  and  f2  defined  on  a  Riemann  surface  R 
are  Teichmuller  equivalent  if  and  only  if  /2  o  /f 1  is  homotopic  to  a  conformal  map. 
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Now,  Proposition  2.4.1  establishes  the  equivalence  of  these  two  characterizations  of 
points  in  Teichmiiller  space. 

Proposition  2.4.1.  Fix  a  Riemann  surface  R,  and  suppose  f\  and  f2  are  maps  from 
R  to  Riemann  surfaces  Ri  and  R2,  respectively.  R\  and  R2  are  equivalent  in  the 
Teichmiiller  space  of  R  if  and  only  if  fi  and  f2  are  Teichmiiller  equivalent. 

Thus,  points  in  the  Teichmiiller  space  of  a  Riemann  surface  R  may  be  considered 
as  either  points  or  as  functions.  There  is  a  natural  metric  on  Teichmiiller  space  as 
a  function  of  how  close  to  conformal  (or  how  quasiconformal)  maps  which  preserve 
the  markings  might  be.  If  we  fix  a  Riemann  surface  R  and  a  marking  E  on  R,  the 
distance  between  two  points  Ri  =  fi{R)  and  R2  =  f2 (R)  in  the  Teichmiiller  space  of 
R  is  given  by 

d(Ri,R>)  =  llog(AT*), 

where  K*  is  the  infimum  of  the  dilatation  of  g2  0  gfl  where  gi  and  g2  are  equivalent 
to  fi  and  f2,  respectively.  This  infimum  is  attained,  by  definition,  by  the  unique 
Teichmiiller  map. 

Remark  2.4.1.  The  Teichmiiller  maps  on  the  space  of  tori  are  linear  transformations. 
That  is,  the  Teichmiiller  map  from  the  torus  T  to  the  torus  T  is  affine. 

2.5  Earthquakes  on  Surfaces 

Now,  we  describe  the  action  of  earthquakes  on  compact  Riemann  surfaces  with 
genus  g  >  2  (i.e.,  n-holed  tori,  n  >  2)  and  their  Teichmiiller  spaces.  To  do  so,  we 
must  describe  the  constructions  on  these  surfaces  necessary  to  the  earthquake  actions. 
That  is,  we  need  measured  geodesic  laminations  consisting  of  simple  closed  geodesics 
on  the  surface  (though  they  will  not  necessarily  be  finite)  and  hyperbolic  shearing 
and  grafting  maps. 

Let  R  be  a  compact  Riemann  surface  with  genus  g  >  2,  (i.e.,  R  is  a  compact, 
hyperbolic  Riemann  surface).  Just  as  in  the  disk,  a  finite  geodesic  lamination  on  R 
is  a  finite,  pairwise  disjoint  collection  of  geodesics  in  the  intrinsic  metric  on  R,  but 
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how  do  we  define  these  geodesics?  Well,  we  know  that  a  geodesic  in  the  Poincare 
model  of  the  hyperbolic  plane  is  an  arc  of  a  circle  which  intersects  dB  orthogonally. 
Further,  we  know  that  for  any  hyperbolic  Riemann  surface,  the  universal  cover  for 
that  surface  is  the  hyperbolic  plane  where  R  is  related  to  the  universal  cover  through 
its  conformal  structure  [14,  15,  29,  30,  45].  So,  given  a  curve  7  C  R,  7  is  a  geodesic  on 
R  if  and  only  if  7  lifts  to  a  geodesic  in  the  hyperbolic  plane  [45].  Thus,  a  lamination 
on  the  Riemann  surface  R  lifts  to  a  geodesic  lamination  in  the  hyperbolic  plane, 
as  illustrated  in  Figure  2.14,  where  we  lift  a  geodesic  lamination  from  a  compact 
surface  of  genus  2  to  the  disk.  The  tiling  is  accomplished  through  the  application  of 
transition  maps  which,  in  this  case,  are  hyperbolic  isometries  with  axes  perpendicular 
to  corresponding  sides  of  the  fundamental  region.  The  lift  of  geodesic  segments  are 
shown  as  dashed  curves.  For  more  detailed  information  on  the  structure  of  geodesics 
and  laminations  on  surfaces,  see  the  work  of  Bonahon  [5,  7].  Take  a  collection  of 
simple  closed  geodesics  in  R,  and  place  on  each  geodesic  in  R  a  real  (for  shearing), 
imaginary  (for  grafting),  or  complex  (for  shearing  then  grafting)  weight.  This  gives  us 
a  measured  geodesic  lamination  (£,  a)  on  R.  Note,  however,  that  while  we  began  with 
a  finite  lamination  on  R,  through  the  lift  we  obtain  a  countably  infinite  lamination 
on  D. 


(a)  A  Fundamental  Region 


(b)  One  Layer  of  Tiling 


(c)  Two  Layers  of  Tiling 


Figure  2.14:  Geodesics  on  Surfaces  Lift  to  Geodesics  on  D 


Now,  we  can  define  an  earthquake  on  R  in  terms  of  the  mechanisms  already 
described  in  Section  2.2  and  Section  2.3.  A  finite  measured  geodesic  lamination  on 
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R  lifts  to  a  measured  geodesic  lamination  on  O  comprising  a  sequence  of  geodesics 
{£n}„GN-  To  this  point,  we  have  only  defined  finite  earthquakes  on  the  disk,  so  how 
do  we  define  an  earthquake  on  a  countably  infinite  set  of  geodesics?  We  do  so  just  as 
a  limit  of  finite  earthquakes. 

On  the  surface  itself,  this  earthquake  action  has  a  nice  geometric  interpretation. 
Shearing  along  a  geodesic  in  the  disk  induced  by  a  geodesic  on  the  hyperbolic  Riemann 
surface  R  is  equivalent  to  simply  cutting  open  the  surface  along  the  geodesic,  twisting 
by  an  amount  prescribed  by  the  weight  on  the  geodesic,  and  gluing  the  two  ends  back 
together,  as  illustrated  in  Figure  2.15.  Similarly,  grafting  along  a  geodesic  in  the  disk 
induced  by  a  geodesic  on  the  surface  R  is  equivalent  to  cutting  open  the  surface  along 
the  geodesic,  inserting  a  hyperbolic  cylinder  with  a  height  prescribed  by  the  weight 
on  the  geodesic,  and  gluing  the  ends  together,  as  illustrated  in  Figure  2.16. 


(a)  Cut  Along  a  Geodesic  and  Twist  (b)  Glue  the  Ends  Together 

Figure  2.15:  Shearing  Action  on  a  Hyperbolic  Surface 

The  question  now  becomes  whether  or  not  the  resulting  object  is  a  Riemann 
surface.  Are  the  objects  illustrated  in  Figure  2.15b  and  Figure  2.16b  Riemann  sur¬ 
faces?  That  is,  is  the  action  of  an  earthquake  on  a  hyperbolic  surface  a  map  in  the 
Teichmiiller  space  of  that  surface? 

To  answer  this  question,  we  now  turn  to  some  technical  definitions  and  results 
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(a)  Cut  Along  a  Geodesic  and  “Separate”  (b)  Insert  a  Cylinder  and  Glue  Together 


Figure  2.16:  Grafting  Action  on  a  Hyperbolic  Surface 


regarding  the  action  of  earthquakes,  culminating,  in  Section  2.6,  in  the  explicit  con¬ 
struction  of  a  conformal  structure  on  the  image  of  a  Riemann  surface  under  an  earth¬ 
quake  (involving  positive  and  negative  shearing  and  positive  grafting) ,  demonstrating 
that  earthquakes  are  maps  in  the  Teichmiiller  space  of  the  hyperbolic  surface  R. 


Definition  2.5.1.  An  orientation-preserving  homeomorphism  (f) :  5B  — >  9B  is  called 

a  k- quasisymmetry  if  and  only  if 


<p(e$x+t'))  -  <p(eix) 
<p(eix)  —  <p{e^x~ *)) 


<  k, 


(2.5.1) 


for  all  x  and  for  all  t  G  (— §,  §)• 

Quasisymmetries  and  quasiconformal  maps  are  closely  related  in  that  the  boundary 
values  of  a  quasiconformal  map  give  a  quasisymmetry,  and  a  quasisymmetry  on  <9B 
induces  a  quasiconformal  map  on  B. 


Definition  2.5.2.  Let  (£,  a)  be  a  measured  geodesic  lamination,  not  necessarily  finite, 
in  B;  say  C  =  {Ln}.  Let  I  be  a  curve  in  B  with  unit  hyperbolic  length  such  that  I  f]Ln 
contains  no  ( non-empty )  open  subsets  for  every  Ln  €  C,  and  I  intersects  any  geodesic 
Ln  €  C  at  most  finitely  many  times.  Define  the  quantity  m(I )  as 

m(I )  -  ^2  (kn  mod  2)  k(£n)|,  (2.5.2) 

neN/ 
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where  N/  =  {n  €  N  :  I  f)Ln  ^  0}  and  kn  is  the  number  of  times  I  intersects  Ln.  The 
transverse  measure  of  the  measured  geodesic  lamination  (£,  a)  is  given  by 

ra  =  sup{m(/)} .  (2.5.3) 

/CD 

In  computing  the  quantity  m(I),  the  first  time  we  cross  a  geodesic,  we  call  that 
orientation  positive,  and  for  each  successive  crossing  of  the  geodesic  in  that  direction 
we  add  the  associated  weight.  Crossing  in  the  other  direction  is  given  a  negative 
orientation,  and  such  a  crossing  subtracts  the  associated  weight  from  our  measure. 
Thus,  crossing  a  geodesic  an  odd  number  of  times  adds  the  associated  weight  to  m(/) 
and  crossing  an  even  number  of  times  adds  zero. 

Now,  we  establish  the  intimate  link  between  quasisymmetries  and  the  transverse 
measure  of  a  measured  lamination  on  B.  One  direction  of  Theorem  2.5.1  was  proven 
by  Thurston  [49],  and  versions  of  the  complete  theorem  have  been  proven  indepen¬ 
dently  by  Gardiner,  Hu,  and  Lakich  [23],  Gardiner  and  Lakic  [24],  Hu  [25],  and  Saric 
[40]. 

Theorem  2.5.1.  Let  (£,  cr)  be  a  measured  geodesic  lamination  on  D  and  let  E  be 
the  earthquake  on  D  induced  by  (£,  a) .  The  continuous  extension  of  E  to  c®  is  a 
quasisymmetry  if  and  only  if  the  transverse  measure  of  (£,  a)  is  bounded. 

Note  that  a  single  measured  geodesic  lamination  may  give  rise  to  more  than  one  earth¬ 
quake,  as  illustrated  in  Example  2.2.2.  As  stated  in  Proposition  2.2.2,  however,  the 
extension  of  these  earthquakes  to  their  boundaries  differ  by  a  Mobius  transformation. 
As  a  conformal  map,  this  Mobius  transformation  has  no  effect  of  the  nature,  in  terms 
of  the  quasisymmetry  constant  k,  of  the  boundary  quasisymmetry.  That  is,  we  are 
justified  in  speaking  of  the  unique  earthquake  E  in  Theorem  2.5.1.  Thus,  an  appro¬ 
priately  measured  lamination  gives  an  earthquake  which  induces  a  quasisymmetry  on 
dB.  This  quasisymmetry  induces,  in  turn,  a  quasiconformal  map  on  B.  Now,  this  is 
how  we  described  the  relationship  between  points  in  Teichmiiller  space. 
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2.6  A  Conformal  Structure  on  the  Images  of  Hyperbolic 
Riemann  Surfaces  under  Complex  Earthquake  Maps 

Let  R  be  a  compact,  hyperbolic  Riemann  surface,  and  suppose  that  we  apply  to  R 
a  finite  complex  earthquake  E  induced  by  a  finite  measured  geodesic  lamination  (£,  a) 
on  R,  where  the  lamination  consists  of  n  geodesics  with  complex  weights  a(Li)  = 
Hi  +  i\i,  1  <  i  <  n,  where  Hi  £  ®  and  A*  €  R  such  that  A*  >  0.  This  earthquake 
thus  involves  both  shearing  maps  and  (positive)  grafting  maps  composed  to  form  the 
complex  earthquake. 

Since  R  is  a  fixed  Riemann  surface,  it  has  a  conformal  structure  associated  with 
it;  say  {< pv}v<=.y  is  the  atlas  for  the  conformal  structure  on  R.  Further,  since  each 
grafting  map  involves  the  insertion  of  a  cylinder  Cj,  we  have  a  conformal  structure 
associated  with  each  such  cylinder  such  that  the  coordinate  charts  map  c,  to  subsets 
of  a  hyperbolic  wedge.  We  do  this  by  mapping  first  to  an  infinite  strip  and  then 
applying  the  exponential  map  ez  to  that  strip.  For  each  of  Cj,  0  <  j  <  m,  where 
m  gives  the  number  of  geodesics  in  the  lamination  C  which  have  non-zero  imaginary 
components  for  their  weights,  let  {</){} k(_  K  be  the  atlas  associated  with  each  such 
cylinder. 

Define  a  pullback  map  P  on  R,  the  image  of  R  under  the  complex  earthquake  E, 
so  that  the  image  of  a  point  r  G  R  under  P  is  the  unique  point  (on  R  or  on  one  of 
the  cylinders  Cj )  whose  image  under  the  complex  earthquake  E  is  r. 

With  this  information  we  now  define  a  conformal  structure  on  the  Riemann  surface 
R.  To  define  this  conformal  structure  we  need  only  describe  a  set  of  coordinate  charts 
mapping  open  regions  in  the  surface  R  to  C  so  that  the  transition  maps  associated 
with  the  coordinate  charts  are  analytic.  It  is  sufficient  to  construct  for  every  point 

- - -  /S 

r  £  R  a  map  from  an  open  neighborhood  U  C  R  to  C  so  that  these  maps  satisfy 
the  analyticity  condition.  It  is  sufficient  to  consider  open  sets  of  the  following  four 
classes. 

1.  U a  is  the  collection  of  all  open  sets  Ua  C  R  such  that  Ua  does  not  intersect  the 
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image  under  the  earthquake  E  of  any  geodesic  Li,  1  <  i  <  n,  and  is  disjoint 
from  every  inserted  cylinder  E(cj),  1  <  j  <  m. 

2.  UB  is  the  collection  of  all  open  sets  Up  C  R  such  that  Up  C  E(cj )  for  some 
j,  1  <  3  <m. 

3.  Ur  is  the  collection  of  all  open  sets  C/7  C  R  such  that  f/7  intersects  E(Li)  for 
exactly  one  value  i,  1  <  i  <  n,  where  a  (Li)  £  R,  (i.e.,  the  only  map  on  that 
geodesic  is  a  shearing  operation). 

4.  UA  is  the  collection  of  all  open  sets  Us  C  R  such  that  Us  intersects  both  P~1(R ) 
and  E(cj)  for  some  j,  1  <  j  <  m. 

An  illustration  of  representative  sets  from  these  classes  is  given  in  Figure  2.17. 


Case  I:  Let  Ua  £  Ua,  and  note  that  P(Ua )  is  an  open  set  in  R.  From  the  conformal 
structure  on  R  and  corresponding  to  this  open  set  P(Ua )  we  have  a  coordinate  chart 
< pa  such  that  ipa{P(Ua))  C  H  =  {z  £  C  :  Im(2:)  >  0}.  We  thus  define  a  coordinate 
chart  ijja  :  R  — »  C  on  the  open  set  Ua  by  ipa  =  (pa°  P-  In  this  way  we  define  a  family 
of  coordinate  charts  on  R  by  =  {il>a}aeA- 

Case  II:  Let  Up  £  Ub ■  In  this  case,  Up  C  Int(£’(cj)),  for  some  1  <  j  <  m,  and 
the  pullback  map  P  takes  this  set  to  an  open  set  on  the  finite  cylinder  Cj  as  shown  in 
Figure  2.18a.  Now,  from  the  conformal  structure  on  Cj  we  have  a  coordinate  chart  <j^ 
corresponding  to  P(Up)  such  that  (fPp(P(Up ))  is  a  subset  of  an  infinite  hyperbolic  strip 
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of  uniform  width  determined  by  the  complex  part  of  the  weight  on  the  appropriate 
geodesic  in  £.  We  then  have  a  grafting  map  gjp  that  opens  the  upper  half-plane 
along  the  imaginary  axis  and  glues  in  the  hyperbolic  strip  as  shown  in  Figure  2.18b. 
Following  this  with  the  conformal  map  M,  that  transforms  the  result  into  the  upper 
half-plane,  we  have  a  map  :  R  ->•  C,  shown  in  Figure  2.18c,  on  the  open  set  Up 
defined  by  x/jp  =  Mj  o^o^oP.  In  this  way  we  define  a  family  of  coordinate  charts 
on  R  by  =  {ipp}peB- 


/' . \ 

(  i  P(U0) 

\ 

(a)  Image  of  Up  C  R  under  P  (b)  Image  of  Up  C  R  under  gjfj  o  ^  o  P 


Figure  2.18:  Construction  of  Coordinate  Charts  for  Open  Sets  Up 

Case  III:  Let  C/7  e  Ur,  and  note  that  P(t/7)  is  a  set  as  shown  in  Figure  2.19a.  From 
the  conformal  structure  on  R  and  corresponding  to  an  open  set  containing  P(U7)  we 
have  a  map  <p^  such  that  t/97(P(I/7))  c  H.  Now,  we  have  disjoint  sets  that  differ 
by  the  application  of  a  hyperbolic  shearing  map,  Sy,  along  the  vertical  axis  shown 
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in  Figure  2.19b.  We  have  thus  defined  a  coordinate  chart  ipy  :  R  -»  C  on  the  open 
set  U, y  by  r/>7  =  S1  o  o  P.  In  this  way,  we  define  a  family  of  coordinate  charts 

=  {^7}7er- 


(c)  Image  of  U7  C  R  under  57  o  ipy  o  p 
Figure  2.19:  Construction  of  Coordinate  Charts  for  Open  Sets  Uy 


Case  IV:  Let  Us  €  U&,  and  note  that  P(U$)  is  divided  between  two  distinct 
Riemann  surfaces  Cj,  as  shown  in  Figure  2.20a,  and  R,  as  shown  in  Figure  2.20b.  To 
define  coordinate  charts  on  set  of  the  form  Us  we  define  a  pair  of  maps  from  these 
surfaces  to  the  upper  half-plane  H.  Beginning  with  the  cylinder,  we  have  a  map  (jxg 
such  that  <f>l(P(E(cj)  D  Us))  is  a  hyperbolic  strip  of  uniform  width  determined  by 
the  complex  part  of  the  weight  on  the  appropriate  geodesic  in  C.  We  then  have  a 
grafting  map  g$  that  opens  the  upper  half-plane  along  the  imaginary  axis  and  glues  in 
the  hyperbolic  strip  as  shown  in  Figure  2.20c.  At  the  same  time,  from  the  conformal 
structure  on  R  we  have  a  map  ips  such  that  (ps(P((E(cj))c  fl  Us))  C  H  as  shown  in 
Figure  2.20d.  Applying  the  grafting  map  g\  we  put  the  pieces  back  together  in  C. 
Now,  we  may  need  to  shear  along  the  hyperbolic  geodesic  whose  image  to  this  point  is 
the  imaginary  axis  in  Figure  2.20d;  this  will  depend  upon  whether  or  not  the  geodesic 
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associated  with  this  grafting  action  had  a  weight  with  non-zero  real  part.  Note  that 
the  shearing  will  take  place  along  that  portion  of  the  set  not  intersecting  E(cj ),  since 
we  define  complex  earthquakes  as  the  composition  of  shearing  and  grafting  maps,  in 
that  order.  Finally,  we  know  that  there  exists  a  conformal  map  Mj  that  takes  the 
opened  half-plane  to  the  upper  half-plane. 

,  .  ,  f  Mj  ogl  0^0  P(r)  if  re  E(Cj)  n  Us  C  R 

Mr)  =  \  _  (2.6.1) 

[  Mj  o  Ss  o  $  o  (ps  o  P(r)  if  r  G  ( E(cj))c  D  Us  C  R 

This  function  now  defines  a  map  ipg  :  R  — ¥  C  on  the  open  set  Us,  as  shown  in 
Figure  2.20f.  Taking  all  such  functions  over  the  set  of  possible  sets  Ug  G  17a,  we  define 
a  family  of  coordinate  charts  =  {^}jeA- 

Define  a  collection  of  open  sets,  U  =  Ua\JUb{J  Ur  (J  Ua ,  and  a  collection  of  maps 
from  R  to  C,  $  =  ^  (J  U 

Proposition  2.6.1.  The  transition  maps  in  the  structure  (17,  \k)  are  analytic. 

Proof.  Let  U  be  the  intersection  of  two  open  sets  in  U.  We  need  only  consider  those 
intersections  for  which  U  G  U.  This  gives  four  cases  for  the  overlap  region  U  in  which 
we  must  verify  the  analyticity  of  the  transition  maps. 

Case  I:  Suppose  that  U  G  Ua-  There  are  four  possible  ways  in  which  U  may  occur 
as  the  result  of  intersections  of  sets  in  U. 

Case  la:  U  =  Uai  f]Ua2,  where  Uai,Ua2  £  Ua-  Associated  with  Uai  and  Ua2  we 
have  coordinate  charts  ipa l,'ipa2  e  ^  a-  Say  ipai  =  ipai  o  P  and  i>a2  =  ipa2  o  P.  Now, 
consider  the  transition  maps  ipa2  o  and  ipai  o  . 

°  =  (<Pa2  °P)°  (¥>ai  °  P)~l 

=  ipa2  O  P  O  P-1  O 

=  T«2  °  Ta] 

Thus,  ipa2  oip-i  :  C  — »•  C  is  analytic,  since  pai,  i  =  1,2,  were  taken  from  the  conformal 
structure  on  R,  making  (pa2  o  analytic.  Similarly,  <pai  o  (p~*  is  analytic. 
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(c)  Image  of  P(Us  p)(I?(cj)))  under  4>l  (d)  Image  of  P(Ug  f](E(cj))c)  under  ipg 


(e)  Image  under  the  Grafting  g3s  (f)  Image  of  Us  C  R  under  ips 

Figure  2.20:  Construction  of  Coordinate  Charts  for  Open  Sets  Us 


Case  lb:  U  =  Uaf]  t/7,  where  Ua  G  Ua  and  U1  G  Ur-  Associated  with  Ua  and  U1 
we  have  coordinate  charts  x\)a  G  and  G  Wp,  respectively.  Say  ipa  =  ipa  o  P  and 
t/>7  =  S7  o  tp1  o  P.  Now,  consider  the  transition  maps  ipa  o  t\)~x  and  ip1  o  ip~l . 

O  ll>~' 1  =  ((Pa  oP)o  (S7  0(p70  P)_1 

=  ^oFoF-'o^o^1 
=  (pa  O  y~l  O  S~l 

Note  that  the  transition  map  ipQ  o  xp~l  is  defined  only  on  the  set  U.  Since  the  set  U 
does  not  itself  intersect  the  geodesic  associated  with  the  shearing  map  57,  on  U  the 
shearing  map  acts  as  either  the  identity  map  or  a  hyperbolic  Mdbius  transformation. 
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In  either  case,  the  action  of  57  is  (locally)  analytic.  Next,  we  note  that  since  <pQ 
and  v?7  are  taken  from  the  conformal  structure  on  R,  ipa  o  is  analytic,  where 
defined.  Thus,  the  transition  map  ^  °  is  the  composition  of  analytic  maps,  and 
is  therefore  analytic.  Similarly,  ip7  o  ip~l  is  also  analytic. 

Case  Ic:  U  —  Uaf]Us,  where  Ua  G  Ua  and  Us  G  Ua-  Associated  with  each 
of  Ua  and  U$  we  have  coordinate  charts  ipa  G  and  ips  G  respectively.  Say 

ipa  =  V’a  °  P  and  ips  of  the  form  given  in  (2.6.1).  We  note  that  since  U  does  not 
intersect  the  geodesic  associated  with  the  grafting  of  the  cylinder,  for  every  r  G  U 
we  have  r  G  (E(cj))c  D  U$  C  R,  where  j  denotes  the  cylinder  associated  with  the 
grafting.  Thus,  on  U  we  have  ips  =  Mj  o S$  o^o^o P.  Now,  consider  the  transition 
maps  ipaoip~l  and  ipsoip~l- 

ipaorpj1  =  (cpa  °  P)  o  (Mj  o  Ss  o  gj  o  cpg  o  P)-1 

=  (pa°PoP~lo  (pj1  o  (gl)  1  o  Sjl  o  M-1 

=  V«0^l°{g3s)  loSp;loMrl 

Note  that  we  define  the  transition  map  ipa  o  ipppl  only  on  the  overlap  region  U.  Since 
the  set  U  does  not  itself  intersect  the  inserted  cylinder  (or  the  boundary  along  which 
we  weld  the  cylinder  to  the  Riemann  surface  R  via  the  grafting  map  (g°6)  '),  on  U 
the  inverse  of  the  grafting  map  acts  as  either  the  identity  map  or  as  multiplication  by 
an  exponential  map.  In  either  case,  the  action  of  (gj)  1  is  (locally)  analytic.  For  the 
same  reason,  on  U  we  have  that  S^1  acts  as  either  the  identity  map  or  a  hyperbolic 
Mobius  transformation,  and  is  therefore  (locally)  analytic.  Next,  we  note  that  since 
the  maps  cpa  and  <ps  were  borrowed  from  the  conformal  structure  on  the  original 
Riemann  surface  R,  <pa  °  is  analytic,  where  defined.  Finally,  the  map  M  was 
chosen  as  the  conformal  map  taking  the  region  shown  in  Figure  2.20e  to  the  upper 
half  plane.  Thus,  the  transition  map  ipa  o  ip~l  is  the  composition  of  analytic  maps  on 
U,  and  is  therefore  analytic.  Similarly,  ip6  o  ip~l  is  also  analytic. 

Case  Id:  U  =  C/7  f)  Us,  where  C/7  G  Ur  and  Us  G  Ua-  Associated  with  each 
of  U7  and  U$  we  have  coordinate  charts  07  G  Tr  and  ips  G  Ta,  respectively.  Say 
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ipy  =  S7o^oF  and  ip#  of  the  form  given  in  (2.6.1).  We  note  that  since  U  intersects 
neither  the  geodesic  associated  with  the  grafting  of  the  cylinder  nor  the  geodesic 
associated  with  the  simple  shearing,  for  every  ret/  we  have  r  G  (E{cj))c  fi  Ug  C  R, 
where  j  denotes  the  cylinder  associated  with  the  grafting.  Thus,  on  U  we  have 
ipg  =  Mj  o  Sg  o  g3s  o  cps  o  P.  Now,  consider  the  transition  maps  ip-yOipj1  and  ips  °  V;7 1- 

o  ipj1  =  (57  o<p1oP)o  {Mj  o  Ss  o  g35  o  (ps  o  py1 

=  S7  O  ^7  o  p  o  p-1  o  cpj1  o  (g3)  1  o  Sg1  o  Mr1 
=  S^ocp^o  (fij1  o  (, g3 )  1  o  Sg1  O  Mr1 

Note  that  we  define  the  transition  map  •07  o  ip5l  only  on  the  overlap  region  U.  Since 
the  set  U  does  not  itself  intersect  the  inserted  cylinder  (or  the  boundary  along  which 
we  weld  the  cylinder  to  the  Riemann  surface  R  via  the  grafting  map  {g3s)  1),  on  U 
the  inverse  of  the  grafting  map  acts  as  either  the  identity  map  or  as  multiplication 
by  an  exponential  map.  In  either  case,  the  action  of  (g^)  1  is  (locally)  analytic.  For 
the  same  reason,  on  U  we  have  that  Sg1  and  57  each  act  as  either  identity  maps  or 
hyperbolic  Mobius  transformations;  they  are  thus  (locally)  analytic.  Next,  we  note 
that  since  the  maps  </?7  and  <pg  were  borrowed  from  the  conformal  structure  on  the 
original  Riemann  surface  R,  ipa  o  pj1  is  analytic,  where  defined.  Finally,  the  map 
M  was  chosen  as  the  conformal  map  taking  the  region  shown  in  Figure  2.20e  to  the 
upper  half  plane.  Thus,  the  transition  map  -07  o  'ipg1  is  the  composition  of  analytic 
maps  on  U,  and  is  therefore  analytic.  Similarly,  ipg  o  ip~l  is  also  analytic. 

Case  II:  Suppose  that  U  e  Ub-  There  are  two  ways  in  which  U  may  occur  as  the 
result  of  intersections  of  sets  in  U. 

Case  Ha:  U  =  Up1  f)  Up2,  where  Up1,Up2  G  Ub-  Associated  with  Up1  and  Up2 
we  have  coordinate  charts  ipp1 , ipp2  G  fg.  Say  ippl  =  Mj  o  g3^  o^oP  and  ipp2  = 
Mj  o  g3^  o  <pi2  o  P.  Note  that  since  both  Up1  and  Up2  must  necessarily  lie  in  the 
same  finite  inserted  cylinder,  the  grafting  maps  g3^  and  g3^,  are  the  same  map.  Now, 
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consider  the  transition  maps  ipp2  o  ipp'  and 

°  =  (Mj  o  g^  o  (f)^  oP)o  (Mj  o  gjji  o  ^  o  P)”1 

=  MjO^o^oP  oP~l  o  (0Ji)  1  o  (^jx)  1  o Mr1 

=  Mj  o  gfh  o  y  o  (^y1  o  fay1  o  M~l 

=  Mi  °  9k  °  (fh  °  Wi)"1)  °  (<)_1  °  MIl 

Now,  notice  that  the  maps  Mj  and  M~l  are  analytic.  Next,  notice  that  the  maps  gi 
and  gp2  have  no  effect  on  the  analyticity  of  the  maps  involved,  since  both  Up1  and  Up2 
are  on  the  interior  of  the  cylinder  E(cj).  These  maps  merely  indicate  that  we  consider 
the  infinite  hyperbolic  strip  in  a  different  context.  The  action  affecting  analyticity 
takes  place  in  the  conformal  structure  on  the  cylinder  Cj,  and  (pjh  o  (4>k)  1  is  analytic 
since  these  maps  were  taken  from  this  conformal  structure.  Thus,  as  a  composition 
of  analytic  maps,  ipp2  o  ipy  is  analytic.  Similarly,  ipp1  o  ipy  is  analytic. 

Case  lib:  U  —  Up  Q  Us,  where  Up  G  UB  and  Us  G  Ua-  Associated  with  each 
of  Up  and  Us  we  have  coordinate  charts  ipp  G  T#  and  ips  G  TA,  respectively.  Say 
ipp  —  Mj  o  g°p  o  o  P  and  ips  of  the  form  given  in  (2.6.1).  Note  that  since  U  does  not 
intersect  any  geodesic  associated  with  the  grafting  of  the  cylinder,  for  every  r  G  U  we 
have  r  G  E(cj)  n  Us  C  R,  where  j  denotes  the  cylinder  associated  with  the  grafting. 
Thus,  on  U  we  have  ips  =  Mj-o^o^oP.  Now,  consider  the  transition  maps  ippoipj1 
and  ips  o  ipp1. 

ippoip-1  =  (Mj  o  g^  o  ftp  o  P)  o  (Mj  o  g*  o  $  o  P)-1 

This  case,  thus  reduces  to  Case  Ha,  and  the  transition  maps  ipp  o  ipj1  and  ips  o  ipp 1 
are  analytic. 

Case  III:  Suppose  that  U  G  Ur.  Then,  U  =  t/7l  p|  U12 ,  where  C/7l ,  Ul2  G  Ur- 
Associated  with  each  of  t/7l  and  Ul2  we  have  coordinate  charts  iplx ,  ipl2  G  'kp.  Say 
ipn  =  S7  o  (/?7l  o  P  and  tpl2  =  S7  o  ipri  o  P.  Note  that  the  earthquake  (shearing)  map 
in  each  of  these  coordinate  charts  is  the  same,  since  this  earthquake  is  defined  by 
intersection  with  a  particular  weighted  geodesic  rather  than  a  particular  set  U.  Now, 
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consider  the  transition  maps  ?/>72  o  ^7  1  and  ipyi  o  ^  1. 

^72  °  =  (Sy  O  ipl2  oP)o  (Sy  o  <^71  O  P)-1 

=  SyOipl2oPo  p-1  o  p"1  o  S'-1 
=  57  O  (py2  O  o  s~ 1 

The  action  of  the  transition  map  here  reduces  to  a  hyperbolic  isometry,  and  is  thus 
analytic. 

Case  IV:  Suppose  that  U  G  UA.  U  =  Us1f]Ug2,  where  Us1,Ug2  G  Pa-  Associated 
with  each  of  U$1  and  Us2  we  have  coordinate  charts  'ips1,i’s2  €  Say  and  ipg2 
are  each  of  the  form  given  in  (2.6.1).  The  set  U  can  now  be  decomposed  into  three 
disjoint  sets,  U  =  Ua\jL\jUp,  where  Ua  G  Ua,  Up  G  Ub,  and  L  =  d  Uaf]d  Up, 
as  shown  in  Figure  2.21.  We  have  already  shown  that  the  transition  maps  ipsx  °  Vv/ 
and  tpg2  o  are  analytic  on  each  of  the  sets  Ua  and  Up.  These  transition  maps, 
though,  are  also  clearly  continuous  on  all  of  U.  Since  the  image  of  L  C  U  in  the  plane 
has  measure  zero,  any  map  which  is  if-quasiconformal  on  the  image  of  U  \  L  is  also 
P-quasiconformal  on  L  [33];  special  cases  of  this  result  are  often  given  in  texts  on 
complex  analysis  [18].  Thus,  since  the  transition  maps  are  analytic  (1-quasiconformal) 
on  the  images  of  Ua  and  Up,  the  transition  maps  are  analytic  on  all  of  U. 


Figure  2.21:  Partition  of  U  into  Ua,  Up,  and  L 
Thus,  the  transition  maps  in  the  structure  (U,  T)  are  analytic. 


□ 
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Since  the  structure  defined  above  describes  a  cover  for  the  surface  R  and  a  set 
of  maps  on  that  cover  so  that  the  transition  maps  are  analytic,  we  have  defined  a 
sufficient  conformal  structure  on  the  surface  R\  we  simply  include  any  other  sets  and 
maps  compatible  with  this  structure. 
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CHAPTER  III 
CIRCLE  PACKING 


3.1  Preliminaries  and  Definitions 

A  circle  packing  is  a  configuration  of  circle  with  a  prescribed  pattern  of  tangencies. 
William  Thurston  conjectured  in  1985  that  these  circle  packings  might  be  used  to 
approximate  the  action  of  conformal  maps  [50] .  These  circle  packings  have  since  been 
widely  studied,  with  applications  in  many  different  areas  of  mathematics.  We  begin 
here  with  some  basic  definitions  and  a  general  discussion  of  circle  packing.  Several 
excellent  resources  are  available  with  much  greater  detail  [20,  42,  45,  46]. 

Definition  3.1.1.  A  vertex  is  called  a  O-simplex.  An  edge  is  called  a  l-simplex. 
A  face  is  called  a  2-simplex. 

We  use  these  simplexes  to  describe  a  combinatorial  object  called  a  2-complex. 

Definition  3.1.2.  A  simplicial  2-complex  JC  is  a  topological  space  represented 
as  a  countable  (possibly  finite )  union  of  0-,  1-,  and  2-simplexes  with  the  following 
properties. 

1.  Every  face  of  a  simplex  in  JC  is  itself  a  simplex  in  JC. 

2.  2-simplexes  of  JC  intersect  either  in  the  empty  set  or  in  a  simplex  of  JC. 

3.  No  0-simplex  (vertex)  of  JC  is  in  more  than  a  finite  number  of  simplexes  of  JC. 
4-  JC  is  connected. 

We  will  consider  a  subset  of  all  possible  simplicial  2-complexes  in  which  the  2- 
simplexes  (faces)  are  combinatorial  triangles. 

Definition  3.1.3.  A  bounded  degree  abstract  triangulation  JC  is  an  abstract 
simplicial  2- complex  which  triangulates  an  orientable  topological  surface  such  that 
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1.  the  set  of  interior  vertices  (vertices  such  that  every  incident  edge  belongs  to  two 
faces )  is  non-empty  and  edge-connected; 

2.  no  interior  edge  (an  edge  belonging  to  two  faces)  in  K,  has  both  vertices  on  the 
boundary; 


3.  no  vertex  in  JC  belongs  to  more  than  two  boundary  edges; 


4.  there  is  an  upper  bound  on  the  degree  of  vertices  in  1C. 


It  is  this  combinatorial  object,  the  abstract  triangulation,  which  encapsulates  the 
“prescribed  pattern  of  tangencies”  in  our  circle  packing.  We  refer  to  these  trian¬ 
gulations  as  abstract  to  emphasize  the  fact  that  in  the  definition  we  have  implied 
no  concrete  geometric  realization.  A  2-complex  and,  by  extension,  the  associated  ab¬ 
stract  triangulation  are  purely  combinatorial  objects;  they  have  no  inherent  geometric 
structure  until  they  are  realized  as  a  circle  packing. 


(a)  An  Abstract  Triangulation 


(c)  Violates  Condition  2. 


(b)  Violates  Condition  1. 


Figure  3.1:  An  Abstract  Triangulation  and  Pathological  Triangulations 


An  example  of  a  simple  valid  abstract  triangulation  is  shown  in  Figure  3.1;  we  also 
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show  how  triangulations  may  fail  to  meet  conditions  (1)  through  (3)  from  Definition 
3.1.3.  In  Figure  3.1b  we  see  two  interior  vertices  which  are  not  edge  connected; 
in  Figure  3.1c  we  see  an  edge  whose  endpoints  are  both  boundary  vertices;  and 
in  Figure  3. Id  we  see  a  vertex  belonging  to  four  boundary  edges.  Note  that  since 
the  triangulations  shown  are  all  finite,  they  must  also  satisfy  the  bounded  degree 
condition. 

Definition  3.1.4.  A  circle  packing  is  a  configuration  of  circles  with  a  specified 
pattern  of  tangencies.  In  particular,  if  /C  is  an  abstract  triangulation  of  a  topological 
surface,  then  a  circle  packing  P  for  K  is  a  configuration  of  circles  such  that 

1.  P  contains  a  circle  Cv  for  every  vertex  v  £  1C; 

2.  if  [u,  v]  is  an  edge  of  1C,  then  Cv  is  externally  tangent  to  Cu; 

3.  if  (v,  u,  w)  is  a  positively  oriented  face  of  1C,  then  (Cv,  CU,CW)  forms  a  positively 
oriented  mutually  tangent  triple  of  circles  in  P. 

Realizing  such  a  configuration  is,  at  its  most  essential,  a  problem  in  computing  the 
necessary  radii  of  each  circle  in  the  packing.  Examples  of  such  algorithms  are  given 
by  Collins  and  Stephenson  [17,  42]  and  Mohar  [37].  A  simple  example  of  an  algorithm 
implemented  in  Mathematica  for  computing  the  radii  is  also  given  in  Appendix  D. 
All  of  these  algorithms  amount  to  numerically  solving  a  boundary  value  problem  in 
which  the  radii  of  the  boundary  circles  are  fixed  and  the  remaining  radii  are  computed 
algorithmically.  A  circle  packing  is  called  univalent  if  the  circles  in  the  packing  have 
mutually  disjoint  interiors.  That  is,  the  packing  is  univalent  if  no  two  circles  intersect 
in  more  than  one  point.  This  univalent  circle  packing  represents  a  geometric  real¬ 
ization  of  the  underlying  abstract  triangulation  K.  Vertices  in  the  triangulation  may 
be  realized  in  this  packing  as  the  centers  (in  some  particular  geometry,  hyperbolic, 
Euclidean,  or  spherical)  of  the  circles,  and  the  edges  as  geodesic  segments  connecting 
the  centers.  This  embedding  is  called  the  carrier  of  the  circle  packing. 
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If  JC  is  embedded  in  C  in  two  different  ways  (e.g.,  by  giving  two  different  sets  of 
values  for  the  radii  of  the  boundary  circles),  there  is  a  natural  piecewise  map  from  the 
carrier  associated  with  one  packing  to  the  other  achieved  by  sending  triangles  of  one 
packing  to  their  counterparts  in  the  other  using  affine  maps.  These  piecewise  affine 
maps  are  referred  to  as  discrete  conformal  maps. 

In  Figure  3.2  we  have  a  graphic  description  of  a  bounded  degree  abstract  trian¬ 
gulation,  K.  This  triangulation  comprises  27  vertices,  with  edges  between  vertices 
indicating  tangencies. 


We  can  produce  a  geometric  realization  (as  a  circle  packing)  of  this  combinatorial 
structure  in  any  number  of  ways.  Two  different  circle  packings,  obtained  using  the 
algorithm  in  Appendix  D  and  displayed  using  the  Mathematica  code  in  Appendix  E, 
are  shown  in  Figure  3.3.  In  each  of  Figure  3.3a  and  Figure  3.3b,  a  triangle  is  shaded 
for  reference;  the  triangle  in  each  case  corresponds  to  a  face  determined  by  the  same 
three  vertices  from  the  triangulation  shown  in  Figure  3.2. 

The  packings  Px  and  P2  in  Figure  3.3  are  embedded  in  C  and  are  normalized  by 
placing  the  center  of  a  circle  (Circle  1)  in  each  packing  at  the  origin  and  placing  the 
center  of  an  adjacent  circle  (Circle  2)  on  the  real  axis.  The  shaded  triangles  illustrate 
how  the  discrete  conformal  map  from  Px  to  P2  is  described.  The  map  which  takes 
the  shaded  triangle  in  3.3a  to  the  shaded  triangle  in  3.3b  is  given  by 

f(z)  =  f(x  +  iy)  =  1.67628a:  +  0.808753y.  (3.1.1) 

It  is  easy  to  verify  that  this  map  is  ft-quasiconformal  with  k  =  2.07267,  since  it 
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(a)  Euclidean  Carrier  for  Pj 


(b)  Euclidean  Carrier  for  P2 


Figure  3.3:  The  Carriers  of  Two  Different  Circle  Packings  of  K, 

maps  gross  circles,  not  just  infinitesimal  circles,  to  ellipses  with  k  =  2.07267  the  ratio 
between  their  major  and  minor  axes. 

Definition  3.1.5.  A  chain  of  circles  in  a  packing  P  for  an  abstract  triangulation 
1C  is  a  collection  of  circles  (CVl ,  CV2 , . . . ,  CVn)  in  P  such  that  Vi  and  vi+i  share  an  edge 
in  K,  for  i  =  l,2,...,n  —  1,  and  Vi  ^  Vj,  if  i  ^  j.  Thus,  a  chain  of  circles  describes 
a  non-self-intersecting  edge  path  in  K.  Similarly,  a  closed  chain  is  a  collection  of 
circles  corresponding  to  a  closed  non-self-intersecting  edge  path  in  1C. 


Figure  3.4:  Chains  of  Circles  and  Flowers 


Examples  of  closed  chains  of  circles  in  the  packing  P2  are  shown  (as  shaded  circles) 
in  Figure  3.4.  A  special  category  of  chains  leads  to  a  building-block  for  our  circle 
packings.  A  flower  consists  of  a  central  circle  and  some  number  of  petal  circles, 
the  closed  chain  of  successively  tangent  neighbors  surrounding  the  central  circle.  An 
example  of  a  flower  from  the  packing  P2  is  shown  (as  shaded  circles)  in  Figure  3.4b. 

One  particularly  useful  circle  packing  is  the  regular  hexagonal  packing,  a  pack¬ 
ing  in  which  every  interior  vertex  has  degree  6  and  each  circle  has  the  same  radius.  An 
example  of  a  regular  hexagonal  packing,  sometimes  also  referred  to  as  a  penny  pack¬ 
ing,  is  shown  in  Figure  3.5a.  Another  packing  commonly  seen  is  the  ball-bearing 
packing,  shown  in  Figure  3.5b. 


(a)  A  Regular  Hexagonal  Packing  (b)  A  Ball-Bearing  Packing 

Figure  3.5:  Two  Important  Types  of  Circle  Packings 


3.2  Discrete  Function  Theory 

The  important  characteristic  of  the  discrete  conformal  maps  induced  by  circle 
packings  is  not  that  they  are  quasiconformal;  the  key  fact,  suggested  by  Thurston  [51] 
and  proven  by  Rodin  and  Sullivan  [39],  is  that  these  maps  are  “nearly  conformal.” 
This  is  the  result  given  in  Theorem  3.2.1,  the  Rodin-Sullivan  Theorem.  Before  we 
state  this  theorem,  however,  we  first  state  some  geometric  results  associated  with 
circle  packing  that  are  interesting  in  themselves  and  required  for  the  proof  of  Theorem 
3.2.1. 
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Lemma  3.2.1  (Length- Area  Lemma).  Let  P  be  a  univalent  packing  in  ID)  and  Cv 
a  circle  in  P  with  Euclidian  radius  r.  Assume  there  exist  m  disjoint  chains  of  circles 
in  P  having  combinatorial  lengths  ni,n2, . . .  ,nm,  such  that  each  chain  separates  Cv 
from  0  and  a  point  on  dO.  Then 


r  < 


(3.2.1) 


i—  1 


As  the  number  of  generations  separating  a  circle  in  a  packing  from  the  boundary 
increases  Lemma  3.2.1,  the  Length- Area  Lemma,  has  the  effect  of  forcing  the  radius 
of  this  circle  to  zero  (in  the  limit).  For  a  more  detailed  discussion  of  Lemma  3.2.1, 
the  Length- Area  Lemma,  see  [39,  42]. 


Lemma  3.2.2  (Ring  Lemma).  Given  a  univalent  flower  ( Cvo ;  CVl ,  CV2 , . . . ,  CVn)  in 
C  there  is  a  lower  bound  Cn,  depending  only  on  n,  on  the  ratio  of  the  radius  rt  of  CVi 
to  the  radius  r0  of  CVo  for  each  i  =  1, 2, . . . ,  n;  that  is 

Cn  <  (3.2.2) 

r  o 

i  =  1,2, . . .  ,  n. 


Lemma  3.2.2,  the  Ring  Lemma,  guarantees  that  central  angles  in  the  carrier  on  a 
flower  are  bounded  away  from  zero  and  7 r.  That  is,  suppose  we  are  given  a  complex 
/C  in  which  the  degree  of  each  vertex,  the  number  of  adjacent  vertices,  is  bounded; 
also  suppose  we  have  two  (different)  packings  Pi  and  P2  associated  with  JC.  The  Ring 
Lemma  guarantees  that  the  quasiconformality  of  the  induced  conformal  map  from  Pi 
to  P2  is  bounded. 


Lemma  3.2.3  (Hexagonal  Packing  Lemma).  There  is  a  sequence  {sn}n  £  N, 
decreasing  to  zero,  with  the  following  property.  Let  Ci  be  a  circle  in  a  univalent 
Euclidean  circle  packing  P,  and  suppose  the  first  n  generations  of  circles  about  c\  are 
combinatorially  equivalent  to  n  generations  of  the  regular  hexagonal  packing  around 
one  of  its  circles.  Then  for  any  circle  c  G  P  tangent  to  C\, 


1  -  — 
ra 


(3.2.3) 
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where  rc  is  the  radius  of  the  circle  c  in  P  and  rCl  is  the  radius  of  the  circle  Ci  in  P. 

The  immediate  value  of  Lemma  3.2.1,  Lemma  3.2.2,  and  Lemma  3.2.3  is  their 
use  in  proving  Theorem  3.2.1,  the  Rodin-Sullivan  Theorem,  one  of  the  fundamental 
results  in  the  study  of  circle  packing. 

Theorem  3.2.1  (Rodin-Sullivan  Theorem).  Fix  a  simply  connected  domain  U  C 
C  and  points  p,q  6  fl.  Let  Pk  be  the  portion  lying  in  Ll  of  the  infinite  regular  hexagonal 
packing  whose  circles  all  have  radius  and  let  KLk  be  the  underlying  complex  for  the 
packing  Pk-  Suppose  Pk  is  a  packing  in  D  for  K,k  with  all  boundary  circles  tangent  to 
3D,  and  let  fk  :  carr(Pk)  — >  carr(Pk)  be  the  induced  discrete  conformal  map.  If  each 
Pk  has  been  normalized  so  that  fk(p)  =  0  and  fk{q)  >  0,  then  {fk}  converges  locally 
uniformly  to  the  unique  Riemann  map  f  :  Q  — >•  D  satisfying  f(p)  =  0  and  f(q)  >  0. 

An  illustration  of  the  use  of  circle  packing  to  approximate  the  action  of  the  unique 
Riemann  map  from  a  simply  connected  domain  to  the  disk  is  shown  in  Example  3.2.1. 

Example  3.2.1.  Let  =  {z  €  C  :  \z\  <  1  and  Re(z)  >  0}.  That  is,  Q  is  the  upper 
half  unit  disk.  We  wish  to  map  Q  to  the  unit  disk  D  with  the  unique  Riemann  map  /, 
normalized  so  that  f(i(V 2  —  1))  =  0  and  /  maps  the  imaginary  axis  to  the  imaginary 
axis.  In  this  case,  we  may  explicitly  construct  the  desired  map:  f(z)  =  —i  f  - 


Figure  3.6:  Unique  Riemann  Map  /  :  — >  D 

The  action  of  this  map  is  illustrated  in  Figure  3.6.  For  reference,  we  have  included  an 
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orthogonal  grid  on  fi,  shown  in  Figure  3.6a,  and  the  image  of  this  grid  in  D>,  shown 
in  Figure  3.6b. 


Figure  3.7:  Approximate  Riemann  Map  with  Circle  Packing  (171  Circles) 


(a)  Circles  in  ft 


(b)  Packing  in  the  Disk 


Figure  3.8:  Approximate  Riemann  Map  with  Circle  Packing  (342  Circles) 


(a)  Circles  in  ft 


(b)  Packing  in  the  Disk 


Figure  3.9:  Approximate  Riemann  Map  with  Circle  Packing  (684  Circles) 
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We  now  illustrate,  in  Figure  3.7,  Figure  3.8,  and  Figure  3.9,  the  approximation 
of  this  map  through  successively  finer  circle  packings,  showing  the  result  given  by 
Theorem  3.2.1,  the  Rodin-Sullivan  Theorem.  Clearly,  the  action  of  the  explicit  map 
described  in  Figure  3.6  is  emerging;  the  image  of  the  orthogonal  grid  in  Figure  3.6 
appears  in  chains  of  circles.  The  shaded  circle  in  each  of  Figure  3.7a,  Figure  3.8a, 
and  Figure  3.9a  represents  the  point  mapped  to  the  origin.  In  the  images  under  / 
shown  in  Figure  3.7b,  Figure  3.8b,  and  Figure  3.9b,  the  shaded  circles  are  centered 
at  the  origin. 

The  requirement  that  the  packings  used  in  Theorem  3.2.1,  the  Rodin-Sullivan 
Theorem,  all  be  of  uniformly  degree  6  is  quite  restrictive.  Since  the  initial  proof 
of  Theorem  3.2.1,  however,  Stephenson  [43]  relaxed  the  degree  6  condition  using 
techniques  of  random  walks,  and  He  and  Rodin  [27]  showed  that  only  a  uniform  bound 
on  the  degree  is  necessary.  To  thus  relax  the  requirement  on  the  combinatorics  of  the 
packing,  we  require  Lemma  3.2.4,  sometimes  referred  to  as  the  Packing  Lemma. 

Lemma  3.2.4  (Packing  Lemma).  Let  {Kn}n€n  be  a  sequence  of  combinatorial 
closed  disks  such  that 

1.  there  exists  a  uniform  bound  on  the  degree  of  the  vertices  in  Kn  for  each  n  E  N, 
and 

2.  the  sequence  {7f„}neN  is  either  a  nested  sequence  which  exhausts  a  parabolic 
combinatorial  disk  or  is  asymptotically  parabolic. 

There  exists  a  sequence  {sm}meN  C  M.  decreasing  to  zero,  with  the  following  property. 
Suppose  that  for  some  n,  u  and  v  are  adjacent  interior  vertices  of  Kn  whose  combi¬ 
natorial  distance  from  dKn  are  both  at  least  m,  and  suppose  that  Pn  and  Pn  are  two 
univalent,  Euclidean  circle  packings  for  Kn.  Then 

<  sm ,  (3.2.4) 

where  ru  and  rv  are  the  radii  of  the  circle  in  Pn  corresponding  to  u  and  v,  and  fu  and 
fv  are  the  radii  of  the  circle  in  Pn  corresponding  to  u  and  v. 
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Essentially,  Lemma  3.2.4,  the  Packing  Lemma,  states  that  for  a  circle  “deep”  in  a 
packing,  the  ratio  of  its  radius  to  any  given  neighbor  is  nearly  the  same  in  the  packings 
Pn  and  Pn;  in  other  words,  the  triangles  in  Pn  and  Pn  are  nearly  similar  triangles.  This 
gives  us  that  away  from  the  boundary,  the  induced  discrete  conformal  map  between 
Pn  and  Pn  is  nearly  conformal.  This  fact  will  play  an  important  role  in  the  proof  of 
our  main  result  in  Chapter  IV. 


3.3  Hex  Refinement 

In  order  to  obtain  the  various  approximation  results  for  circle  packing  and  dis¬ 
crete  analytic  function  theory,  we  need  a  method  to  refine  given  circle  packings.  The 
primary  requirement  in  any  such  refinement  is  that  we  maintain  some  uniform  con¬ 
trol  over  the  degree  of  the  complexes  generated  by  the  refinement  algorithm  since 
we  require  that  Lemma  3.2.2,  the  Ring  Lemma,  applies  at  each  successive  level  of 
refinement.  The  hex  refinement  method  developed  by  Bowers  and  Stephenson  [11]  is 
especially  nice. 

Definition  3.3.1.  If  1C  is  a  2-complex,  the  hex  refinement  of  K,  is  the  complex 
formed  by  adding  a  vertex  to  each  edge  and  adding  an  edge  between  any  two  vertices 
lying  on  the  same  face,  as  shown  in  Figure  3.10. 


Figure  3.10:  Hex  Refinement 

Note  that  hex  refinement,  and  refinement  in  general,  is  really  a  combinatorial  pro¬ 
cess,  refining  the  combinatorics  of  the  complex  K\  one  must  repack  the  new  complex 


51 


obtained  by  refining  K  in  order  to  realize  the  effect  of  the  refinement  in  a  circle  pack¬ 
ing.  Hex  refinement  has  a  number  useful  characteristics  described  by  Bowers  and 
Stephenson  [11]  and  summarized  in  Proposition  3.3.1. 

Proposition  3.3.1  (Bowers  and  Stephenson).  Any  new  interior  vertices  added  to 
a  2-complex  JC  by  hex  refinement  have  degree  6,  while  the  degree  of  the  original  vertices 
remain  unchanged.  If  fC  is  embedded  in  C  in  such  a  way  that  the  edges  correspond 
to  Euclidean  line  segments,  then  its  hex  refinement  may  be  realized  by  adding  line 
segments  joining  the  midpoint  of  each  edge  within  every  face.  In  this  case,  each  face 
in  1C  is  subdivided  into  four  new  faces,  each  similar  to  the  original  face  in  which  it  is 
contained  and  having  edges  one-half  as  long. 

Notice  that  refining  only  one  edge  in  a  complex  is  not  permitted,  since  this  would 
result  in  a  complex  that  is  not  a  triangulation;  the  faces  bordering  the  refined  face 
will  have  an  extra  vertex  along  the  common  edge  they  share  with  the  refined  face, 
giving  combinatorial  quadrilaterals  rather  than  triangles.  We  can,  however,  locally 
refine  only  those  triangles  in  the  complex  which  present  some  difficulty  with  respect 
to  desired  characteristics  of  the  complex,  then  correct  the  introduced  problems  on 
adjacent  faces  by  adding  a  single  edge  from  a  vertex  to  the  midpoint  of  the  opposite 
side,  as  shown  in  Figure  3.11. 


Figure  3.11:  Hex  Refinement  and  Correction  on  Adjacent  Triangle 

This  process  of  hex  refinement  of  individual  faces  (or  subsets  of  the  triangles  in 
a  triangulation)  and  the  addition  of  edges  to  absorb  extra  vertices  can  be  used  to 
locally  refine  an  abstract  triangulation  in  order  to  improve  the  discrete  conformal 
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approximation  in  troublesome  areas.  Local  refinement  will  play  a  key  role  in  the 
discrete  conformal  approximation  of  earthquakes  in  Chapter  IV. 

3.4  Combinatorial  Welding 

We  now  describe  the  process  of  combinatorial  welding  developed  by  Williams 
[53,  54].  Suppose  two  abstract  triangulations  K\  and  /C2  are  embedded  in  a  surface 
S  and  suppose  h  :  Bx  C  dfCx  — »  B2  C  d)C2  is  a  homeomorphism.  The  map  h  will 
be  used  to  attach  the  triangulations  1C  1  and  /C2  along  the  subsets  of  their  respective 
boundaries  Bx  and  B2  and  form  a  new  triangulation  K.  In  Chapter  IV,  h  will  be 
taken  as  a  discrete  version  of  the  map  defining  an  earthquake  along  a  geodesic,  and 
the  attaching  described  here  will  act  as  a  combinatorial  shearing  and/or  combinatorial 
grafting  in  the  construction  of  combinatorial  earthquakes. 

If  the  map  h  respects  the  combinatorial  structures  of  /C  1  and  /C2,  that  is,  if  h  sends 
vertices  and  edges  in  Bx  and  B2  to  vertices  and  edges,  respectively,  the  combinatorial 
welding  process  is  trivial.  We  simply  identify  vertices  and  edges  in  Bx  c  dfC  1  with 
their  images  in  B2  C  dlC2,  as  shown  in  Figure  3.12.  Note  that  the  action  of  h  on  the 
vertices  in  Bx  is  sufficient  to  determine  the  combinatorial  action  of  h  on  all  of  Bx. 


Figure  3.12:  Welding  Action  when  h  Respects  the  Combinatorics  of  K 1  and  fC2 

Much  more  interesting  is  the  case  in  which  h  does  not  respect  the  combinatorics  of 
1C  1  and  /C2  in  this  way.  In  general,  the  images  of  vertices  in  Bx  under  that  action  of  h 
will  not  be  vertices  in  B2.  In  order  that  the  map  h  be  well-defined  in  these  cases,  it  is 
necessary  to  modify  the  triangulations  K,  1  and  1C2  so  that  h  will  respect  the  modified 
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combinatorial  structure.  This  modification  is  accomplished  using  a  refinement  similar 
to  that  shown  in  Figure  3.11.  Specifically,  for  each  boundary  vertex  v  £  Bx  we  refine 
B2  by  adding  a  vertex  to  B2  embedded  at  the  point  h(v).  Likewise,  for  each  boundary 
vertex  w  E  B2  we  refine  B\  by  adding  a  vertex  to  B\  embedded  at  the  point  h~l  (w). 
When  vertices  are  thus  inserted,  however,  /Cx  and  JC2  will  not  necessarily  remain 
triangulations.  To  maintain  the  characteristics  of  a  triangulation,  we  must  further 
refine  each  of  K\  and  fC2.  If  we  add  a  vertex  h(v)  or  h~1(w)  to  an  edge  [a,  b]  of  a 
triangle  (a,  b ,  c),  we  refine  the  original  complex  (of  JC\  or  )C2)  by  adding  an  edge  from 
c  to  the  new  vertex.  An  example  of  this  process  is  shown  in  Figure  3.13. 


Figure  3.13:  Welding  Action  with  Combinatorial  Refinement  of  K\  and  JC2 

Note  that  this  refinement  process  is  well-defined  and  produces  new  triangulations 
since,  by  definition,  each  face  in  the  triangulations  may  have  at  most  one  edge  con¬ 
tained  in  the  boundary.  We  therefore  add  vertices  to  only  one  edge  in  any  given  face 
which  might  interfere  with  the  refinement  process. 

It  is  important  that  in  the  course  of  executing  these  combinatorial  weldings  and 
then  embedding  the  resulting  combinatorics  in  a  circle  packing  to  realize  a  geometry 
that  we  maintain  control  over  the  combinatorial  and  geometric  characteristics  of  the 
process.  Specifically,  we  must  ensure  that 

1.  the  degree  of  the  complex  resulting  from  the  combinatorial  welding  is  bounded, 
and 
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2.  the  new  edges  added  to  the  carriers  of  K, \  and  /C 2  do  not  result  in  triangles  with 
arbitrarily  small  angles. 

Condition  (1)  will  then  guarantee  that  we  can  apply  Lemma  3.2.2,  the  Ring 
Lemma,  to  the  packing  associated  with  the  new  complex  which  results  from  the  com¬ 
binatorial  welding,  thus  giving  bounds  on  the  angles  in  the  new  packing.  Condition 
(2)  will  ensure  the  existence  of  similar  bounds  for  the  packing  prior  to  welding  and 
repacking;  when  we  artificially  insert  edges  in  the  carrier  of  the  initial  packing  so  that 
there  is  a  bijection  between  triangles  in  the  two  carriers,  we  must  not  create  triangles 
with  arbitrarily  small  angles.  Thus,  with  these  two  conditions  satisfied  we  may  con¬ 
struct  quasiconformal  discrete  maps  between  the  carriers  in  the  surfaces  associated 
with  the  initial  packing  and  the  post-welding  packing.  These  technical  considera¬ 
tions  are  easily  overcome  through  the  use  of  local  hex  refinement  (Condition  (1))  and 
through  “rounding  off”  the  refinement  to  existing  vertices  if  the  new  vertex  added  in 
the  welding  refinement  is  too  close  to  an  existing  vertex  (Condition  (2)).  The  details 
of  these  processes,  and  verification  of  their  validity,  are  given  in  [53,  54].  Once  we  have 
welded  the  two  sides  of  the  complex  together  in  this  way,  we  may  pack  the  resulting 
complex.  Now,  we  have  a  bijection  between  the  triangles  in  the  carrier  of  the  original 
packing  (after  it  has  been  suitably  refined)  and  the  new  packing  which  results  once 
the  shearing  operation  is  completed;  we  use  these  triangles  to  construct  a  piecewise 
affine,  and  thus  quasiconformal,  map  from  the  original  packing  to  the  new  packing. 
Thus,  the  discretization  of  the  map  along  the  boundary  leads  to  a  discretization  of 
the  map  on  the  entire  surface.  This  now  leads  us  to  Proposition  3.4.1,  proven  in  [54]. 

Proposition  3.4.1  (Williams).  Combinatorial  attachment  by  h  induces  a  piecewise 
linearization  h  of  h.  If  h  is  bilipschitz,  then  the  linearization  h  will  be  bilipschitz  as 
well. 
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3.5  Density  of  Packable  Surfaces 

Demonstrations  of  the  existence  of  circle  packings  have  been  given  variously  by 
Thurston  [51],  Minda  and  Rodin  [36],  and  Beardon  and  Stephenson  [4].  Brooks  [12] 
showed  that  compact  packable  surfaces  are  dense  in  moduli  space,  and  Bowers  and 
Stephenson  [8,  10]  extended  this  result  to  include  surfaces  of  finite  analytic  type. 
Other  results  on  the  density  of  packable  surfaces  have  been  given  by  Barnard  and 
Williams  [1]  and  Williams  [54].  The  results  thus  demonstrated  which  are  germane  to 
this  research  are  summarized  in  Theorem  3.5.1. 

Theorem  3.5.1.  Let  K,  be  an  abstract  triangulation  of  a  surface  of  type  ( g,n,m ). 
There  exists  a  unique  surface  in  moduli  space  which  supports  a  packing  for  fC.  A  com¬ 
plex,  along  with  a  choice  of  marking,  then  determines  a  unique  point  in  Teichmuller 
space.  Moreover,  the  collection  of  all  packable  surfaces  is  dense  in  Teichmuller  space. 

The  density  of  merely  packable  surfaces  as  given  in  Theorem  3.5.1  is  extremely 
useful  and  powerful.  Of  equal  interest  is  the  density  of  other  subsets  of  packable 
surfaces  and  the  characteristics  of  convergent  sequences  of  packable  surfaces. 

Proposition  3.5.1.  Let  R  be  a  compact  Riemann  surface  that  is  not  packable.  If 
{Rn}™-_i  is  any  sequence  of  compact,  packable  Riemann  surfaces  such  that  Rn  —>  R 
in  the  Teichmuller  metric  where  each  admits  a  packing  Pn  with  degree  uniformly 
bounded  through  the  sequence,  then  the  radii  in  the  circle  packings  {Pn}™=  i  tend  to 
zero  as  n  —¥  oo. 

Proof.  First  note  that  since  the  circle  packings  {Pn}™=1  are  assumed  to  be  of  bounded 
degree,  Lemma  3.2.1  guarantees  that  the  radius  of  one  circle  in  a  packing  is  arbitrarily 
small  if  and  only  if  the  radii  of  all  circles  in  that  packing  are  small.  Proceeding  by  way 
of  contraposition,  suppose  that  the  radii  of  circles  in  Pn  are  bounded  away  from  zero 
by  some  fixed  positive  constant  for  each  neN.  For  every  n  £  N,  Rn  is  a  compact 
Riemann  surface,  and  therefore  has  finite  area.  Since  Rn  — »•  R  and  R  is  compact,  there 
exists  a  uniform  upper  bound  on  the  area  of  the  Riemann  surfaces  Rn.  Hence,  there 
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exists  a  uniform  bound  on  the  total  number  of  circles  in  each  packing  Pn.  There  are 
thus  only  a  finite  number  of  abstract  triangulations  possible  among  the  circle  packings 
{Pn}^Li-  By  the  pigeonhole  principle,  there  exists  a  subsequence  {-R^l^Lx  C  {Rn)^i 
such  that  for  each  k  £  N  the  circle  packing  realized  in  Pnk  corresponds  to  a  single 
abstract  triangulation  JC,  (i.e.,  there  exists  a  subsequence  of  packings  each  having  the 
same  combinatorial  structure).  Since  the  abstract  triangulation  uniquely  determines 
a  surface  for  which  that  triangulation  is  realized  as  a  packing,  we  find  that  {Rnk}(f=] 
is  a  constant  sequence  of  Riemann  surfaces.  Now,  Rn  -4  R  as  n  -4  oo,  so  every 
subsequence  of  {Rn}((Li  must  likewise  converge  to  R;  that  is  Rnk  -4  R  as  k  — >  oo. 
Thus  R  is  an  element  of  the  sequence  {-RnK^Lii  in  fact,  R  is  equivalent  to  infinitely 
many  elements  of  this  sequence.  Therefore,  R  is  a  compact,  packable  Riemann  surface. 
Thus,  by  our  contrapositive  argument,  if  a  sequence  of  packable  surfaces,  Rn,  with 
degree  uniformly  bounded  through  the  sequence  converge  to  a  Riemann  surface  that 
is  not,  itself,  packable,  then  the  mesh  of  the  packings  on  the  surfaces  goes  to  zero.  □ 

We  now  introduce  another  class  of  compact  Riemann  surfaces  called  equilateral 
surfaces  which  will  be  useful  in  constructing  sequences  of  packable  Riemann  surfaces 
while  maintaining  control  over  the  combinatorics  of  their  underlying  complexes. 

Definition  3.5.1.  Suppose  S  denotes  a  compact,  orientable  topological  surface  and 
let  K  denote  a  triangulation  of  the  surface  S.  If  we  paste  together  equilateral  triangles 
(triangles  conformally  equivalent  to  equilateral  triangles)  in  the  pattern  of  K  to  impose 
a  piecewise  affine  structure  on  S,  this  affine  structure  defines  a  conformal  structure 
on  S,  guaranteeing  that  S  is  a  Riemann  surface.  Riemann  surfaces  thus  constructed 
are  called  equilateral  surfaces. 

Just  as  has  been  shown  for  compact  surfaces  by  Brooks  [12]  and  for  surfaces  of  finite 
analytic  type  by  Bowers  and  Stephenson  [8,  10],  if  we  fix  a  genus  g  >  0  this  seemingly 
very  restrictive  class  of  (compact)  equilateral  surfaces  of  genus  g  is  dense  in  Teichmiiler 
space  of  genus  g,  as  shown  by  Belyi  [3].  We  state  this  result  as  Theorem  3.5.2. 
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Theorem  3.5.2  (Belyi).  If  S  is  a  Riemann  surface  of  genus  g  >  0,  the  set  of 
equilateral  surfaces  of  genus  g  is  countable  and  dense  in  the  Teichmuller  space  of  S. 

This  result  guarantees  the  existence  of  sequences  of  packable  Riemann  surfaces  with 
underlying  combinatorial  structures  over  which  we  may  exercise  a  significant  degree 
of  control.  This  ability  to  control  the  underlying  combinatorics  allows  us  the  ability 
to  manipulate  the  geometry  of  our  packings  in  that  we  may  force  the  circles  in  a 
sequence  of  packable  (and  packed)  surfaces  to  decrease  in  size. 

Corollary  3.5.1.  Let  S  be  a  Riemann  surface  and  let  R  be  an  arbitrary  point  in  the 
Teichmuller  space  of  S.  There  exists  a  sequence  of  points  {R.n}  in  the  Teichmuller 
space  of  S  such  that  Rn  —¥  R  in  the  Teichmuller  metric  as  n  — *  oo,  Rn  is  packable 
for  every  n,  and  the  radii  of  the  circles  in  Pn  go  to  zero  as  n  — >  oo,  where  Pn  is  the 
unique  packing  on  the  surface  Rn  for  every  n. 

Proof.  First  note  that  since  the  set  of  equilateral  surfaces  on  S  is  dense  in  the  Te¬ 
ichmuller  space  of  S,  there  exists  a  sequence  of  equilateral  surfaces  { E *}  in  the  Te- 
ichmiiller  space  of  S  such  that  — »•  R  as  n  — >  oo.  Corresponding  to  the  equilateral 
surface  E\  we  have  a  triangulation  K\.  This  triangulation  corresponds  to  a  unique 
packable  surface  R{,  distinct  from  E\.  An  example  of  an  equilateral  surface  and  the 
packable  surface  it  induces  is  shown  in  Figure  3.14.  (This  figure  appeared  in  [11]  and 
is  used  with  permission  of  the  authors.) 

Now,  refine  the  triangulation  K\  using  hex  refinement  to  create  a  new  triangula¬ 
tion  Kf.  Note  that  refining  this  triangulation  has  no  effect  on  the  structure  of  the 
surface  E\\  that  is,  the  equilateral  surface  E\  corresponding  to  Kf  is  the  same  as  E\ , 
since  hex  refinement  on  an  equilateral  triangle  divides  that  face  into  4  new  equilat¬ 
eral  triangles.  So,  K\  corresponds  to  a  unique  packable  surface  R\,  distinct  from  E\ 
(and  El).  Continuing  in  this  manner,  we  construct  a  sequence  of  packable  surfaces 
{RT)m=i-  As  m  — >  oo,  the  triangles  in  the  carrier  of  the  packing  P{n  on  R™  converge 
to  equilateral  triangles,  and  the  sequence  {R™}^=1  converges  to  a  surface  confor¬ 
mally  equivalent  to  E\.  Thus,  R™  — »■  E\  as  m  — >  oo  in  the  Teicmuller  metric.  Note, 
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Figure  3.14:  An  Equilateral  Surface  and  the  Induced  Packable  Surface 


however,  that  the  degree  of  each  packing  in  the  sequence  is  bounded  by 

the  same  constant,  since  the  hex  packing  does  not  change  the  degree  of  any  vertex  in 
the  initial  packing  and  every  new  vertex  has  degree  6  since  all  vertices  in  the  packing 
are  interior  vertices.  Now,  a  standard  application  of  Lemma  3.2.1,  the  Length-Area 
Lemma,  guarantees  that  the  radii  of  the  circles  in  the  sequence  {P™}^=i  go  to  zero  as 
k  -*•  oo.  We  repeat  this  process  for  each  equilateral  surface  P*,  generating  for  each  a 
sequence  of  packable  surfaces  {P™}*=1  such  that  R™  — *■  P*  in  the  Teichmiiller  metric 
and  the  radii  of  the  circles  in  the  sequence  of  packings  {P™}™=1  corresponding  the 
packable  surfaces  go  to  zero  as  m  — >•  oo. 

R\  Rl  . 

P?  PI  . 

El  E\  •••  R 

Now,  for  each  n  €  N,  choose  a  surface  Rn  from  the  sequence  {P”!}^=1  such  that 
the  Teichmiiller  distance  from  Rn  to  P*  is  less  than  2~n  and  the  maximum  radius  of 
any  circle  in  the  packing  on  Rn  (in  the  intrinsic  metric  on  the  surface)  is  similarly  less 
than  2~n.  We  now  have  a  sequence  {Pn}^l;1  such  that  Rn  R  in  the  Teichmiiller 


metric  as  n  -4  oo,  Rn  is  packable  for  every  n,  and  the  radii  of  the  circles  in  Pk  go  to 
zero  as  n  -4  oo,  where  Pn  is  the  unique  packing  on  the  surface  Rn  for  every  n.  □ 
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CHAPTER  IV 
PACKING  EARTHQUAKES 


Throughout  this  chapter  let  (£,  a)  be  a  finite  measured  geodesic  lamination  of 
D  associated  with  a  set  of  points  S  C  ID.  Suppose  C  consists  of  n  disjoint,  simple, 
closed  geodesics  {L®}"=1.  Note  that  throughout  this  chapter  we  suppose  that  the  real 
(shearing)  part  of  the  weight  on  any  geodesic  Ll  may  assume  any  real  value,  while  the 
imaginary  (grafting)  part  may  assume  only  non-negative  real  values.  This  indicates 
that  the  shearing  maps  we  consider  may  shear  to  either  the  left  or  the  right,  and  the 
grafting  maps  we  consider  involve  only  the  insertion  of  cylinders  of  positive  height. 
Much  of  the  construction  below  is  taken  from  [54];  we  reproduce  it  here  in  the  interest 
of  making  this  exposition  as  self-contained  as  possible. 

4.1  Discrete  Laminations 

Let  {Pn}kLi  be  a  sequence  of  packings  in  D  with  uniformly  bounded  degree  and 
mesh  decreasing  to  zero.  Let  Dk  denote  the  Euclidean  carrier  of  Pk  for  each  k  €  N. 
Note  that  the  construction  of  a  “combinatorial”  earthquake  Ek  of  this  Euclidean 
carrier  as  an  approximation  of  the  finite  earthquake  E  is  complicated  by  the  fact  that 
Dk  is  not  equal  to  D  for  any  fixed  k  G  N. 

Note  that  any  given  geodesic  in  the  lamination  C  will  pass  indiscriminately  through 
the  points  and  triangles  of  the  Euclidean  carrier  Dk.  Thus,  shearing  and  grafting 
Dk  along  the  geodesics  of  the  lamination  is  not  directly  possible.  In  order  to  ap¬ 
proximate  the  action  of  the  earthquake  E,  therefore,  we  construct  a  combinatorial 
lamination  of  Dk  that  respects  the  combinatorial  structure  and  converges  to  the  finite 
geodesic  lamination  £  as  k  — >■  oo.  This  combinatorial  lamination  comprises  combi¬ 
natorial  “geodesics”  {£k}™=l  corresponding  to  the  true  hyperbolic  geodesics  (L*}”=1, 
constructed  in  the  following  manner. 

We  consider  packings  that  are  hyperbolic  in  D;  that  is,  we  consider  packings  that 
fill  the  disk.  Let  a\  and  b\  be  the  endpoint  of  the  geodesic  Ll  C  3D.  We  construct  t\ 
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as  an  edge  path  which,  in  the  limit,  connects  a\  and  b\,  whose  Hausdorff  distance  to 
L1  is  minimal,  and  which  satisfies  the  following  three  conditions. 

1.  Every  vertex  in  l\  must  have  a  neighbor  in  both  components  of  Dk\tk. 

2.  The  combinatorial  geodesics  {£lk}  are  disjoint,  except  possibly  at  their  common 
endpoints. 

3.  Any  hyperbolic  geodesic  perpendicular  to  Ll  intersects  i\  in  exactly  one  point. 

It  is  not  immediately  clear  that  we  can  in  general  construct  such  a  rigidly  character¬ 
ized  edge  path  through  the  combinatorics  of  Dk.  Lemma  4.1.1  guarantees  that  such 
constructions  are,  in  fact,  possible. 


Lemma  4.1.1  (Williams).  By  taking  k  sufficiently  large  and  possibly  modifying 
Dn  near  L\  uie  can  find  paths  t\  so  that  conditions  (1)  through  (3)  are  satisfied. 
Moreover,  i\  — >  L1  uniformly  as  k  — >  oo. 


In  Figure  4.1  we  see  successively  finer  packings  in  the  disk  and  increasingly  more 
accurate  piecewise  linear  combinatorial  approximations  of  the  hyperbolic  geodesic  L. 


(a)  Combinatorial  Geodesic 
with  19  Circles 


(b)  Combinatorial  Geodesic 
with  127  Circles 


(c)  Combinatorial  Geodesic 
with  547  Circles 


Figure  4.1:  Approximating  a  Geodesic  L  with  Combinatorial  Geodesics  £k 


62 


4.2  Hyperbolic  Projections 

Having  described  a  discrete  analog  of  the  finite  measured  geodesic  lamination 
(£,  a),  we  now  must  describe  discrete  analogs  for  the  hyperbolic  shearing  maps  and 
grafting  maps  associated  with  each  geodesic.  We  cannot,  of  course,  use  the  shearing 
and  grafting  maps  associated  with  the  lamination  (£,  a)  directly;  rather,  we  construct 
hyperbolic  projections  p\  which  take  the  combinatorial  geodesic  i\  to  L\  Thus, 
by  pre-and  post-composition  with  the  map  p\,  we  may  construct  an  appropriate 
transformation  to  give  a  discrete  shearing  and/or  grafting  along  l\  using  the  shearing 
and/or  grafting  map  on  L*. 

Notice  that  the  collection  of  hyperbolic  geodesics  perpendicular  to  the  geodesic 
Ll  in  the  lamination  fill  O.  By  the  construction  of  the  discrete  geodesic  £*k,  each 
such  perpendicular  geodesic  p  intersects  i\  exactly  once.  We  may  thus  define  a  map 
p\  :  Dk  fl  D>  — »■  D  by  the  requirement  that  p\\n  be  the  unique  hyperbolic  Mobius 
transformation  with  axis  p  and  translation  length  equal  to  the  hyperbolic  distance 
from  t\  fl  p  to  Ll  n  p.  The  action  of  such  a  map  on  a  segment  of  a  discrete  geodesic 
in  D  corresponding  to  the  geodesic  between  1  and  i  is  shown  in  Figure  4.2.  Note  that 
this  map  is  clearly  a  quasisymmetry. 


Figure  4.2:  Action  of  a  Hyperbolic  Projection  plk 

We  have  constructed  discrete  geodesic  approximations  llk  to  the  geodesic  L1  and 
the  corresponding  hyperbolic  projections  p\,  and  Lemma  4.1.1  guarantees  that  i\  — ¥ 
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Ll  as  k  — *  oo.  We  now  ask  how  the  hyperbolic  projections  p],  behave  as  k  oo.  The 
answer  is  given  by  Proposition  4.2.1. 

Proposition  4.2.1  (Williams).  The  sequence  {pk}  is  uniformly  bilipschitz  on  B 
and  converges  uniformly  to  the  identity  map  on  D. 

4.3  Discrete  Shearing  Maps 

With  the  hyperbolic  projections  {pk}  we  now  define  discrete  shearing  operations 
{h\}  which  approximate  the  hyperbolic  shearing  maps  {h1}  such  that  h'k  — >■  hl  as 
k  —t  oo.  For  each  k  €  N,  define  h\  by 

K  =  (Pk)~lohio(pi).  (4.3.1) 


mm 


(a)  A  Combinatorial  Geodesic  (b)  Local  Refinement,  Shearing,  and 

Welding 

Figure  4.3:  Combinatorial  Shearing  and  Welding 

By  Lemma  4.3.1,  this  map  behaves  in  the  manner  we  require,  simulating  the  action 
of  a  hyperbolic  shearing  and  converging  to  that  hyperbolic  shearing. 

Proposition  4.3.1  (Williams).  The  discretized  maps  h\  are  bilipschitz  on  tk  with 
bilipschitz  constants  bounded  independently  of  i  and  k,  and  hk  converges  uniformly  to 
the  hyperbolic  shearing  map  hl  onD  as  k  — »  oo. 

4.4  Discrete  Grafting  Maps 

As  we  did  for  discrete  shearing  maps  in  Section  4.3,  we  may  use  the  hyperbolic 
projections  {pk}  to  define  discrete  grafting  operations  {gk}  which  approximate  the 
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hyperbolic  grafting  maps  {g1}  such  that  g\  — >  gl  as  k  — >  oo.  First,  we  need  to 
understand  the  action  of  the  grafting  map,  gl.  The  action  of  this  map  is  to  open  the 
disk  along  the  geodesic  L\  creating  two  copies  of  the  geodesic.  We  then  parameterize 
each  copy  of  the  geodesic  and  the  boundary  of  an  infinite  hyperbolic  strip  by  arc  length 
and  identify  corresponding  points  in  the  parameterizations.  The  only  difference  in 
the  discrete  grafting  is  that  the  discrete  geodesic  tk  and  Ll  do  not  agree.  So,  for  each 
fceN,  define  g\  by 

d  =  (pir1  °  s‘ °  (pi)-  (4.4.1) 


(a)  A  Combinatorial  Geodesic 


(b)  Local  Refinement,  Grafting,  and 
Welding 


Figure  4.4:  Combinatorial  Grafting  and  Welding 


Proposition  4.4.1.  The  discretized  maps  g\  are  bilipschitz  on  i\  with  bilipschitz  con¬ 
stants  bounded  independently  ofi  and  k,  and  g\  converges  uniformly  to  the  hyperbolic 
grafting  map  gl  on  O  as  k  — >•  oo. 

Proof.  The  fact  that  the  maps  g\  are  bilipschitz  on  follows  directly  from  Proposition 
3.4.1.  As  for  convergence,  Proposition  4.2.1  guarantees  that  p\  converges  uniformly 
to  the  identity.  Thus,  g\  converges  uniformly  to  gl  as  k  — >•  oo.  □ 

We  now  define  a  combinatorial  earthquake  as  the  composition  of  a  discrete  grafting 
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map  with  a  discrete  shearing  map,  Ek  =  gk  o  hk,  just  as  we  defined  continuous 
earthquakes  as  the  composition  of  grafting  and  shearing  maps,  E  =  g  o  h. 

4.5  Convergence  of  Discrete  Earthquake  Maps  on  Hyperbolic 

Riemann  Surfaces 

Thus  far,  we  have  described  a  family  of  transformations  on  the  Teichmiiller  space 
of  a  hyperbolic  Riemann  surface  and  a  process  whereby  we  may  discretely  mimic  the 
action  of  those  transformations  through  circle  packing.  In  Theorem  4.5.1  we  show 
that  this  ’’mimicry”  is  in  fact  a  convergent  numerical  method. 

Theorem  4.5.1.  Let  R  be  a  compact  hyperbolic  Riemann  surface  with  a  finite  mea¬ 
sured  geodesic  lamination  (£,  <r).  Let  {Pk}  be  a  sequence  of  finite  bounded  degree 
packings  with  mesh  decreasing  to  zero  corresponding  to  Riemann  surfaces  {Rk}  such 
that  Rk  — >  R  as  k  — >•  oo  in  the  Teichmiiller  metric.  Then  the  surfaces  Ek(Rk )  =  Rk 
induced  by  the  discrete  earthquake  maps  Ek  converge  to  the  surface  E{R)  =  R  under 
the  earthquake  map  E  induced  by  (£,  a) . 

Proof.  Let  R  be  a  compact  hyperbolic  Riemann  surface,  and  let  (£,  a)  be  a  finite 
measured  geodesic  lamination  on  R.  In  a  sense,  the  proof  of  this  result  reduces  to 
showing  that  the  diagram  (4.5.1)  commutes. 


{i?,  Cl,  C2,  . . . 

E 

>  Cm  /  5 

•  R 

fk\fik 

fk 

(4.5.1) 

{Rki  cl*.j  c2fc>  •  • 

■,CmJ 

"  Rk 

Note  that  requiring  the  mesh  of  the  sequence  of  bounded  degree  packings  {Pk}  to 
approach  zero  as  k  — >  oo  is  equivalent  to  requiring  that  the  radii  of  the  circles  in  the 
packings  approach  zero.  Thus,  we  are  guaranteed  the  existence  of  such  a  sequence 
by  Corollary  3.5.1.  Further,  the  construction  of  the  surfaces  in  the  proof  of  Corollary 
3.5.1  allows  the  freedom  to  control  the  degree  of  vertices  adjacent  to  the  geodesics 
along  which  we  shear  and  graft. 
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The  Teichmiiller  distance  between  the  packable  surfaces  Rk  and  R  is  going  to  zero, 
so  for  each  k  E  N  there  exists  a  map  :  Rk  — >  R  that  respects  the  markings  on 
Rk  and  R  and  is  1  +  e^-quasiconformal,  where  ek  — »  0  as  k  — ¥  oo,  (i.e.,  as  k  — >■  oo, 
/*—>/,  where  /  is  conformal).  More  precisely,  the  map  fk  consists  of  a  map  from 
the  disk  (or  other  model  of  the  hyperbolic  plane)  to  the  disk  through  a  map  between 
the  surfaces  and  their  respective  conformal  structures.  Where  no  ambiguity  arises, 
however,  we  will  speak  about  the  map  as  fk'-  Rk — >■  R- 

Let  {cj}™=l,  be  the  cylinders  grafted  into  the  Riemann  surface  R  by  the  finite 
earthquake  E.  For  each  j  =  1,2, ...,m,  let  {cJi}fceN  be  a  sequence  of  cylinders 
admitting  a  circle  packing  of  finite,  uniformly  bounded  degree  so  that  c3k  — >  Cj  as 
k  — »  oo  in  the  Teichmiiller  metric.  Further,  require  that  the  mesh  of  the  packing 
on  Cjk  goes  to  zero  as  k  -¥  oo  for  each  j  =  1, 2, . . . ,  m.  Thus,  for  each  k  £  N  there 
exists  a  map  fjk  :  Cjk  — »  Cj  that  is  1  +  £jk -qu asiconform al ,  where  £jk  — >  0  as  k  — >  oo, 
j  =  1, 2, . . . ,  m,  (i.e.,  as  k  -*  oo,  fjk  — v  fj,  where  f3  is  conformal). 

Let  R  be  the  image  of  R  under  the  action  of  the  earthquake  E,  and  for  each  k  eM 
let  Rk  be  the  image  of  Rk  under  the  action  of  a  combinatorial  earthquake  Ek.  We 
construct  these  combinatorial  earthquakes  by  describing  laminations  on  the  surfaces 
Rk  corresponding  to  the  lamination  on  R.  Since  on  a  hyperbolic  Riemann  surface 
there  is  exactly  one  simple  closed  geodesic  in  each  non-zero  homotopy  class  [14,  29], 
corresponding  to  each  geodesic  in  the  lamination  £  on  R  there  is  a  unique  geodesic 
on  Rk  in  the  corresponding  homotopy  class.  Taking  the  collection  of  these  geodesics, 
we  have  a  finite  measured  geodesic  lamination  (£*,  a)  on  Rk.  Notice  that  we  have  left 
the  weights  a  unchanged;  we  simply  assign  to  each  geodesic  in  £fc  the  weight  on  the 
corresponding  geodesic  in  R.  Now,  we  have  finite  measured  geodesic  laminations  on 
each  Rk,  k  €  N,  and  we  may  carry  out  the  discrete  shearing  and  welding  operations 
of  combinatorial  earthquakes. 

Around  each  geodesic  Ll,i  =  1,2, ...,n,  in  the  lamination  £  on  R  we  place  a 
collar,  a  neighborhood  of  uniform  width  isometric  to  a  cylinder,  so  that  no  two  such 
neighborhoods  intersect.  Call  this  collar  C\  Since  the  geodesics  in  the  finite  lamina- 
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tion  are  disjoint  by  construction,  the  Collar  Theorem  [14]  guarantees  that  there  exists 
such  a  collar  about  each  geodesic.  Similarly,  about  each  geodesic  Lk,i  =  1,2, ...  ,n 
in  the  lamination  Ck  on  Rk  there  exists  a  collar  Ck  such  that  the  collection  of  collars 
on  Rk  is  disjoint.  We  require  that  the  collar  Clk  on  Rk  have  width  less  than  \  the 
width  of  the  corresponding  Cl  on  R.  Now,  since  Rk  — >■  R  in  the  Teichmiiller  metric, 
the  maps  between  Rk  and  R  are  going  to  conformal  maps  and  the  markings  on  Rk 
are  converging  to  the  markings  on  R.  Thus,  for  k  sufficiently  large,  the  image  of 
the  collar  Ck  under  the  map  fk  is  contained  in  the  collar  C\  and,  by  construction, 
this  containment  is  proper,  (i.e.,  fk(Ck)  £  Cl  for  all  k  sufficiently  large).  Further, 
since  the  mesh  of  the  packings  on  Rk  goes  to  zero,  for  k  sufficiently  large,  the  discrete 
geodesic  i\  corresponding  to  Llk,i  =  1, 2, . . . ,  n,  lies  strictly  within  the  interior  of  the 
collar  Ck. 

Consider  the  images  of  the  collars  Ck  under  the  action  of  the  discrete  earthquake 
Ek.  If  the  geodesic  Lk  about  which  we  built  the  collar  Ck  has  real  weight,  and 
thus  induces  only  a  discrete  shearing  action  on  the  discrete  geodesic,  then  there  is 
no  difficulty  in  defining  what  we  will  call  a  quasicollar  on  Rk  as  the  image  of  C{ 
under  Ek.  The  image  of  the  discrete  geodesic  tk  and  the  image  of  the  geodesic  Lk 
are  contained  in  this  quasicollar.  This  is  shown  in  Figure  4.5,  where  the  geodesic 
Vk  is  shown  as  a  solid  gray  curve,  the  discrete  geodesic  i\  is  shown  as  a  solid  black, 
piecewise  linear  curve,  and  the  collar  C\  is  denoted  by  dashed  curves. 

The  case  in  which  the  weight  on  L\  has  non-zero  imaginary  part  is  somewhat 
more  problematic.  The  action  of  the  earthquake  in  this  case  will  divide  the  collar 
into  two  disjoint  sets.  Since  the  grafting  (imaginary)  component  of  the  earthquake 
involves  the  insertion  of  a  non-trivial  cylinder,  however,  we  may  extend  these  disjoint 
regions  at  least  one  generation  of  triangles  deep  into  the  interior  of  the  cylinder  to 
create  quasicollars  about  the  ends  of  the  inserted  cylinder.  (Note  that  since  the  mesh 
of  the  packing  is  going  to  zero,  the  region  occupied  by  one  generation  of  triangles  is 
getting  arbitrarily  small.)  Since  the  inserted  cylinder  has  non-zero  height,  we  may 
construct  these  quasicollars  so  that  they  are  disjoint. 
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(a)  Geodesic,  Discrete  Geodesic,  and  Collar  on  Rk 


(b)  Images  under  a  Discrete  Shearing  Map 


Figure  4.5:  Geodesics,  Discrete  Geodesics,  and  Quasicollars 


Define  a  pullback  map  Pk  :  Rk  ->•  {Rk,  clfc, . . . ,  cmk }  by  letting  Pk(r)  be  the  unique 
pre-image  r  e  Rk  of  the  point  f  G  Rk  under  the  combinatorial  earthquake  Ek.  Note 

-X*S. 

that  since  the  degree  of  the  vertices  in  the  complex  on  Rk  is  bounded  and  since,  by 
construction,  the  angles  in  the  packing  on  Rk  are  bounded  away  from  zero,  the  map 
Pk  is  1  +  e/c-quasicon  formal.  Further,  at  any  point  isolated  from  the  combinatorial 
geodesics,  as  k  — ¥  oo  the  number  of  generations  of  circles  on  the  packed  surfaces 
between  that  point  and  the  combinatorial  geodesics  goes  to  infinity;  thus,  at  any  such 
point,  Lemma  3.2.4,  the  Packing  Lemma,  guarantees  that  the  dilatation  of  the  maps 
Pk  goes  to  1  as  k  — >■  oo,  (i.e.,  as  k  — >  oo,  on  points  isolated  from  the  combinatorial 

xv  xs 

geodesics,  Pk  — >  P,  where  P  is  conformal). 

Since  each  of  R,  Rk,  and  Rk  are  Riemann  surfaces,  they  each  have  a  conformal 

xs 

structure;  in  the  case  of  Rk  and  Rk,  they  inherit  their  conformal  structures  from 
their  packings.  Let  {<pv}ve T,  Wkv}verk>  and  be  the  collection 

xv  xv 

of  coordinate  charts  in  the  conformal  structure  on  R,  R,  Rk,  and  Rk,  respectively. 

xv 

Note  that  the  structure  on  R  will  be  the  structure  described  in  Section  2.6.  Now, 
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we  consider  a  collection  of  maps  fk'-Rk — >•  R,  each  defined  on  open  sets  in  Rk.  Let 
U  c  Rk. 

Case  I:  Suppose  U  C  Rk  is  an  open  set  such  that  U  does  not  intersect  the  interior 

of  any  quasicollar,  and  U  does  not  intersect  any  cylinder  inserted  by  the  earthquake 

— «. 

Ek .  Associated  with  this  open  set  in  Rk  we  have  a  coordinate  chart  ip kv-  Now,  the 
map  P~l  o  fk  o  Pk  o  ip~^  takes  an  open  subset  of  the  unit  disk  B  corresponding  to 
the  pair  (U,ipkv)  to  an  open  subset  Ua  C  R  so  that  Ua  G  Ua,  as  described  in  Section 
2.6.  Corresponding  to  Ua,  we  take  from  the  conformal  structure  on  R  described  in 
Section  2.6  a  coordinate  chart  ipa.  Thus  we  have  the  following: 

ipa  o  P~l  O  fk  o  Pk  o  ip~'  =  (fiaoPoP~l  o/fcOPfcO-0-1 

=  Va°fk  oPfcO-0-1. 

Since  <pa  and  ipjp*  are  conformal  (or  1-quasiconformal)  and  the  maps  fk  and  Pk  are 
(1  +  £fc)-quasiconformal  and  (1  +  ik ) -quasiconforrn al ,  respectively,  the  map  thus  de¬ 
scribed  is  (l+£fc)(l+£0)-quasiconformal,  by  Lemma  2.4.1.  Note  that  sk  and  ik  are  in¬ 
dependent  of  the  set  U,  since  they  depend  only  on  the  maps  fk  and  Pk,  which  are  glob¬ 
ally  (1  -I-  £fc)-quasiconformal  and  (1  +  efc)-quasiconformal,  respectively.  By  construc¬ 
tion,  U  is  isolated  in  Rk  from  the  images  of  the  geodesics  defining  the  earthquake  Ek, 
and  £fe,  ik  — >■  0  aa  k  — >  oo.  Thus,  as  k  — »  oo,  the  map  0Q° P-1  °  fk° Pk° '■  Rk  R 
approaches  a  conformal  map. 

Case  II:  Suppose  U  C  Rk  is  an  open  set  such  that  U  does  not  intersect  the  interior 
of  any  quasicollar,  and  U  is  contained  in  the  interior  of  some  cylinder  inserted  by  the 
earthquake  Ek .  Associated  with  this  open  set  in  Rk  we  have  a  coordinate  chart  ip 
Now,  the  map  P_1  ofjk  °Pkoip^  takes  an  open  subset  of  the  unit  disk  B  corresponding 
to  the  coordinate  pair  (U,  ipkv)  from  the  conformal  structure  on  Rk  to  an  open  subset 

/S 

Up  C  R  so  that  Up  €  Ubi  as  described  in  Section  2.6.  Corresponding  to  Up,  we  take 
from  the  conformal  structure  on  R  described  in  Section  2.6  a  coordinate  chart  ipp. 
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Thus  we  have  the  following: 


$PoP  1  °  fhoPk°^K  =  Mi  o  Si  °  </>p  o  P  o  P  1  °  fh  °pk° 

=  Mj  0  9p  °  <f>p  °  fjk  °pk°  'ipk!  ■ 

Since  Mj  o  ^  o  ^  and  are  conformal  (or  1-quasiconformal)  and  the  maps  fJk 
and  Pk  are  (1  +  e^J-quasiconformal  and  (1  +  efc)-quasiconformal,  respectively,  the 
map  thus  described  is  (1  +  £jk){  1  +  e/k)-quasiconformal.  Again,  we  note  that  the 
constants  describing  the  deviation  from  conformality,  £jk  and  ik,  are  independent  of 
the  choice  of  U.  By  construction,  U  is  isolated  in  Rk  from  the  images  of  the  geodesics 
defining  the  earthquake  Ek,  and  £jk,£k  — >  0  as  k  — »  oo.  Thus,  as  k  — »  oo,  the  map 
ipp  o  P~L  o  fjk  oPko  (pkj  :  Rk  — ^  R  approaches  a  conformal  map. 

Case  III:  Suppose  that  U  is  a  small  open  subset  of  Rk  such  that  U  intersects 
Ek(£k)  for  some  i  =  1,2, ...  ,n,  and  U  intersects  the  boundary  of  the  quasicollar  Ck 

associated  with  Lk.  Corresponding  to  U  we  have  a  coordinate  chart  tpkl/.  Now,  we 

have  two  cases  to  consider  depending  on  whether  the  weight  on  the  geodesic  hl  has 
non-zero  imaginary  part. 

Case  Ilia:  Suppose  Im(a(L*))  =  0;  that  is,  the  earthquake  action  on  the  geodesic 
U  involves  only  a  shearing  action.  This  gives  that  the  map  Pk  takes  U  to  a  subset  of 
Rk  such  that  Pk{U)  ^  Ck  and  Pk(U)f]tk  /  0,  as  shown  in  Figure  4.6. 


(a)  An  Open  Set  U  C  Rk  (b)  Image  of  U  Under  Pk 

Figure  4.6:  Pullback  of  Geodesics,  Discrete  Geodesics,  and  Quasicollars 
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Note  that  in  Figure  4.6b  we  have  included  the  refining  edges  associated  with  the 
shearing  so  that  there  is  a  bijection  between  individual  triangles  in  Figure  4.6a  and 
Figure  4.6b. 

Now,  we  wish  to  send  this  set  Pk(U)  to  a  subset  of  R,  but  there  is  a  geometric 
difficulty  to  be  resolved.  The  “halves”  of  Pk(U )  have  as  one  portion  of  their  boundaries 
a  segment  of  the  discrete  geodesic  l\.  We  need  these  segments  of  their  boundaries  to 
lie  on  the  actual  geodesic  L\  in  Rk.  The  problem  is  easily  resolved,  at  the  expense  of 
admitting  some  quasiconformality  in  our  eventual  transformation,  by  the  application 
of  the  hyperbolic  projection  pk  described  in  Section  4.2,  as  illustrated  in  Figure  4.7. 


(a)  The  Pullback  of  U,  Pk(U) 


(b)  The  Projection  of  Pk(U),  p\{Pk(U)) 


Figure  4.7:  Projection  of  Geodesics,  Discrete  Geodesics,  and  Quasicollars 


This  pk  map  on  the  set  Pk{U )  is  1  +  e*-quasiconformal,  but  Proposition  4.2.1  guar¬ 
antees  that  plk  converges  uniformly  to  the  identity;  so  ek  — >  0  as  k  — >  oo.  Further, 
the  constant  ek  is  independent  of  the  choice  of  U.  This  independence  is  a  result  of 
two  facts.  First,  the  projection  associated  with  any  given  geodesic  L1  and  the  ac¬ 
companying  discrete  geodesic  £k  has  a  quasiconformality  constant  uniform  across  the 
entire  geodesic.  Second,  there  are  only  finitely  many  geodesics  for  which  we  must  use 
these  projections.  Taking  the  maximum  such  quasiconformal  constant  as  ek,  we  have 
independence  of  the  choice  for  U. 

Corresponding  to  the  set  pk(Pk(U))  C  Rk,  or  more  precisely  to  some  open  subset  of 
Rk  containing  p\(Pk{U))  C  Rk,  we  have  a  coordinate  chart  <pkv  which  takes  pk(Pk{U)) 
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to  a  subset  of  a  fundamental  region  for  Rk  in  the  Poincare  unit  disk,  B.  The  image 
of  plk(Pk{U))  pi  L\  under  this  map  is  thus  an  arc  of  a  hyperbolic  geodesic  in  ID).  There 
exists  a  disk  automorphism  Mk  :  ID)  — >-  ]D>  such  that  Mk  (cpkv(p{(Pk(U))))f]L[)  is 
a  small  arc  of  a  hyperbolic  geodesic  in  B  that  is  the  projection  of  the  geodesic  L1 
in  R.  We  may  choose  this  arc  of  the  projection  of  Ll  sufficiently  small  that  the  set 
Mk (ipkv (pl(Pk(U))))  is  contained  in  a  single  fundamental  region  for  R. 

Now,  a  second  difficulty  arising  from  the  projection  operation  must  be  addressed. 
If  we  were  to  apply  a  shearing  map  on  the  lift  of  the  geodesic  I)  which  separates  the 
“halves”  of  MjJ (cpkv  (p\ (Pk (U))))  using  the  weight  associated  with  L\  the  sides  will 
not  necessarily  match  along  the  geodesic,  as  illustrated  in  Figure  4.8. 


Figure  4.8:  Misalignment  under  Pullback,  Projection,  and  Shearing 

The  difference  along  the  geodesic  is  a  quasisymmetry  on  one  half  of  the  intersection 
with  L\  however,  and  extends  to  a  1 +efc-quasiconformal  map,  Qk,  in  the  correspond¬ 
ing  half  of  Mk(ipkv (plk(Pk(U))))  [24,  33].  Recall,  however,  that  the  quasisymmetry 
and  associated  quasiconformal  map  arise  as  a  result  of  the  need  for  a  hyperbolic  pro¬ 
jection  plk  to  correct  the  difference  between  the  geodesic  L\  and  the  discrete  geodesic 
£\.  Lemma  4.1.1  guarantees  that  £\  — >  L\  uniformly  as  k  —>  oo.  Thus,  Ek  — ►  0  as 
k  —>  oo.  Further,  since  the  quasiconformality  of  Mp  depends  only  on  the  quasisym¬ 
metry  associated  with  the  projection  plk,  ek  is  independent  of  the  choice  of  U. 

Corresponding  to  the  set  Q% (Mp (pkv{Pkipk(U)))))  C  B,  we  have  a  map  from  the 
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conformal  structure  on  R,  ip~l ,  taking  this  set  to  an  open  set  on  R  intersecting  the 
geodesic  Ll.  (Again,  the  coordinate  chart  really  comes  from  an  open  set  containing 
the  set  of  interest.)  Notice  that  as  a  result  of  applying  the  correction  Qk,  the  two 
halves  of  this  image  (on  either  side  of  the  lift  to  D  of  the  geodesic  Ll)  now  match 
the  action  of  the  shearing  caused  by  P_1.  We  now  apply  P-1,  to  obtain  a  set  V  C 
R.  Corresponding  to  this  open  set  we  have  a  coordinate  chart  from  the  conformal 
structure  on  R  defined  in  Section  2.6.  Since  V  €  UT,  as  described  in  Section  2.6,  we 
have  a  coordinate  chart  corresponding  to  V.  This  gives  the  following: 

■07  °  P~l  °  TV 1  °  Qk  °  Mk  °  TV  O  pi  oPko 

=  S7ocpyoPo  P'1  o  (p-1  oQVo  Mk  o  (pkv  o  p\  oPko 
=  S1oip1o  cp-1  oQuko  Mk  o  <pkv  op\oPko  yj"1. 

This  map  is  (1  +  ek)(  1  +  erfe)(l  +  £fe)-quasiconformal,  where  this  quasiconformality  is 
independent  of  the  choice  for  U.  Now,  consider  the  dilatation  at  a  point  u  e  U.  If  u 
is  not  contained  in  the  set  Ek(lk),  then  Lemma  3.2.4,  the  packing  Lemma,  guarantees 
that  the  dilatation  at  u  goes  to  1  since  each  of  ek,ik,ik  ->  0  as  k  -»  oo.  Thus,  this 
map  from  Rk  to  R  converges  to  a  conformal  map,  except  perhaps  on  the  image  of  the 
discrete  geodesic  i\  under  the  earthquake  Ek.  But  this  set  has  measure  zero,  and  the 
map  is  continuous.  Thus,  the  limit  map  is  conformal  on  all  of  U  [33]. 

Case  Illb:  Suppose  Im(a(L1))  >  0;  that  is,  the  earthquake  action  on  the  geodesic 
Ll  is  a  grafting  action  (and  potentially  a  shearing  action  as  well) .  The  situation  here 
is  somewhat  more  complicated  than  the  previous  case  since  the  map  Pk  takes  the 
portion  of  U  on  the  inserted  cylinder  to  a  cylinder  Cjk  and  the  remaining  points  of  U 
to  a  subset  of  Rk.  We  will  deal  with  the  maps  on  these  two  subsets  of  U  separately. 

First,  consider  that  subset  of  U  which  is  taken  to  Rk.  Call  this  subset  U\.  Pk(U\ )  ^ 
Ck  and  Pk(Ui)  ^  0.  We  wish  to  send  this  set  Pk{U\ )  to  a  subset  of  R,  but  there 
is  a  geometric  difficulty  to  be  resolved.  Pfc(Pi)  has  as  one  portion  of  its  boundary 
a  segment  of  the  discrete  geodesic  tk,  we  need  this  segment  of  the  boundary  to  lie 
on  the  actual  geodesic  L\  in  Rk.  As  before,  this  problem  is  easily  addressed,  at  the 
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expense  of  admitting  some  quasiconformality  into  our  eventual  transformation,  by 
the  application  of  the  hyperbolic  projection  p\  described  in  Section  4.2.  This  map 
on  the  set  Pk{U\)  is  1  +  e^-quasiconformal,  but  Proposition  4.2.1  guarantees  that  p\ 
converges  uniformly  to  the  identity;  so  ek  — >  0  as  k  — >  oo.  As  in  Case  Ilia,  ek  is 
independent  of  the  choice  of  U. 

Corresponding  to  the  set  pk(Pk(Ui))  C  Rk,  or  more  precisely  to  some  open  sub¬ 
set  of  Rk  containing  plk(Pk(Uij)  C  Rk,  we  have  a  coordinate  chart  <pkv  which  takes 
pk(Pk(Ui))  to  a  subset  of  a  fundamental  region  for  Rk  in  the  Poincare  unit  disk,  O.  The 
image  of  pk{Pk{U\))  f]Llk  under  this  map  is  thus  an  arc  of  a  hyperbolic  geodesic  in  D. 
There  exists  a  disk  automorphism  Mk  :  D  ->■  D  such  that  Mk(<pkv (pl(Pk(Ui))))  f]Pk) 
is  a  small  arc  of  the  hyperbolic  geodesic  in  D  that  results  from  the  lift  of  the  geodesic 
Ll  in  R.  We  may  choose  Mk  such  that  the  set  Mk{ipkv {p\{Pk(Ui))))  is  an  open  set 
intersecting  this  geodesic  in  a  fundamental  region  for  R. 

Again,  we  will  require  a  correction  to  compensate  for  the  effect  of  the  hyperbolic 
projection  p\.  Call  this  correction  Qk,  and  note  that  this  map  is  a  quasisymmetry 
on  the  intersection  with  the  lift  of  the  geodesic  L1  and  is  (1  +  ek )  -qu asiconforni al ,  as 
in  Case  Ilia.  Again,  the  quasiconformality  constant  ik  is  independent  of  the  choice 
for  U. 

Corresponding  to  the  set  Q% (Mf  (<Pkv{Pk(Pk(Ui)))))  C  D,  we  have  a  map  from  the 
conformal  structure  on  R,  <Pv\  taking  this  set  to  an  open  set  on  R  intersecting  the 
geodesic  L\  (Again,  the  coordinate  chart  really  comes  from  an  open  set  containing 
the  set  of  interest.)  We  now  apply  P~l,  to  obtain  a  set  V\  C  R  adjoining  the  edge  of 
an  inserted  cylinder  Cj. 

Now,  consider  that  subset  of  U  which  is  taken  into  Cjk.  Call  this  subset  U2- 
Pk(U2)  C  ch  and  Pk{U2)  intersects  the  boundary  of  the  cylinder  Cjk.  We  take  this  set 
into  the  cylinder  c,  via  the  1  +  e*-quasiconformal  map  fjk ,  where  ek  is  independent  of 
the  choice  of  U,  since  we  may  choose  lk  as  the  maximum  of  the  constants  associated 
with  the  finitely  many  cylinders  {cj}.  We  then  apply  the  map  P~l  to  obtain  a  set 
V2CR. 
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The  way  we  have  constructed  these  maps,  in  particular  the  application  of  the 
map  pi,  ensures  that  the  set  V  =  V\  \JV2  C  R  is  an  element  of  the  collection  UA 
as  described  in  Section  2.6.  Corresponding  to  this  set  V  we  have  a  coordinate  chart 
ips  from  the  conformal  structure  on  R  given  in  Section  2.6.  This  gives  the  following 

. — V  xN 

piecewise  map  Uk  :  Rk  ^  R- 

J ips  °  °  (PZ1  o  Qk  o  Mk  °  <Pkv  °p{  oPko  ifueUiCRk 

n  k{u)  =  < 

[i>6  o  P~l  o  fjh  oPko  (u)  if  u  e  U2  c  Rk. 

The  first  of  these  maps,  for  u  €  L\,  is  (1  +  ek)(l  +  ek){l  +  £fc)-quasicon for rn al .  The 
second  map,  for  u  G  U?,  is  (1  +  £jk){  1  +  e/t)-quasiconformal.  As  we  have  shown,  each 
of  these  constants  is  independent  of  the  choice  of  U.  Now,  consider  the  dilatation  at 
a  point  u  eU.  If  u  is  not  contained  in  the  set  Ek(£[),  then  Lemma  3.2.4,  the  packing 
Lemma,  guarantees  that  the  dilatation  at  u  goes  to  1  since  each  of  £k,  ek,  ik,  £jk  — »  0 

^  xN 

as  k  ->  oo.  Thus,  this  piecewise  map  Uk  from  Rk  to  R  converges  to  a  conformal 
map  II,  except  perhaps  on  the  image  of  the  discrete  geodesic  t\  under  the  earthquake 
Ek.  But  this  set  has  measure  zero,  and  the  map  is  continuous.  Thus,  the  map  II  is 
conformal  on  all  of  U  [33]. 

Thus,  as  k  -»  oo  the  maps  from  Rk  to  R  are  becoming  conformal,  and,  by  con- 

x»v  XV. 

struction,  the  markings  on  surfaces  R  and  Rk  are  consistent.  Therefore,  the  surfaces 

XV  XV 

Rk  converge  to  a  surface  conformally  equivalent  to  R  with  equivalent  markings;  but 
this  is  equivalent  to  saying  Rk  — >■  R  (in  the  Teichmiiller  metric)  as  k  — >  oo,  and  the 
result  is  shown.  □ 
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CHAPTER  V 

EARTHQUAKES  ON  EUCLIDEAN  SURFACES 

In  the  work  of  Thurston  [51],  McMullen  [35],  and  others  [16,  21,  23,  24,  25,  40,  54], 
earthquake  maps  are  confined  to  hyperbolic  Riemann  surfaces.  To  this  point,  we  have 
similarly  limited  the  scope  of  our  attention,  examining  only  the  approximation  of  these 
complex  earthquakes,  shearing  and  grafting  maps  on  hyperbolic  Riemann  surfaces. 
We  now  turn  our  attention  to  the  question  of  whether  or  not  we  can  meaningfully 
define  earthquake  maps  on  the  Teichmiiller  space  of  Euclidean  Riemann  surfaces.  In 
this  chapter  we  address  four  basic  issues  related  to  earthquakes  on  tori.  Initially,  we 
have  a  question  of  definition;  we  must  define  those  earthquakes.  Given  this  definition, 
we  have  three  problems  we  must  address.  First,  we  must  determine  whether  or  not 
the  earthquakes  thus  defined  are  transformations  of  points  in  the  Teichmiiller  space 
of  tori.  That  is,  is  the  image  of  a  torus  under  a  Euclidean  earthquake  another  torus? 
Second,  we  must  determine  whether  or  not  the  Euclidean  earthquakes  thus  defined 
exhibit  properties  similar  to  those  of  the  hyperbolic  earthquakes.  Finally,  can  we 
approximate  the  action  of  these  Euclidean  earthquakes  using  circle  packing? 

5.1  Defining  Earthquakes  on  the  Torus 

From  the  perspective  of  shearing  and  grafting  maps  on  the  Riemann  surface,  we 
define  (Euclidean)  earthquakes  on  tori  in  a  manner  analogous  to  the  manner  in  which 
those  maps  were  defined  on  hyperbolic  surfaces.  That  is,  we  define  a  shearing  by 
opening  the  torus  along  a  geodesic,  twisting  by  a  prescribed  amount,  and  gluing  the 
surface  back  together  along  the  seam  of  the  geodesic.  Similarly,  we  define  a  grafting 
by  opening  the  torus  along  a  geodesic,  inserting  a  cylinder  of  some  prescribed  height, 
and  gluing  the  surface  together. 

To  define  such  an  action,  we  first  need  to  discuss  finite  measured  geodesic  lam¬ 
inations  of  tori.  Recall  that  in  constructing  earthquakes  on  hyperbolic  surfaces  we 
used  the  Poincare  unit  disk  and  the  upper  half  plane  (with  appropriate  metrics)  to 
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describe  the  earthquake  action.  We  have  a  similar  construction  here,  where  we  use 
the  Euclidean  (rather  than  hyperbolic)  upper  half-plane  H  as  a  lifting  cover  for  the 
Euclidean  surface,  the  torus.  We  restrict  ourselves  to  the  upper  half-plane  so  that 
each  point  in  the  Teichmuller  space  of  tori  is  uniquely  defined  (after  normalization) 
by  a  single  point  uel  [29,  30].  We  now  define  a  finite  Euclidean  geodesic  lamination 
on  a  torus. 

Definition  5.1.1.  A  finite  geodesic  lamination  C  of  the  Euclidean  torus  T  is  a 
collection  of  finitely  many  disjoint  geodesics  on  T  which  lifts  to  a  collection  of  disjoint 
( and  thus  parallel)  Euclidean  lines  in  the  upper  half-plane  H. 

Note  that  the  structure  of  the  Euclidean  plane  dictates  that  in  a  finite  geodesic 
lamination  of  the  plane  consisting  of  more  than  a  single  geodesic,  those  geodesics  must 
all  be  parallel  (in  a  Euclidean  sense).  This  limits  substantially  the  forms  that  such  a 
lamination  may  take  in  the  Euclidean  case;  we  have  less  freedom  in  constructing  these 
laminations  than  in  the  hyperbolic  case  since  we  have  only  one  “point  at  infinity.” 
Just  as  in  the  hyperbolic  case,  we  place  a  measure  on  the  lamination  C. 

Definition  5.1.2.  A  finite  geodesic  lamination  C  of  the  Euclidean  plane  C  is  a 
collection  of  finitely  many  disjoint  Euclidean  lines  (geodesics).  In  this  way,  C  divides 
H  into  disjoint  strips. 

Recall  that  the  construction  of  complex  earthquakes  on  hyperbolic  surfaces  re¬ 
quired  the  use  of  hyperbolic  isometries.  In  the  Euclidean  complex  plane,  the  isome¬ 
tries  we  need  are  simply  translations  z  +  a,  a  £  C.  A  shearing  map  arises  from  a 
translation  with  axis  parallel  to  the  geodesic  along  which  we  wish  to  shear,  and  a 
grafting  map  arises  from  a  translation  with  axis  perpendicular  to  the  geodesic  along 
which  we  intend  to  graft.  In  grafting,  we  insert  finite  cylinders,  and  in  the  plane  this 
is  represented  by  the  insertion  of  a  finite  strip  into  the  gap  created  by  the  piecewise 
application  of  the  grafting  translation  (with  axis  perpendicular  to  the  geodesic).  Note 
that  we  may  select  this  finite  strip  (essentially  a  Euclidean  parallelogram)  so  that  its 
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generators  (and  thus  the  marking)  are  consistent  with  the  generators  (and  thus  the 
marking)  that  defines  the  Teichmiiller  space  under  examination. 

We  model  the  actions  of  shearing  and  grafting  on  a  torus  on  the  fundamental 
regions  for  the  tori  in  the  upper  half-plane,  just  as  we  used  the  upper  half-plane  (and 
the  Poincare  unit  disk)  to  model  hyperbolic  earthquakes.  These  four  actions,  left 
and  right  shearing  and  positive  and  negative  grafting,  are  illustrated  in  Figure  5.1. 
We  then  define  complex  Euclidean  earthquakes,  earthquakes  with  complex  weights, 
as  the  composition  of  Euclidean  shearing  and  then  Euclidean  grafting  operations. 
Recall  that  in  the  hyperbolic  case  we  claimed  that  the  term  complex  earthquake 
was  illustrative  in  that  the  earthquake  action  was  described  as  multiplication  by  a 
complex  number  in  the  local  model  on  H.  A  similar  idea  holds  in  the  case  of  Euclidean 
earthquakes;  rather  than  multiplication  by  a  complex  number,  we  now  transform 
defining  points  in  the  fundamental  region  through  the  addition  of  a  complex  number. 


Im(z)  Im(z)  Irn(z) 


(a)  Fundamental  Region  and  (b)  Right  Shearing  (c)  Left  Shearing 

a  Geodesic 

Im(z)  Ira(z) 


(d)  Positive  Grafting  (e)  Negative  Grafting 

Figure  5.1:  Earthquake  Actions  on  the  Fundamental  Region  of  a  Torus 

Note  that  in  constructing  these  Euclidean  earthquake  maps,  we  place  some  re¬ 
strictions  on  the  weights.  First,  consider  the  imaginary  part  (grafting  component)  of 
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the  complex  weights.  In  [35],  positive  grafts  of  any  size  were  permitted,  but  negative 
grafts  (amounting  to  removal  of  portions  of  the  surface)  required  some  restrictions 
for  obvious  geometric  reasons.  We  must  restrict  negative  grafting  weights  similarly 
in  the  Euclidean  case.  There  are  two  reasons  for  this  restriction.  First,  we  do  not 
want  to  remove  so  much  in  a  negative  graft  along  a  particular  geodesic  that  we  re¬ 
move  another  geodesic.  Second,  we  must  not  remove  all  of  the  torus  remaining  in  one 
direction. 


5.2  A  Conformal  Structure  on  the  Images  of  Tori  under 
Complex  Earthquake  Maps 

We  can  now  turn  to  the  first  of  our  outstanding  problems,  Problem  5.2.1. 

Problem  5.2.1.  Do  earthquakes  on  Euclidean  surfaces,  and  in  particular  compact 
tori,  result  in  Riemann  surfaces?  That  is,  are  shearing  and  grafting  maps  defined  on 
tori  transformations  in  the  Teichmiiller  space  of  tori? 

Suppose  that  we  apply  a  finite  complex  Euclidean  earthquake  E  induced  by  a  finite 
measured  geodesic  lamination  (£,  a)  to  the  torus  T  (a  compact  Euclidean  Riemann 
surface),  where  the  lamination  consists  of  n  geodesics  with  complex  weights  a(Li)  = 
Pq  +  iXq,  1  <  q  <n,  where  pq  G  R.  and  A,Gl  such  that  Xq,  if  negative,  is  sufficiently 
small  so  that  the  execution  of  the  negative  grafting  operation  does  not  remove  another 
geodesic  in  the  lamination  and  does  not  reduce  the  torus  to  a  degenerate  surface.  This 
earthquake  thus  involves  both  shearing  maps  and  grafting  maps  composed  to  form 
the  complex  earthquake. 

Since  T  is  a  fixed  Riemann  surface,  it  has  a  conformal  structure  associated  with  it; 
say  {ipv}ver  1S  the  atlas  for  the  conformal  structure  on  T.  Further,  since  each  grafting 
map  with  positive  imaginary  weight  involves  the  insertion  of  a  cylinder  Cj,  we  have  a 
conformal  structure  associated  with  each  such  cylinder  so  that  the  coordinate  charts 
map  Cj  to  subsets  of  a  Euclidean  strip.  For  each  of  Cj ,  0  <  j  <m,  where  m  gives  the 
number  of  geodesics  in  the  lamination  £  which  have  positive  imaginary  components 
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for  their  weights,  let  {<f> i}keK  be  t^ie  a^as  associated  with  each  such  cylinder. 

Define  a  pullback  map  P  on  T  so  that  the  image  of  a  point  r  €  T  under  P  is  the 
unique  point  (on  T  or  on  one  of  the  cylinders  c,)  whose  image  under  the  complex 
earthquake  E  is  r. 

With  this  information  we  now  define  a  conformal  structure  on  the  Riemann  surface 
T,  the  image  of  T  under  the  complex  earthquake  E.  To  define  this  conformal  structure 
we  need  only  describe  a  set  of  coordinate  charts  mapping  open  regions  in  the  surface  T 
into  C  so  that  the  transition  maps  associated  with  the  coordinate  charts  are  analytic. 
It  is  sufficient  to  construct  for  every  point  r  e  T  a  map  from  an  open  neighborhood 
U  C  T  into  C  so  that  these  maps  satisfy  the  analyticity  condition.  We  consider  open 
sets  of  the  following  five  classes. 

1.  U a  is  the  collection  of  all  open  sets  Ua  C  T  such  that  TJa  does  not  intersect  the 
image  under  the  earthquake  E  of  any  geodesic  Li,  1  <  i  <  n,  and  is  disjoint 
from  every  inserted  cylinder  Cj,  1  <  j  <  m. 

2.  U n  is  the  collection  of  all  open  sets  Up  C  T  such  that  Up  C  E(cj)  for  some 

j,  1  <j<m. 

3.  Ur  is  the  collection  of  all  open  sets  t/7  C  T  such  that  U1  intersects  E(Li)  for 
exactly  one  value  i,l  <  i  <  n,  where  a  (Li)  e  R,  (i.e.,  the  only  map  on  that 
geodesic  is  a  shearing  operation). 

/■S 

4.  Ua  is  the  collection  of  all  open  sets  Us  CT  such  that  Us  intersects  both  P~1(R) 
and  E(cj)  for  some  j,  1  <  j  <  m. 

5.  Us  is  the  collection  of  all  open  sets  U$  C  T  such  that  U^  intersects  the  image 
of  the  boundary  of  some  cylinder  removed  by  the  action  of  a  negative  grafting 
action. 

An  illustration  of  representative  sets  from  these  classes  is  given  in  Figure  5.2. 
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Figure  5.2:  Open  Sets  to  Define  a  Conformal  Structure  on  the  Torus  T 


Case  I:  Let  Ua  G  UA,  and  note  that  P(Ua)  is  an  open  set  in  T.  From  the  conformal 
structure  on  T  and  corresponding  to  this  open  set  P(Ua)  we  have  a  coordinate  chart 
tpa  such  that  (pa(P(Ua ))  C  C.  We  thus  define  a  coordinate  chart  ipa  :  T  -»  C  on  the 
open  set  Ua  by  rpa  —  <Pa°  P-  In  this  way  we  define  a  family  of  coordinate  charts  ^A 
on  f  by  =  {^a}aeA. 


Figure  5.3:  Construction  of  Coordinate  Charts  for  Open  Sets  Up 


Case  II:  Let  Up  <G  UB-  In  this  case,  Up  C  Int(£'(cj)),  for  some  1  <  j  <  m,  and 
the  pullback  map  P  takes  this  set  to  an  open  set  on  the  finite  cylinder  Cj  as  shown  in 
Figure  5.3a.  Now,  from  the  conformal  structure  on  Cj  we  have  a  coordinate  chart  (f>p 
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corresponding  to  P{Up)  such  that  <frp(P(Up))  is  a  subset  of  an  infinite  Euclidean  strip 
of  uniform  width  determined  by  the  complex  part  of  the  weight  on  the  appropriate 
geodesic  in  C.  We  then  have  a  grafting  map  that  opens  the  plane  along  the 
corresponding  geodesic  in  C  and  glues  in  the  strip  as  shown  in  Figure  5.3b.  This  gives 
us  a  map  ipp  :  T  — >  C  on  the  open  set  Up  defined  by  ipp  —  o  <pJ^  o  P.  In  this  way  we 
define  a  family  of  coordinate  charts  ^B  on  T  by  $B  =  {^p\p€B- 


(a)  Image  of  Uy  C  T  under  P  (b)  Image  of  I/7  C  T  under  ipy  o  P 


(c)  Image  of  C/7  C  T  under  S1oip1oP 

Figure  5.4:  Construction  of  Coordinate  Charts  for  Open  Sets  t/7 

Case  III:  Let  f/7  £  C/r,  and  note  that  P(f/7)  is  a  set  as  shown  in  Figure  5.4a.  From 
the  conformal  structure  on  T  and  corresponding  to  an  open  set  containing  P(P7) 
we  have  a  map  (p1  such  that  <y97(P(£/7))  C  C.  Now,  we  have  disjoint  sets  that  differ 
by  the  application  of  a  Euclidean  shearing  map  (a  piecewise  translation),  Sy,  along 
the  axis  determined  by  the  geodesic  shown  in  Figure  5.4b.  We  have  thus  defined  a 
coordinate  chart  t/>7  :  T  — »  C  on  the  open  set  U1  by  -07  =  SjOip^o  P.  In  this  way,  we 
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(c)  Image  of  P(Ug  f)(E(cj)))  under  <f)jg  (d)  Image  of  P(Ug  n(S(Cj))c)  under  ipg 


Figure  5.5:  Construction  of  Coordinate  Charts  for  Open  Sets  Ug 


Case  IV:  Let  Ug  £  C/a,  and  note  that  P(Ug)  is  divided  between  two  distinct 
Riemann  surfaces  Cj,  as  shown  in  Figure  5.5a,  and  T,  as  shown  in  Figure  5.5b.  To 
define  coordinate  charts  on  set  of  the  form  Ug  we  define  a  pair  of  maps  from  these 
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surfaces  to  the  upper  half-plane  H.  Beginning  with  the  cylinder,  we  have  a  map 
such  that  (fjg(P(E(cj )  fl  Us))  is  a  Euclidean  strip  of  uniform  width  determined  by  the 
complex  part  of  the  weight  on  the  appropriate  geodesic  in  C  as  shown  in  Figure  5.5c. 
At  the  same  time,  from  the  conformal  structure  on  T  we  have  a  map  ips  such  that 
ips(P((E(cj))c  fl  Us))  C  H  as  shown  in  Figure  5.5d.  Applying  the  grafting  map  gJ5  we 
put  the  pieces  back  together  in  H.  Now,  we  may  need  to  shear  along  the  hyperbolic 
geodesic  whose  image  to  this  point  is  the  imaginary  axis  in  Figure  5.5d;  this  will 
depend  upon  whether  or  not  the  geodesic  associated  with  this  grafting  action  had 
a  weight  with  non-zero  real  part.  Note  that  the  shearing  will  take  place  along  that 
portion  of  the  set  not  intersecting  E(cj),  since  we  define  complex  earthquakes  as  the 
composition  of  shearing  and  grafting  maps,  in  that  order. 


f  si  o  o  p(t)  if  re  E(cj)  nuscf 
j^o^'o^o  P(t)  if  r  €  (E(cj))c  n  Us  C  T 


This  function  now  defines  a  map  ips  :T  C  on  the  open  set  Us,  as  shown  in  Figure 
5.5e.  Taking  all  such  functions  over  the  set  of  possible  sets  Us  €  C/a,  we  define  a 
family  of  coordinate  charts  \?a  =  A- 

Case  V:  Let  U ^  €  Us,  and  note  that  P(U $)  is  a  set  as  shown  in  Figure  5.6a.  From 
the  conformal  structure  on  T  and  corresponding  to  an  open  set  containing  P(UI j)  we 
have  a  map  cp^  such  that  < p^(P(U f))  C  C.  Now,  we  have  disjoint  sets  that  differ  by  the 
application  of  a  negative  Euclidean  grafting  map  (a  piecewise  translation),  g%,  along 
the  axis  perpendicular  to  the  geodesic  shown  in  Figure  5.6b  and  a  negative  Euclidean 
shearing  map  (piecewise  translation)  along  the  axis  shown  in  Figure  5.6b.  We  have 
thus  defined  a  coordinate  chart  :  T  — >  C  on  the  open  set  by  ip£  =  S^og^ocp^oP. 
In  this  way,  we  define  a  family  of  coordinate  charts 

Define  a  collection  of  open  sets,  U  =  Ua  U  Ub  U  Ut  U  Ua  U  C/e,  and  a  collection 
of  maps  from  T  to  C,  W  =  ^ a  IJ  U  U  U  Given  this  collection  of  open 
sets  and  corresponding  coordinate  charts  on  the  image  T,  Proposition  5.2.1  provides 
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(c)  Image  ofU^CT  under  S^og^o^oP 
Figure  5.6:  Construction  of  Coordinate  Charts  for  Open  Sets  U $ 

an  answer  to  Problem  5.2.1;  we  prove  that  T  is,  in  fact  a  Riemann  surface.  This  leads 
us  immediately  to  Corollary  5.2.1,  that  Euclidean  earthquakes  are  transformations 
on  the  Teichmiiller  space  of  tori. 

Proposition  5.2.1.  The  transition  maps  in  the  structure  ( U ,  T)  are  analytic. 

Proof.  Let  U  be  the  intersection  of  two  open  sets  in  U.  We  need  only  consider  those 
intersections  for  which  U  £  U.  This  gives  five  cases  for  the  overlap  region  U  in  which 
we  must  verify  the  analyticity  of  the  transition  maps. 

Case  I:  Suppose  that  U  G  Ua ■  There  are  nine  ways  in  which  U  may  occur  as  the 
result  of  intersections  of  sets  in  U. 

Case  la:  U  =  Uai  f]Ua2,  where  Uai , Ua2  G  Ua-  Associated  with  Uai  and  Ua2  we 
have  coordinate  charts  ipai,ifa2  G  Say  ipai  =  (pai  o  P  and  ipa2  =  <pa2  o  P.  Now, 
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consider  the  transition  maps  ipa2  o  ipax  and  ipai  o  ipax. 

^a2  O  lj!~l  =  (ipa2  oP)o  (ipai  O  P)-1 

=  (pa2oPo  P~l  O 

=  V«2  ° 

Thus,  ipa20'll;ai  :  — >  C  is  analytic,  since  ipai,  i  =  1,2,  were  taken  from  the  conformal 

structure  on  T,  making  ipa2  o  (p~x  analytic.  Similarly,  (pai  o  ip~x  is  analytic. 

Case  lb:  U  =  Uaf]  C/7,  where  Ua  G  U a  and  C/7  G  Ur-  Associated  with  Ua  and  IJ1 
we  have  coordinate  charts  ipa  G  a  and  ip-,  G  Tr,  respectively.  Say  ipa  =  pao  P  and 
ip-,  =  S-,  o  ip-,  o  P.  Now,  consider  the  transition  maps  ipa  o  ip~x  and  ip-,  o  ip~l. 

ipa  °  #7 1  =  (Va  °P)°  (S-y  otp^o  P )_1 

=  <pa  o  P  O  p-1  o  (p-1  oS-1 
=  Vcl  O  1  O  S -1 

Note  that  the  transition  map  ipa  o  ip-1  is  defined  only  on  the  set  U.  Since  the  set 
U  does  not  itself  intersect  the  geodesic  associated  with  the  shearing  map  S7,  on  U 
the  shearing  map  acts  as  either  the  identity  map  or  as  a  translation.  In  either  case, 
the  action  of  S7  is  (locally)  analytic.  Next,  we  note  that  since  <pa  and  p-f  are  taken 
from  the  conformal  structure  on  T,  <pa  o  i p~l  is  analytic,  where  defined.  Thus,  the 
transition  map  ipaOip*1  is  the  composition  of  analytic  maps,  and  is  therefore  analytic. 
Similarly,  ip 7  o  ip~l  is  also  analytic. 

Case  Ic:  U  =  Uaf)U$,  where  Ua  G  Ua  and  U$  G  Ua-  Associated  with  each 
of  Ua  and  Ug  we  have  coordinate  charts  ipa  G  ^  a  and  ipg  G  respectively.  Say 
ipa  =  <pa  °  P  and  ips  of  the  form  given  in  (5.2.1).  We  note  that  since  U  does  not 
intersect  the  geodesic  associated  with  the  grafting  of  the  cylinder,  for  every  r  G  U 
we  have  r  G  ( E(cj))c  n  Ug  C  T,  where  j  denotes  the  cylinder  associated  with  the 
grafting.  Thus,  on  U  we  have  ipg  =  Sg  o  o  ips  o  P.  Now,  consider  the  transition 
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maps  ipaOipg1  and  ips  °  ipa  1  • 

Ipa  O  Ipg1  =  ((pa  oP)o  (SS  ogjocpsO  P)”1 

=  (pa  O  P  o  P-1  o  (pj1  O  (gfl  1  o  Sg1 
=  Voc  o  1  °  (s£)  1  o  ^_1 

Note  that  we  define  the  transition  map  ipa  o  ipg1  on  only  on  the  overlap  region  U. 
Since  the  set  U  does  not  itself  intersect  the  inserted  cylinder  (or  the  boundary  along 
which  we  weld  the  cylinder  to  the  Riemann  surface  T  via  the  grafting  map  (<^)  ' ) ,  on 
U  the  inverse  of  the  grafting  map  acts  as  either  the  identity  map  or  as  a  translation. 
In  either  case,  the  action  of  (gj)  1  is  (locally)  analytic.  For  the  same  reason,  on  U 
we  have  that  Sg 1  acts  as  either  the  identity  map  or  a  translation,  and  is  therefore 
(locally)  analytic.  Next,  we  note  that  since  the  maps  ipa  and  g>g  were  borrowed  from 
the  conformal  structure  on  the  original  Riemann  surface  T,  g>ao<pppl  is  analytic,  where 
defined.  Thus,  the  transition  map  ipa  o  ipg1  is  the  composition  of  analytic  maps  on 
U,  and  is  therefore  analytic.  Similarly,  tpg  o  ij)al  is  also  analytic. 

Case  Id:  U  =  U1  f)  Us,  where  t/7  €  Ur  and  Us  €  U&.  Associated  with  each 
of  Uf  and  Us  we  have  coordinate  charts  rp1  G  Tr  and  ips  G  Ta,  respectively.  Say 
ipry  =  57  o  <p7  o  P  and  ips  of  the  form  given  in  (5.2.1).  We  note  that  since  U  intersects 
neither  the  geodesic  associated  with  the  grafting  of  the  cylinder  nor  the  geodesic 
associated  with  the  simple  shearing,  for  every  r  eU  we  have  r  G  (E(cj))c  HUs  CT, 
where  j  denotes  the  cylinder  associated  with  the  grafting.  Thus,  on  U  we  have 
ips  =  Ss  o  gj  o  (ps  o  P.  Now,  consider  the  transition  maps  ipy  o  ipj 1  and  ips  o  ip~l. 

ip'r  °  iPs1  =  (St  °  V-y  o  P)  o  (Ss  o  gi  o  (ps  o  P)-1 

=  57o^7oPoF-1o<^-1  o  (^)_1  oSg1 
=  S1oip1o  tp~l  o(g3g)  1  o  Sg 1 

Note  that  we  define  the  transition  map  ip7  o  ipg1  only  on  the  overlap  region  U.  Since 
the  set  U  does  not  itself  intersect  the  inserted  cylinder  (or  the  boundary  along  which 
we  weld  the  cylinder  to  the  Riemann  surface  T  via  the  grafting  map  (g^)  ^  on  U 
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the  inverse  of  the  grafting  map  acts  as  either  the  identity  map  or  as  a  translation.  In 
either  case,  the  action  of  (g3s)  1  is  (locally)  analytic.  For  the  same  reason,  on  U  we 
have  that  Sg1  and  S7  each  act  as  either  identity  maps  or  translations;  they  are  thus 
(locally)  analytic.  Next,  we  note  that  since  the  maps  cp7  and  tps  were  borrowed  from 
the  conformal  structure  on  the  original  torus  T,  ipa  o  ipj1  is  analytic,  where  defined. 
Thus,  the  transition  map  ip7  o  ip#1  is  the  composition  of  analytic  maps  on  U ,  and  is 
therefore  analytic.  Similarly,  tps  °  ipZ1  is  also  analytic. 

Case  Ie:  U  =  Uaf]  U^,  where  Ua  G  Ua  and  C/f  G  C/3.  Associated  with  Ua  and  C/f 
we  have  coordinate  charts  tpa  a  and  ipt  G  T3,  respectively.  Say  tpa  =  <pa  o  P  and 
tpt  =  o  o  (p£  o  P.  Now,  consider  the  transition  maps  ipa  o  tpp1  and  ip£  o  tp~l. 

^  0  '•Pa1  =  (St  °9tonoP)°  (v*  °  py1 
=  S(O^O^oPo  P_1  o  <p~l 

=  St  O  o  ipt  o  (f-1 

Note  that  we  define  the  transition  map  ipt  o tpal  only  on  the  set  U.  Since  U  does  not 
itself  intersect  either  the  geodesic  associated  with  U^,  on  U  the  shearing  and  grafting 
maps  act  as  either  identity  maps  or  as  translations.  In  either  case,  the  action  of  g^ 
and  St  are  (locally)  analytic.  Next,  we  note  that  since  the  maps  <pa  and  (fit  were 
borrowed  from  the  conformal  structure  on  the  original  torus  T,  ipt  o  p~l  is  analytic, 
where  defined.  Thus,  the  transition  map  ipt  o  ip”1  is  the  composition  of  analytic  maps 
on  U,  and  is  therefore  itself  analytic.  Similarly,  ipQ  o  tpj1  is  also  analytic. 

Case  If:  U  =  C/7  P|  C/^,  where  C/7  G  Ur  and  C/g  G  C/3.  Associated  with  C/7  and  C/^ 
we  have  coordinate  charts  tp7  G  \&r  and  ipt  G  $3,  respectively.  Say  tp7  =  S7  o  ip7  o  P 
and  'ipt~9t,0(P(,°  p-  Now,  consider  the  transition  maps  ip7  o  xpj1  and  ipt  o  tp~l . 

ipt  o  ^7 1  =  (St  o  gt  o  ipt  o  P)  o  (S7  o  (p7  o  P)_1 
=  S^o^o^oPoP-'o  ip"1  o  S~l 
=  St°gt°<Pt°  ¥>7 1  0  Sj1 

Note  that  we  define  the  transition  map  tpt  o  tp~l  only  on  the  set  U.  Since  U  does  not 
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itself  intersect  either  the  geodesic  associated  with  the  set  t/7  or  the  geodesic  associated 
with  the  set  U^,  on  U  the  maps  S’"1,  g%,  and  S%  each  act  as  either  identity  maps  or 
translations.  In  either  case,  the  action  of  each  on  U  is  (locally)  analytic.  Finally,  we 
note  that  since  the  maps  py  and  p%  were  borrowed  from  the  conformal  structure  on 
the  original  torus  T,  p^  o  p~l  is  analytic,  where  defined.  Thus,  the  transition  map 
ip%  o  ip-1  is  the  composition  of  analytic  maps  on  U,  and  is  therefore  itself  analytic. 
Similarly,  ipy  o  ippp1  is  also  analytic. 

Case  Ig:  U  =  Us  f]U^,  where  Us  G  C/a  and  U^  G  U~,.  Associated  with  Us  and  U% 
we  have  coordinate  charts  ips  €  Ta  and  ip%  €  respectively.  Say  ip£  =  °  ° 

o  P  and  ips  is  of  the  form  given  in  (5.2.1).  Since  U  intersects  neither  the  geodesic 
associated  with  Us  nor  the  geodesic  associated  with  C/^,  for  every  r  G  U  we  have 
r  G  ( E(cj))c  C\Us  C  T,  where  j  denotes  the  cylinder  associated  with  the  grafting. 
Thus,  on  U  we  have  ips  =  Ss  °  °  ps  °  P.  Now,  consider  the  transition  maps  ip%  o  ipj1 
and  ips  o  ip^1. 

ips  o  ip ^  -  (Ss  o  g*  o  ps  o  P)  o  (Se  og^op^o  P)"1 
=  Ss  o  g3s  o  ps  0  P  °  P~l  o  o  g~l  o  S^1 
=  Ssogjopso  p^1  o  g~l  o  S'"1 

Again,  the  transition  map  ips  o  ip-1  is  defined  only  on  the  overlap  region  U.  Since 
this  set  intersects  neither  the  geodesic  associated  with  Us  nor  the  geodesic  associated 
with  t/je,  on  C7  the  maps  Ss,  S^1,  gJ6,  and  g^1  each  reduces  to  either  the  identity  map 
or  a  translation,  and  thus  each  is  (locally)  analytic.  Since  the  maps  ps  and  p%  were 
borrowed  from  the  conformal  structure  on  the  original  torus  T ,  ps  °  P^1  is  analytic, 
where  defined.  Thus,  the  transition  map  ips  o  ip^1  is  the  composition  of  analytic  maps 
on  U,  and  is  therefore  itself  analytic.  Similarly,  ip£  o  ip^1  is  also  analytic. 

Case  Ih:  U  —  U^f]U^2,  where  U^ ,  b\2  G  Us.  Associated  with  U^  and  U&  we 
have  coordinate  charts  ip^  G  dA,  z  =  1,2.  Say  ip£.  =  S^  o  g^.  o^.oP,i  =  1, 2.  Now, 
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consider  the  transition  maps  ip^  o  ip^  and  ip^2  o  ip^. 

Ai  °  tib  =  (%  °  9n  1  °  <Pb  °P)°  (<%  O  o  pf2  o  P)-1 
=  %  O  06  o^oFo  p-1  o  p"1  o  5”1  O  p-1 

=  %  °  °  °  °  °  'S'f/ 

The  transition  map  ip^  o  ip^1  is  defined  only  on  the  overlap  region  U.  Since  this  set 
intersects  neither  the  geodesic  associated  with  nor  the  geodesic  associated  with 
U& ,  on  U  the  maps  ,  S^1 ,  ,  and  g^1  each  reduces  to  either  the  identity  map  or 
a  translation,  and  thus  each  is  (locally)  analytic.  Since  the  maps  and  y?£2  were 
borrowed  from  the  conformal  structure  on  the  original  torus  T,  ip^  o  is  analytic, 
where  defined.  Thus,  the  transition  map  ip^  is  the  composition  of  analytic  maps 
on  U,  and  is  therefore  itself  analytic.  Similarly,  ip%2  o  ip^1  is  also  analytic. 

Case  II:  Suppose  that  U  eUb-  There  are  two  ways  in  which  U  may  occur  as  the 
result  of  intersections  of  sets  in  U. 

Case  Ha:  U  =  Up1f'\Up2,  where  UpliUp2  £  Ub-  Associated  with  Up1  and  Up2 
we  have  coordinate  charts  ‘tpp1,'fpp2  £  ^b-  Say  ip ^  =  g^  o^oP  and  ipp2  —  g^2  o 
<p3p2  o  P.  Note  that  since  both  Up1  and  Up2  must  necessarily  lie  in  the  same  finite 
inserted  cylinder,  the  grafting  maps  g ^  and  g^2 ,  are  the  same  map.  Now,  consider 
the  transition  maps  ipp2  o  ip and  ippt  o  ip^. 

i’fh  °  =  {gp2  °  <t>p2  oP)°  (4i  °  °  p)_1 

=  Sj2°  ^2  °P0P'l0K,)’1°W’1 

=  i°  (<4°  WJ"1)  °i4i)1 

Now,  notice  that  the  maps  g ^  and  gj^2  have  no  effect  on  the  analyticity  of  the  maps 
involved,  since  both  Up1  and  Ufa  are  on  the  interior  of  the  cylinder  E(cj).  The  action 
affecting  analyticity  takes  place  in  the  conformal  structure  on  the  cylinder  Cj ,  and 
(pp2  o  (<^)  1  is  analytic  since  these  maps  were  taken  from  this  conformal  structure. 
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Thus,  as  a  composition  of  analytic  maps,  ipp2  o  ipp^  is  analytic.  Similarly,  ipp1  o  ip^  is 
analytic. 

Case  lib:  U  =  Up  f)  Ug,  where  Up  G  Ub  and  Ug  G  U&.  Associated  with  each 
of  Up  and  Ug  we  have  coordinate  charts  'ipp  E  ^e  and  ipg  E  Ta,  respectively.  Say 
ipp  =  gpoftpoP  and  ipg  are  of  the  form  given  in  (5.2.1).  Note  that  since  U  is  does  not 
intersect  any  geodesic  associated  with  the  grafting  of  the  cylinder,  for  every  r  G  U  we 
have  t  G  E(Cj)  C\Ug  C  T,  where  j  denotes  the  cylinder  associated  with  the  grafting. 
Thus,  on  U  we  have  ipg  =  gj  o  (p]6  o  P.  Now,  consider  the  transition  maps  ipp  o  ipp1 
and  ipg  o  ip^1. 

ippoipfi1  =  (jJo^oPjogofjOP)"1 

This  case,  thus  reduces  to  Case  Ha,  and  the  transition  maps  tpp  o  ipp1  and  ipg  o  ipp1 
are  analytic. 

Case  III:  Suppose  that  U  G  Ur.  Then,  U  =  C/7l  f)  Ul2 ,  where  Ull ,  U12  G  C/r- 
Associated  with  each  of  C/7l  and  Ul2  we  have  coordinate  charts  Vbi  >  ^72  G  'hp.  Say 
'07l  =  S~f  o  cp7l  o  P  and  ipT2  =  57  o  (p12  o  P.  Note  that  the  earthquake  (shearing)  map 
in  each  of  these  coordinate  charts  is  the  same,  since  this  earthquake  is  defined  by 
intersection  with  a  particular  weighted  geodesic  rather  than  a  particular  set  U.  Now, 
consider  the  transition  maps  ipl2  o  ippp  and  ip^  0  V^1- 

^72  0  ^  =  (Sy  0  <Pi2  °P)°  (Sy  0  <Pyi  0  P)~l 
=  Sj  o  ipl2  o  p  o  p~l  o  ^  o  Spl 
=  A7  O  Vl2  o  p"1  o  S'”1 

The  action  of  the  pre-composition  and  post-composition  by  grafting  and  shearing 
maps  has  no  effect  on  the  analyticity;  these  maps  each  essentially  undo  the  local 
behavior  of  the  other.  Thus,  the  transition  map  ip12  o  ipp1  is  analytic.  Similarly, 
ip7 1  o  ippp  is  also  analytic. 

Case  IV:  Suppose  that  U  G  C/a.  U  =  Ug1f]Ug2,  where  Ug1,Ug2  G  C/a-  Associated 
with  each  of  Ug1  and  Ug2  we  have  coordinate  charts  ipgt ,  ipg2  G  ^a-  Say  tpg1  and  ipg2 
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are  each  of  the  form  given  in  (5.2.1).  We  may  decompose  the  set  U  into  three  disjoint 
sets,  U  =  Ua  [jL[jUp,  where  Ua  G  Ua,  Up  G  Ub,  and  L  =  d  Uaf]d  Up,  as  shown  in 
Figure  5.7.  We  have  already  shown  that  the  transition  maps  ips1  o  ipj1  and  0$2  °  V0* 
are  analytic  on  each  of  the  sets  Ua  and  Up.  These  transition  maps,  though,  are  also 
clearly  continuous  on  all  of  U.  Since  the  image  of  L  C  U  in  the  plane  has  measure  zero, 
any  map  which  is  A-quasiconformal  on  the  image  of  U  \  L  is  also  A-quasiconformal 
on  L  [33];  special  cases  of  this  result  are  often  given  in  texts  on  complex  analysis  [18]. 
Thus,  since  the  transition  maps  are  analytic  (1-quasiconformal)  on  the  images  of  UQ 
and  Up,  the  transition  maps  are  analytic  on  all  of  U. 


Figure  5.7:  Partition  of  U  into  Ua,  Up,  and  L 

Case  V:  Suppose  that  U  G  Us-  Then,  U  =  U^f]U^2,  where  U^ ,  U&  G  LA. 
Associated  with  each  of  U^  and  U^2  we  have  coordinate  charts  0£1,'0£2  G  'Fh-  Say 
=  S0  o  0^0  P  and  0^2  =  S&  o  g^2  o  ip^2  o  P.  Note  that  since  the  geodesic 
intersected  by  each  of  U^  and  U&  is  the  same,  we  have  =  S&  and  =  g^2.  Now, 
consider  the  transition  maps  o  and  °  0^/- 

06  ° 002X  =  (5&  oflfi  ° %  °  f)  °  (S&  °  %  °  i*)-1 

=  ■%  o  g^  ocp^oPo  P-1  o  cp-1  o  0"1  o  Sg1 
=  %  °  9b  o  (fib  o  ip-1  o  g-1  o  S’-1 
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The  maps  and  <p&  were  borrowed  from  the  conformal  structure  on  T,  so  o  ip7* 
is  analytic.  The  action  of  the  pre-composition  and  post-composition  by  grafting  and 
shearing  maps  has  no  effect  on  the  analyticity;  these  maps  each  essentially  undo  the 
local  behavior  of  the  other.  Thus,  the  analyticity  of  the  transition  map  ip^  o  ip7*  is 
shown.  Similarly,  tp^  o  xp~^ 

Thus,  the  transition  maps  in  the  structure  ( U ,  d>)  are  analytic.  □ 

Corollary  5.2.1.  Complex  earthquakes  are  maps  on  the  Teichmuller  space  of  tori. 

Proof.  Proposition  5.2.1  guarantees  that  the  image  T  of  a  torus  T  under  a  complex 
earthquake  is  a  new  torus.  That  is,  the  earthquake  has  mapped  a  point  in  the 
Teichmuller  space  of  tori  to  a  new  point  in  that  space.  □ 

5.3  Density  of  Euclidean  Earthquakes 

Having  shown  that  given  a  point  T  in  the  Teichmuller  space  of  tori  and  a  finite 
complex  Euclidean  earthquake  on  that  surface,  the  earthquake  produces  a  new  surface 
T  in  the  Teichmuller  space  of  tori,  we  turn  to  Problem  5.3.1. 

Problem  5.3.1.  Is  there  an  analogue  to  Thurston’s  Earthquake  Theorem  for  tori? 
That  is,  can  we  produce  any  surface  in  the  Teichmuller  space  of  tori  by  applying  an 
earthquake  to  the  torus  T?  If  so,  is  it  possible  to  do  so  under  the  restrictions  of 
Thurston’s  theorem,  using  exactly  one  class  of  these  transformations  (positive  shear¬ 
ing,  negative  shearing,  positive  grafting,  negative  grafting)?  Finally,  is  the  earthquake 
in  some  sense  unique? 

Taking  the  components  of  Problem  5.3.1  one  at  a  time,  in  Theorem  5.3.1  we  show 
that  we  may,  in  fact,  apply  an  earthquake  to  a  given  torus  T  and  obtain  as  the  image 
any  given  surface  in  the  Teichmuller  space  of  tori. 

Theorem  5.3.1.  Let  T  and  T  be  arbitrary  fixed  points  in  the  Teichmuller  space  of 
tori.  There  exists  a  finite  lamination  (£,  a)  on  T  and  an  associated  earthquake  E 
onT  induced  by  this  lamination  such  that  E(T)  =  T.  Moreover,  there  also  exists  a 


94 


sequence  of  measured  geodesic  laminations  {(£fc,  crk)}™=ko  on  T  and  earthquake  maps 
mz* o  on  T  induced  by  these  laminations  such  that  Ek(T)  — >•  T  as  k  — >  oo  and 
the  generators  of  the  surfaces  Ek(T )  converge  uniformly  to  the  canonical  generators 
off. 

Proof.  Let  T  and  T  be  arbitrary  fixed  points  in  the  Teichmiiller  space  of  tori.  In  this 
space,  we  normalize  so  that  one  generator  for  the  fundamental  group  of  all  tori  in 
the  group  is  z  +  1.  Thus  the  fundamental  regions  for  the  surfaces  T  and  T  are  each 
determined  by  the  complex  numbers  u  and  u  in  H,  respectively,  as  shown  in  Figure 
5.8  [29,  30]. 


Im(z) 


Re(z) 


(a)  Fundamental  Region  for  T 


Im(z) 


Figure  5.8:  Fundamental  Regions  in  the  Normalized  Teichmiiller  Space  of  Tori 


To  demonstrate  the  result  we  explicitly  describe  a  sequence  of  finite  measured 
geodesic  laminations  {(£*,  Ofc)}fc'U0  on  the  fundamental  region  for  T.  We  construct 
these  laminations  so  that  the  image  of  T  under  each  associated  earthquake  is  the 
surface  T  and  the  sequence  of  images  of  T  have  generators  converging  uniformly  to 
the  canonical  generators  of  T.  For  each  k  >  k0  G  N,  let 

=  {Lk})= i  =  +  *)  +  :  0  <  t  <  x}  • 

Examples  of  such  laminations  are  shown  in  Figure  5.9. 

To  complete  the  construction  of  the  finite  measured  geodesic  lamination,  we  need 
only  define  the  measure  <7fe.  Let  =  |(Re(w)— Re(w)),  and  let  =  |(Im(d))— Im(cj)). 
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(a)  Ci  on  T 


(b)  £2  on  T 


(c)  £5  on  T 


Figure  5.9:  Laminations  on  a  Fundamental  Region  for  a  Torus 


We  will  require  that  k  be  large  enough  that 

si  <  1  (5.3.1) 

and 

sl>- 

Clearly,  we  can  always  chose  k  large  enough  that  these  relations  are  satisfied,  since 
slsl^o  as  k  — >•  00.  Let  k0  be  the  smallest  natural  number  such  that  (5.3.1)  and 
(5.3.2)  are  satisfied.  Now,  let  k  >  A:0  and  define  <7*  :  Ck  -4  C  by  (Jk{Ljk)  =  Sk  +  iSk 
for  each  j  —  1, 2, . . . ,  k.  The  effect  of  the  constraints  given  in  (5.3.1)  and  (5.3.2)  is 
two-fold.  First,  inequality  (5.3.1)  ensures  that  we  never  shear  so  far  that  the  resulting 
fundamental  region  is  “split,”  and  no  shearing  will  twist  the  surface  a  full  revolution. 
This  requirement  is  certainly  not  necessary,  since  shearing  lengths  greater  than  one 
will  simply  indicate  that  the  generators  change  homotopy  classes.  We  impose  the 
restriction  here  for  convenience  and  clarity  of  the  construction.  The  second  constraint 
is  similar  in  its  effect,  ensuring  that  the  fundamental  regions  remain  contiguous. 
Notice  that  if  0  <  Im(a>)  <  Im(w),  and  the  grafting  actions  defined  by  the  lami- 


Im(w) 


Re(cu) 


(5.3.2) 
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nations  thus  described  are  negative,  then  for  j  =  1, 2, . . . ,  k  —  1,  we  have 

(iriIm(w)  -  _  ^  =  _  - Im(w)) 

>  — -^(Im(w)  -  Im(o;)) 

K 

=  -^(Imfij)  —  Im(w)) 
k 

>  0. 

Thus  the  weights  have  been  constructed  so  that  the  earthquake  is  always  well-defined. 
That  is,  we  avoid  the  complications  which  arise  if  too  much  of  the  surface  is  removed 
in  the  execution  of  a  grafting  associated  with  a  negative  imaginary  weight;  removing 
“too  much”  may  interfere  with  the  action  of  other  geodesics  (by  removing  them)  or 
reduce  the  torus  to  a  degenerate  surface. 

We  have  defined  a  sequence  {(£*.,  <7fc)}£lfco  of  finite  measured  geodesic  laminations 
on  the  fundamental  region  for  T.  Associated  with  each  lamination  (£fc,  ak)  we  have  an 
earthquake  Ek.  Applying  the  finite  earthquakes  associated  with  these  laminations,  we 
obtain  a  sequence  of  tori,  j  Tk  }  ;  more  precisely,  we  have  a  sequence  of  fundamental 
regions,  each  determining  a  torus.  An  illustration  of  the  action  of  the  earthquakes  on 
T  determined  by  the  laminations  described  is  shown  in  Figure  5.10,  where  we  have 
k0  =  2.  The  resulting  regions  are  clearly  fundamental  regions  for  the  image  surface 
T,  so  the  first  part  of  the  theorem  is  shown. 


Im(z)  Im(z)  Im(z) 


(a)  T^T2  (b)  T  ^  f3  (c)  T  ^  f4 

Figure  5.10:  Convergent  Euclidean  Earthquakes  on  a  Torus  T 
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Now,  as  points  in  the  (normalized)  Teichmiiller  space  of  tori  the  sequence  {t*,} 

t  J  k=ko 

is  a  constant  sequence.  We  wish  to  show  that  generators  in  this  sequence  of  tori,  which 

/■N, 

are  not  constant,  converge  uniformly  to  the  canonical  generators  of  the  given  torus  T, 
z+1  and  z  +  h.  Define  a  “distance”  between  the  surface  J5*(T )  =  Tk  and  E(T)  =  T 
in  terms  of  their  respective  fundamental  regions  by 


d(Ek(T),T)  =  sup  \Ek(r)  —  E(r)\ 

re  T 


=  max\Ek(r)  -  E(r)\ 

=  bJmr  - 


From  this,  uniform  convergence  of  the  fundamental  regions  as  k  — y  oo  is  clear  since 


lim  (  sup  \Ek(r)  —  E(t)\ 

Ko<k — ^oo  \t£T 


=  ^%/3(Irn(w)  -  Im(a>))2  -  (Re(w)  -  Re(a;))2 


=  0. 

Thus,  the  fundamental  regions  for  the  tori  converge  uniformly  to  the  fun- 

l  J  k—ko 

damental  region  of  T  defined  by  the  canonical  generators  z  +  1  and  z  +  ui;  that  is,  the 
generators  converge  uniformly  to  the  canonical  generators  of  T  as  desired.  Thus,  the 
result  is  shown.  □ 


We  notice  that  in  the  construction  used  to  prove  Theorem  5.3.1  we  had  a  sequence 
of  distinct  earthquakes  on  the  torus  T,  each  producing  T  as  the  image  of  T  under 
the  earthquake  action.  This  leads  immediately  to  the  fact  that  the  earthquakes 
transforming  one  surface  into  another  are  not  unique  in  the  Euclidean  case  as  they 
were  in  the  hyperbolic  case,  when  we  allow  all  of  left  and  right  shearing  and  positive 
and  negative  grafting.  The  construction  used  in  Theorem  5.3.1  requires  both  real 
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and  imaginary  parts  (shearing  and  grafting  maps)  taken  from  both  the  positive  and 
negative  reals.  In  order  to  translate  u  to  Cj  (for  arbitrary  w)  with  shearing  and  grafting 
maps,  we  must  have  the  freedom  to  move  along  both  axes  of  translation  and  in  both 
directions  along  each.  We  are  thus  led  to  Conjecture  5.3.1. 

Conjecture  5.3.1.  Let  two  arbitrary  points,  T  and  T,  in  the  Teichmuller  space 
of  tori  be  given.  There  does  not  exist,  in  general,  a  sequence  of  measured  geodesic 
laminations  {(Ck,  crk)}^=1  ofT  and  earthquake  maps  {Ek}°k °=1  on  T  induced  by  these 
laminations  such  that  Ek(T)  — »  T  as  k  — »  oo,  where  ok  :  A,  —5 >  A  C  C,  for  A  equal  to 
any  one  of 


{z  £  C  :  Re(z)  >  0,  Im(z)  =  0} , 

(5.3.3) 

{26C:  Re(z)  <  0,  Im(z)  =  0}  , 

(5.3.4) 

{z  £  C  :  Re(z)  —  0,  Im(z)  >  0}  ,  or 

(5.3.5) 

{z  E  C  :  Re(z)  =  0,  Im(z)  <  0} 

(5.3.6) 

for  every  k  £  N. 

This  conjecture,  together  with  Theorem  5.3.1  and  the  failure  of  uniqueness  in  the 
earthquakes  used  to  transform  one  torus  into  another  completes  our  consideration  of 
Problem  5.3.1. 

As  an  aside,  an  interesting  element  of  the  proof  is  the  constant  ko  chosen  to  satisfy 
the  inequalities  (5.3.1)  and  (5.3.2).  Defined  as  it  is  by  the  relationship  between  the 
real  and  imaginary  parts  of  u  and  w,  this  constant  depends  upon  the  Teichmuller 
distance  between  T  and  T.  That  is,  k0  provides  an  easily  calculated  indication  of  the 
Teichmuller  distance  between  the  two  surfaces. 

5.4  Discrete  Complex  Earthquake  Maps  on  Tori 
Since  Euclidean  complex  earthquake  maps  on  tori  are  easily  constructible,  we  do 
not  have  the  pressing  need  for  a  numerical  method  to  reproduce  their  actions  as  in 
the  case  of  complex  earthquakes  on  hyperbolic  surfaces.  In  the  interest  of  a  complete 
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description  of  earthquakes  on  compact  Riemann  surfaces  and  a  complete  character¬ 
ization  of  the  action  and  convergence  of  discrete  earthquakes  on  those  surfaces,  we 
have  Problem  5.4.1. 

Problem  5.4.1.  Can  we  approximate  the  action  of  complex  earthquakes  on  tori  with 
discrete  conformal  maps  ( circle  packing),  as  we  did  for  earthquakes  on  hyperbolic 
Riemann  surfaces  in  Chapter  IV? 

An  (affirmative)  answer  to  the  question  thus  posed  is  given  in  Theorem  5.4.1.  This 
theorem  is  an  analogue  of  Theorem  4.5.1,  in  which  we  proved  the  convergence  of 
discrete  earthquakes  on  hyperbolic  surfaces.  Note  that  we  construct  these  discrete 
Euclidean  earthquakes  in  exactly  the  same  manner  as  we  constructed  discrete  earth¬ 
quakes  on  hyperbolic  surfaces. 

Theorem  5.4.1.  Let  T  be  a  compact  Euclidean  Riemann  surface  of  genus  g  =  1  (a 
torus)  with  a  finite  measured  geodesic  lamination  (£,  a).  Let  {Pfc}  be  a  sequence  of 
finite  bounded  degree  packings  with  mesh  decreasing  to  zero  corresponding  to  Riemann 
surfaces  {Tk}  such  that  Tk  —>  T  as  k  -*•  oo  in  the  Teichmuller  metric.  Then  the 
surfaces  Ek{Tk )  =  Tk  induced  by  the  discrete  Euclidean  earthquake  maps  Ek  converge 
to  the  surface  E(T)  =  T  under  the  Euclidean  earthquake  map  E  induced  by  (£, a). 

Proof.  Let  T  be  a  compact  Euclidean  Riemann  surface  (a  torus),  and  let  (£,  a)  be 
a  finite  measured  geodesic  lamination  on  T.  Note  that  requiring  the  mesh  of  the 
sequence  of  bounded  degree  packings  {P*,}  to  approach  zero  as  k  — >■  oo  is  equivalent 
to  requiring  that  the  radii  of  the  circles  in  the  packings  approach  zero.  Thus,  we 
are  guaranteed  the  existence  of  such  a  sequence  by  Corollary  3.5.1.  Further,  the 
construction  of  the  surfaces  in  the  proof  of  Corollary  3.5.1  allows  the  freedom  to 
control  the  degree  of  vertices  adjacent  to  the  geodesics  along  which  we  shear  and 
graft. 

The  Teichmuller  distance  between  the  packable  surfaces  Tk  and  T  is  going  to  zero, 
so  for  each  k  e  N  there  exists  a  map  fk  :  Tk  — >■  T  that  respects  the  marking  on  T 
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and  is  1  +  efc-quasiconformal,  where  ek  — >  0  as  k  — >  oo,  (i.e.,  as  k  — »  oo,  fk  — >  /, 
where  /  is  conformal).  Note  that  we  choose  a  consistent  set  of  generators  (markings) 
for  each  surface,  so  that  the  Teichmuller  space  is  well  defined.  More  precisely,  the 
map  fk  consists  of  a  map  from  the  Euclidean  plane  to  itself  through  a  map  between 
the  surfaces  and  their  respective  conformal  structures.  Where  no  ambiguity  arises, 
however,  we  will  speak  about  the  map  as  fk  :Tk  T. 

Let  {cj}™=1,  be  the  cylinders  grafted  into  the  Riemann  surface  T  by  the  finite 
earthquake  E.  For  each  j  =  1,2,..., m,  let  {cjfc}fceN  be  a  sequence  of  surfaces  admit¬ 
ting  a  circle  packing  of  finite  uniformly  bounded  degree  so  that  Cjk  — >  Cj  as  k  — »  oo 
in  the  Teichmuller  metric.  Further,  require  that  the  mesh  of  the  packing  on  Cjk  goes 
to  zero  as  k  — >  oo  for  each  j  =  1,2 ,...,m.  Thus,  for  each  k  £  N  there  exists  a 
map  fjk  :  cjk  — Y  Cj  that  is  1  +  -quasiconformal ,  where  e.Jk  — »  0  as  k  — >  oo,  (i.e.,  as 
k  -*  oo,  fjk  — >  fj ,  where  fj  is  conformal). 

Let  T  be  the  image  of  T  under  the  action  of  the  earthquake  E,  and  for  each  k  £  N 
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let  Tk  be  the  image  of  Tk  under  the  action  of  a  combinatorial  earthquake  Ek.  We 
construct  these  combinatorial  earthquakes  by  describing  laminations  on  the  surfaces 
Tk  corresponding  to  the  lamination  on  T.  Let  k  £  N.  Corresponding  to  each  geodesic 
on  Ll  on  T,  we  choose  a  geodesic  L\  on  Tk  in  the  same  homotopy  class  as  U  such 
that  the  Hausdorff  distance  between  L\  and  fkl{Ll)  is  minimized.  For  k  sufficiently 
large,  the  geodesics  {Lj:  }"=1  thus  constructed  will  be  pairwise  disjoint.  Taking  the 
collection  of  these  geodesics,  we  have  a  finite  measured  geodesic  lamination  (£fc,  a) 
on  Tk.  Notice  that  we  have  left  the  weights  a  unchanged;  we  simply  assign  to  each 
geodesic  in  Ck  the  weight  on  the  corresponding  geodesic  in  T.  Now,  we  have  finite 
measured  geodesic  laminations  on  each  Tk,  k  £  N,  and  we  may  carry  out  the  discrete 
shearing  and  welding  operations  of  combinatorial  earthquakes.  We  also  construct  a 
set  of  secondary  geodesics  on  Tk  describing  the  second  side  of  cylinders  to  be  removed 
by  a  negative  grafting  action.  On  the  torus  T  there  is  a  geodesic  Ll  corresponding  to 
the  end  of  the  cylinder  to  be  removed  opposite  Ll.  This  is  illustrated  in  Figure  5.11. 
By  construction,  the  set  of  geodesics  and  secondary  geodesics  are  pairwise  disjoint, 
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and  no  geodesic  or  secondary  geodesic  is  contained  in  the  portion  of  the  torus  (to  be 
removed)  between  Ll  and  L\  We  then  construct  secondary  geodesics  L\  just  as  we 
constructed  the  geodesics  Vk. 


Cylinder  to  be 
Removed 


Figure  5.11:  Secondary  Geodesics  on  T 

Let  T  be  the  image  of  T  under  the  action  of  the  earthquake  E,  and  for  each  k  €  N 
let  Tk  be  the  image  of  Tk  under  the  action  of  the  combinatorial  earthquake  Ek. 

Around  each  geodesic  L1  and  each  secondary  geodesic  L\  i  =  1, 2, . . . ,  n,  on  T  we 
place  a  collar,  a  neighborhood  of  uniform  width,  so  that  no  two  such  neighborhoods 
intersect.  Call  these  collars  Cl  and  C*.  Since  the  geodesics  in  the  finite  lamination  are 
disjoint  by  construction,  we  know  that  there  exists  such  a  collar  about  each  geodesic. 
Similarly,  about  each  geodesic  L\  and  each  secondary  geodesic  L\,i  =  1,2, .. .  ,n  in 
the  lamination  Ck  on  Tk  there  exist  collars  Clk  and  Ck  such  that  the  collection  of  these 
collars  on  Tk  is  pairwise  disjoint.  We  require  that  the  collars  C\  and  Clk  on  Tk  have 
width  less  than  |  the  width  of  the  corresponding  collars  on  T.  Now,  since  Tk  —>  T 
in  the  Teichmviller  metric,  the  maps  between  Tk  and  T  are  going  to  conformal  maps 
and  the  markings  on  Tk  are  converging  to  the  markings  on  T.  Thus,  for  k  sufficiently 
large,  the  images  of  the  collars  Clk  and  Ck  under  the  map  fk  are  contained  in  the 
collars  Cl  and  C\  respectively,  (i.e.,  fk(Ck)  £  Cl  for  all  k  sufficiently  large).  Further, 
since  the  mesh  of  the  packings  on  Tk  goes  to  zero,  for  k  sufficiently  large,  the  discrete 
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geodesics  tk  and  i\  corresponding  to  Lk  and  Lk,  i  =  1, 2, . . . ,  n,  lie  strictly  within  the 
interior  of  the  collars  Clk  and  Ck,  respectively. 

Consider  the  images  of  the  collars  Clk  under  the  action  of  the  discrete  earthquake 
Ek .  If  the  geodesic  Lk  about  which  we  built  the  collar  Ck  has  real  weight,  and 
thus  induces  only  a  discrete  shearing  action  on  the  discrete  geodesic,  then  there  is  no 
difficulty  in  defining  what  we  will  call  a  quasicollar  on  T*,  as  the  image  of  Ck  under  Ek- 
The  image  of  the  discrete  geodesic  i\  and  the  image  of  the  geodesic  L\  are  contained 
in  this  quasicollar.  This  is  shown  in  Figure  5.12,  where  the  geodesic  Lk  is  shown  as 
a  solid  gray  curve,  the  discrete  geodesic  i\  is  shown  as  a  solid  black,  piecewise  linear 
curve,  and  the  collar  Ck  is  denoted  by  dashed  curves. 


(a)  Geodesic,  Discrete  Geodesic,  and  Collar  on  Tk  (b)  Images  under  a  Discrete  Shearing  Map 
Figure  5.12:  Geodesics,  Discrete  Geodesics,  and  Quasicollars 


The  case  in  which  the  weight  on  Lk  has  non-zero  imaginary  part  is  somewhat 
more  problematic.  The  action  of  the  earthquake  when  the  imaginary  component  is 
positive  will  divide  the  collar  into  two  disjoint  sets.  Since  the  grafting  (imaginary) 
component  of  the  earthquake  involves  the  insertion  of  a  cylinder,  however,  we  may 
extend  these  disjoint  regions  into  the  interior  of  the  cylinder  in  such  a  way  that  we 
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create  quasicollars  about  the  ends  of  the  inserted  cylinder.  Since  the  inserted  cylinder 
has  non-zero  height,  we  may  construct  these  quasicollars  so  that  they  are  disjoint. 
The  action  of  the  earthquake  when  the  imaginary  component  of  the  weight  is  negative 
removes  a  portion  of  the  surface,  but  our  construction  is  such  that  this  removal  simply 
creates  a  new  quasicollar,  and  the  image  of  the  discrete  geodesic  i\  and  the  image 
of  the  geodesic  L\  are  contained  in  this  quasicollar.  When  the  imaginary  part  of  the 
weight  is  negative,  the  action  of  the  combinatorial  earthquake  removes  the  surface 
between  the  discrete  geodesics,  removing  one  side  of  the  collar  about  each  geodesic. 
The  remaining  portions  of  each  collar  join  together  along  the  seam  created  by  joining 
the  two  discrete  geodesics,  to  create  a  quasicollar. 

Define  a  pullback  map  Pk  :  Tk  — *■  {Tk,  cik, . . . ,  cmjc}  by  letting  Pk(f)  be  the  unique 
pre-image  r  G  Tfc  of  the  point  f  e  Tk  under  the  combinatorial  earthquake  Ek.  Note 
that  since  the  degree  of  the  vertices  in  the  complex  on  Tk  is  bounded  and  since,  by 
construction,  the  angles  in  the  packing  on  Tk  are  bounded  away  from  zero,  the  map 
Pk  is  1  +  fSfc-quasiconformal.  Further,  at  any  point  isolated  from  the  combinatorial 
geodesics,  as  k  — »  oo  the  number  of  generations  of  circles  on  the  packed  surfaces 
between  that  point  and  the  combinatorial  geodesics  goes  to  infinity;  thus,  at  any  such 
point,  Lemma  3.2.4,  the  Packing  Lemma,  guarantees  that  the  dilatation  of  the  maps 
Pk  goes  to  1  as  k  — »  oo,  (i.e.,  as  k  ->  oo,  on  points  isolated  from  the  combinatorial 
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geodesics,  Pk  —*■  P,  where  P  is  conformal). 

We  showed  in  Theorem  5.2.1  that  T  is  a  Riemann  surface.  Thus,  since  each  of 
T,  T,  Tk,  and  Tk  are  Riemann  surfaces,  they  each  have  a  conformal  structure;  in  the 
case  of  Tk  and  Tk,  they  inherit  their  conformal  structures  from  their  packings.  Let 
Wk and  {‘lPkJkveNk  be  the  collection  of  coordinate  charts 
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in  the  conformal  structure  on  T,  T,  Tk,  and  Tk,  respectively.  Note  that  the  structure 
on  T  will  be  the  structure  described  in  Section  5.2.  Now,  we  consider  a  collection  of 
maps  fk  :  Tk  — >  T ,  each  defined  on  open  sets  in  Tk.  Let  U  C  Tk. 

yy 

Case  I:  Suppose  U  C  Tk  is  an  open  set  such  that  U  does  not  intersect  the  interior  of 
any  quasicollar,  and  U  does  not  intersect  any  cylinder  inserted  by  the  earthquake  Ek- 
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Associated  with  this  open  set  in  Tk  we  have  a  coordinate  chart  ipki/  from  the  conformal 
structure  on  Tk.  Now,  the  map  P_1  o  fkoPko  ipk  1  takes  an  open  subset  of  the  plane 
C  corresponding  to  the  pair  (U,ipkv)  to  an  open  subset  Ua  C  R  so  that  Ua  €  Ua ,  as 
described  in  Section  5.2.  Corresponding  to  Ua,  we  take  from  the  conformal  structure 
on  T  described  in  Section  5.2  a  coordinate  chart  tpQ.  Thus  we  have  the  following: 

°  P~l  °fk°Pk°  =  <Pa°  P°  P'1  °  fk°  Pk° 

=  <fa  O  fk  O  Pk  O  -0"1 . 

Since  <pa  and  are  conformal  (or  1-quasiconformal)  and  the  maps  fk  and  Pk  are 
(1  +  £fc)-quasiconformal  and  (1  +  ik  )-quasiconformal ,  respectively,  the  map  thus  de¬ 
scribed  is  (l+£fc)(l+£fc)-quasiconformal,  by  Lemma  2.4.1.  Note  that  ek  and  ik  are  in¬ 
dependent  of  the  set  U,  since  they  depend  only  on  the  maps  fk  and  Pk,  which  are  glob¬ 
ally  (1  +  £fc)-quasiconformal  and  (1  +  £*,)  -  qu  asi  conform  al ,  respectively.  By  construc¬ 
tion,  U  is  isolated  in  Rk  from  the  images  of  the  geodesics  defining  the  earthquake  Ek, 
and  £k,ik  ->  0  as  k  — >  oo.  Thus,  as  k  ->  oo,  the  map  ipa° P~  ° fk° Pk° <Pk„  :  Rk  R 
approaches  a  conformal  map. 

Case  II:  Suppose  U  C  Tk  is  an  open  set  such  that  U  does  not  intersect  the  interior 
of  any  quasicollar,  and  U  is  contained  in  the  interior  of  some  cylinder  inserted  by  the 
earthquake  Ek.  Associated  with  this  open  set  in  Tk  we  have  a  coordinate  chart  xpk„. 
Now,  the  map  P~1°fjk°Pk0'^k}  takes  an  open  subset  of  the  plane  C  corresponding  to 
the  pair  (U,ipkl/)  to  an  open  subset  Up  C  R  so  that  Up  6  Ub,  as  described  in  Section 
5.2.  Corresponding  to  Up,  we  take  from  the  conformal  structure  on  T  described  in 
Section  5.2  a  coordinate  chart  ipp.  Thus  we  have  the  following: 

■0/3  o  P-1  o  fjk  oPko  =  gjo^joFo  P-1  o  /Jfc  oPko 

=  9p  °  ftp  °  4  oPk° 

Since  and  tpk*  are  conformal  (or  1-quasiconformal)  and  the  maps  4  an(f 

Pk  are  (1  +  £jk  )-quasiconformal  and  (1  +  £/c)-quasiconformal,  respectively,  the  map 
thus  described  is  (1  +  £jk)(l  +  £fc)-quasiconformal.  Again,  we  note  that  the  constants 
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describing  the  deviation  from  conformality,  £jk  and  e^.are  independent  of  the  choice 
of  U.  By  construction,  U  is  isolated  in  Tk  from  the  images  of  the  geodesics  defining 
the  earthquake  Ek,  and  £jk,£k  ~ >  0  as  k  — >  oo.  Thus,  as  k  —>  oo,  the  map  ipp  o  o 
fjk  °  Pk  o  (pk^  :  Tk  — >■  T  approaches  a  conformal  map. 

Case  III:  Suppose  that  U  is  a  small  open  subset  of  Tk  such  that  U  intersects 
Ek{tk)  for  some  i  —  1, 2, . . . ,  n,  and  U  intersects  the  boundary  of  the  quasicollar  Clk 
associated  with  L\.  Corresponding  to  U  we  have  a  coordinate  chart  ipku-  Now,  we 
have  three  cases  to  consider  depending  on  the  imaginary  part  of  the  weight  on  the 
geodesic  Lb 

Case  Ilia:  Suppose  Im(a(L1))  =  0;  that  is,  the  earthquake  action  on  the  geodesic 
L1  involves  only  a  shearing  action.  This  gives  that  the  map  Pk  takes  U  to  a  subset  of 
Tk  such  that  Pk{U )  C{  and  Pk(U)  [\Pk  ^  0,  as  shown  in  Figure  5.13. 


(a)  An  Open  Set  U  C  Tk 


(b)  Image  of  U  Under  Pk 


Figure  5.13:  Pullback  of  Geodesics,  Discrete  Geodesics,  and  Quasicollars 


Note  that  in  Figure  5.13b  we  have  included  the  refining  edges  associated  with  the 
shearing  so  that  there  is  a  bijection  between  individual  triangles  in  Figure  5.13a  and 
Figure  5.13b. 

Now,  we  wish  to  send  this  set  Pk{U)  to  a  subset  of  T,  but  there  is  a  geometric 
difficulty  to  be  resolved.  The  “halves”  of  Pk{U )  have  as  one  portion  of  their  boundaries 
a  segment  of  the  discrete  geodesic  t\.  We  need  these  segments  of  their  boundaries  to 
lie  on  the  actual  geodesic  L\  in  Tk.  The  problem  is  easily  resolved,  at  the  expense  of 
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admitting  some  quasiconformality  in  our  eventual  transformation,  by  the  application 
of  a  Euclidean  projection  pk,  as  illustrated  in  Figure  5.14. 


Figure  5.14:  Projection  of  Geodesics,  Discrete  Geodesics,  and  Quasicollars 


This  pi  map  on  the  set  Pk(U)  is  1  +  e^-quasiconformal,  but  Proposition  4.2.1  guar¬ 
antees  that  p\  converges  uniformly  to  the  identity;  so  ek  — >■  0  as  k  -»  oo.  Further, 
the  constant  ek  is  independent  of  the  choice  of  U.  This  independence  is  a  result  of 
two  facts.  First,  the  projection  associated  with  any  given  geodesic  Ll  and  the  ac¬ 
companying  discrete  geodesic  i\  has  a  quasiconformality  constant  uniform  across  the 
entire  geodesic.  Second,  there  are  only  finitely  many  geodesics  for  which  we  must  use 
these  projections.  Taking  the  maximum  such  quasiconformal  constant  as  ek,  we  have 
independence  of  the  choice  for  U. 

Corresponding  to  the  set  pk(Pk(U))  C  Tk,  or  more  precisely  to  some  open  subset  of 
Tfc  containing  plk(Pk(U))  C  T/t,  we  have  a  coordinate  chart  pkv  which  takes  pk(Pk(U)) 
to  a  subset  of  a  fundamental  region  for  Tk  in  the  plane  C.  The  image  ofplk(Pk(U))  fj  Lk 
under  this  map  is  thus  an  arc  of  a  Euclidean  geodesic  (a  line)  in  C.  There  exists  a 
conformal  map  (comprising  translation, rotation,  and  scaling)  Mk  :  C  — »  C  such 
that  Mk  {pkv(p\(Pk{U))))  f]Lk)  is  a  small  arc  of  a  Euclidean  geodesic  in  C  that  is 
the  projection  of  the  geodesic  D  in  T.  We  may  choose  this  arc  of  the  projection 
of  Ll  sufficiently  small  that  the  set  Mk  (pkv(plk(Pk(U))))  is  contained  in  a  single 
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fundamental  region  for  T. 

A  second  difficulty  arising  from  the  projection  operation  must  be  addressed.  If 
we  were  to  apply  a  shearing  map  on  the  lift  of  the  geodesic  Ll  which  separates  the 
“halves”  of  Mk {y>kv{p\{Pk{U))))  using  the  weight  associated  with  L\  the  sides  will 
not  necessarily  match  along  the  geodesic,  as  illustrated  in  Figure  5.15. 


Figure  5.15:  Misalignment  under  Pullback,  Projection,  and  Shearing 

The  difference  along  the  geodesic  is  a  quasisymmetry  on  one  half  of  the  intersection 
with  L\  however,  and  extends  to  a  1  +£fc-quasiconformal  map,  in  the  correspond¬ 
ing  half  of  (<pkv(p%{Pk{U))))  [24,  33].  Recall,  however,  that  the  quasisymmetry 
and  associated  quasiconformal  map  arise  as  a  result  of  the  need  for  a  Euclidean  pro¬ 
jection  plk  to  correct  the  difference  between  the  geodesic  Lk  and  the  discrete  geodesic 
l\.  Lemma  4.1.1  guarantees  that  tk  — >  L\  uniformly  as  k  — >  oo.  Thus,  ek  0  as 
k  —>  oo.  Further,  since  the  quasiconformality  of  Mk  depends  only  on  the  quasisym¬ 
metry  associated  with  the  projection  pk,  £k  is  independent  of  the  choice  of  U. 

Corresponding  to  the  set  Qk  (Mf  (pkv{pl{Pk{U)))))  C  C,  we  have  a  map  from  the 
conformal  structure  on  T,  ip~l,  taking  this  set  to  an  open  set  on  T  intersecting  the 
geodesic  V.  (Again,  the  coordinate  chart  really  comes  from  an  open  set  containing 
the  set  of  interest.)  Notice  that  as  a  result  of  applying  the  correction  Qk ,  the  two 
halves  of  this  image  (on  either  side  of  the  lift  to  C  of  the  geodesic  Ll)  now  match 
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the  action  of  the  shearing  caused  by  P  1.  We  now  apply  P  1,  to  obtain  a  set  V  C 
T.  Corresponding  to  this  open  set  we  have  a  coordinate  chart  from  the  conformal 
structure  on  T  defined  in  Section  5.2.  Since  V  G  Ur,  as  described  in  Section  5.2,  we 
have  a  coordinate  chart  corresponding  to  V.  This  gives  the  following: 

^7  °  F_1  °  Vv1  °Qk°  Mk  °  Vkv  o  p\  oPkO 

=  S1o^p1oPoP-lo  ip~l  oQ%o  M%  O  ipkv  o  p\  oPko  <p~l 
=  s^ocp^o  ip-1  oQ"o  Mg  O  <pkv  op\  oPko  (p-\ 

This  map  is  (1  +  e*)(l  +  £fc)(l  +  ) -quasiconfor  mal ,  where  this  quasiconformality  is 

independent  of  the  choice  for  U.  Now,  consider  the  dilatation  at  a  point  u  G  17.  If  u 
is  not  contained  in  the  set  Ek(£k ),  then  Lemma  3.2.4,  the  packing  Lemma,  guarantees 
that  the  dilatation  at  u  goes  to  1  since  each  of  £k,ek,ik  — >  0  as  k  — >  oo.  Thus,  this 
map  from  Tk  to  T  converges  to  a  conformal  map,  except  perhaps  on  the  image  of  the 
discrete  geodesic  l\  under  the  earthquake  Ek.  But  this  set  has  measure  zero,  and  the 
map  is  continuous.  Thus,  the  limit  map  is  conformal  on  all  of  U  [33]. 

Case  Illb:  Suppose  Im(a(L1))  >  0;  that  is,  the  earthquake  action  on  the  geodesic 
Ll  is  a  positive  grafting  action  (and  potentially  a  shearing  action  as  well).  The 
situation  here  is  somewhat  more  complicated  than  the  previous  case  since  the  map 
Pk  takes  the  portion  of  U  on  the  inserted  cylinder  to  a  cylinder  Cjh  and  the  remaining 
points  of  U  to  a  subset  of  Tk.  We  will  deal  with  the  maps  on  these  two  subsets  of  U 
separately. 

First,  consider  that  subset  of  U  which  is  taken  to  Tk.  Call  this  subset  U\.  Pk(Ui) 

C'l  and  Pk(Ui)  f]£k  0.  Now,  we  wish  to  send  this  set  Pk{U\)  to  a  subset  of  T,  but 
there  is  a  geometric  difficulty  to  be  resolved.  Pk{U{)  has  as  one  portion  of  its  boundary 
a  segment  of  the  discrete  geodesic  we  need  this  segment  of  the  boundary  to  lie 
on  the  actual  geodesic  Lk  in  Tk.  As  before,  this  problem  is  easily  addressed,  at  the 
expense  of  admitting  some  quasiconformality  into  our  eventual  transformation,  by 
the  application  of  a  Euclidean  projection  p\.  This  map  on  the  set  Pk{U\)  is  1  +  ek- 
quasiconformal,  but  Proposition  4.2.1  guarantees  that  pk  converges  uniformly  to  the 


109 


identity;  so  ek  -»  0  as  k  -»  oo.  As  in  Case  Ilia,  ek  is  independent  of  the  choice  of  U. 

Corresponding  to  the  set  p\(Pk(Ui ))  C  Tk,  or  more  precisely  to  some  open  sub¬ 
set  of  Tjfc  containing  p\(Pk{Ui ))  C  Tk,  we  have  a  coordinate  chart  ^  which  takes 
plk(Pk(Ui))  to  a  subset  of  a  fundamental  region  for  Tk  in  the  plane  C.  The  image  of 
Pk(Pk(Ui))  f]Ll  under  this  map  is  thus  an  arc  of  a  Euclidean  geodesic  in  C.  There 
exists  a  conformal  map  (comprising  a  translation,  rotation,  and  scaling)  M^1  :  C  — >  C 
such  that  M^(<pkv(pi(Pk(Ui))))riLi)  is  an  arc  of  the  Euclidean  geodesic  in  C  that 
is  the  lift  of  the  geodesic  Ll  in  T.  We  may  choose  this  arc  small  enough  that  the  set 
Mkl {pkv{Pk(Pk{Ui))))  is  contained  in  a  single  fundamental  region  for  T. 

Again,  we  will  require  a  correction  to  compensate  for  the  effect  of  the  Euclidean 
projection  p\.  Call  this  correction  Qk,  and  note  that  this  map  is  a  quasisymmetry 
on  the  intersection  with  the  lift  of  the  geodesic  U  and  is  (1  +  Sfc)-quasiconformal,  as 
in  Case  Ilia.  This  quasiconformality  constant  e*,  is  independent  of  the  choice  for  U. 
Corresponding  to  the  set  Q% (Mf1(^fcu(p|(Ffc([/1)))))  c  C,  we  have  a  map  from  the 
conformal  structure  on  T,  ip~l,  taking  this  set  to  an  open  set  on  T  intersecting  the 
geodesic  L\  (Again,  the  coordinate  chart  really  comes  from  an  open  set  containing 
the  set  of  interest.)  We  now  apply  P  l,  to  obtain  a  set  Vj  C  T  adjoining  the  edge  of 
an  inserted  cylinder  Cj. 

Now,  consider  that  subset  of  U  which  is  taken  into  Cjk.  Call  this  subset  U2. 
Pk{U2)  C  Cjk  and  Pk(U2)  intersects  the  boundary  of  the  cylinder  Cjk.  We  take  this  set 
into  the  cylinder  Cj  via  the  1  +  e^-quasiconformal  map  fjk ,  where  ek  is  independent  of 
the  choice  of  U,  since  we  may  choose  ek  as  the  maximum  of  the  constants  associated 
with  the  finitely  many  cylinders  {c,}.  We  then  apply  the  map  P-1  to  obtain  a  set 

v2cf. 

The  way  we  have  constructed  these  maps,  in  particular  the  application  of  the 
map  Qk,  ensures  that  the  set  V  =  Vi\jV2  C  T  is  an  element  of  the  collection  IJ& 
as  described  in  Section  5.2.  Corresponding  to  this  set  V  we  have  a  coordinate  chart 

from  the  conformal  structure  on  T  given  in  Section  5.2.  This  gives  the  following 
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piecewise  map  Uk  :  Tk  ^  T: 


n*(«) 


^  $po  P  1oiPv1oQ^o  Mk'  O  ipkv  otfkoPko  Ipk'(u)  if  u  <E  Pi  C  Tk 
tip  °  -P'1  °  fjk  °Pk°  ^k!  (u)  if  u  e  U2  C  ffc. 


The  first  of  these  maps,  for  u  €  Pi,  is  (1  +  £fc)(l  +  £fc)(l  +  £fc)-quasiconformal.  The 
second  map,  for  u  €  C/2,  is  (1  +  £jj(l  +  £fc)-quasiconformal.  As  we  have  shown,  these 
constants  are  independent  of  the  choice  of  U.  Now,  consider  the  dilatation  at  a  point 
■u  G  P.  If  u  is  not  contained  in  the  set  Ek(tk),  then  Lemma  3.2.4,  the  packing  Lemma, 
guarantees  that  the  dilatation  at  u  goes  to  1  since  each  of  £k,£k,ik-,  Sjk  — »  0  as  k  —¥  oo. 
Thus,  this  piecewise  map  life  from  Tk  to  T  converges  to  a  conformal  map  II,  except 
perhaps  on  the  image  of  the  discrete  geodesic  t\  under  the  earthquake  Ek.  But  this 
set  has  measure  zero,  and  the  map  is  continuous.  Thus,  the  map  II  is  conformal  on 
all  of  U  [33], 

Case  IIIc:  Suppose  Im(er(I/))  <  0;  that  is,  the  earthquake  action  on  the  geodesic 
D  is  a  negative  grafting  action  (and  potentially  a  shearing  action  as  well).  Here,  the 
map  Pfc  divides  the  set  U  into  two  disjoint  sets,  separated  by  an  inserted  cylinder; 
essentially,  the  first  step  in  this  action  is  the  inverse  of  the  action  described  in  Case 
Illb.  As  in  Case  Illb,  we  will  consider  these  sets  separately. 

Consider  the  subset  of  U  whose  image  borders  the  discrete  geodesic  i\\  call  this 
subset  U\.  Pk(U\)  C  Ck  and  Pk(Ui)f]Pk  /  0.  Now,  we  wish  to  send  this  set  Pfe(Li) 
to  a  subset  of  T,  but  there  is  a  geometric  difficulty  to  be  resolved.  Pfc  (Pi)  has  as  one 
portion  of  its  boundary  a  segment  of  the  discrete  geodesic  Pk,  we  need  this  segment 
of  the  boundary  to  lie  on  the  actual  geodesic  Lk  in  Tk.  As  before,  this  problem  is 
easily  addressed,  at  the  expense  of  admitting  some  quasiconformality  into  our  eventual 
transformation,  by  the  application  of  a  Euclidean  projection  p\.  This  map  on  the  set 
Pfc  (Pi)  is  1  +  efc-quasiconformal,  but  Proposition  4.2.1  guarantees  that  pk  converges 
uniformly  to  the  identity;  so  ek  — >  0  as  k  — >  oo.  As  in  Case  Ilia,  ek  is  independent  of 
the  choice  of  P. 

Corresponding  to  the  set  p\(Pk(Ui))  C  Tk,  or  more  precisely  to  some  open  sub- 
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set  of  Tk  containing  pk(Pk(Ui))  C  Tk,  we  have  a  coordinate  chart  <pkvi  which  takes 
p\(Pk{Ui))  to  a  subset  of  a  fundamental  region  for  Tk  in  the  plane  C.  The  image  of 
pk(Pk(Ui))f]Lk  under  this  map  is  thus  an  arc  of  a  Euclidean  geodesic  in  C.  There 
exists  a  conformal  map  (comprising  a  translation,  rotation,  and  scaling)  M^1  :  C  — >  C 
such  that  MlkJl  (pkvi  (p^(Pfe(171))))  f]  L\)  is  an  arc  of  the  Euclidean  geodesic  in  C  that 
is  the  lift  of  the  geodesic  U  in  T.  We  may  choose  this  arc  small  enough  that  the  set 
Mkl(tpkv  (p\(Pk{Ui))))  is  contained  in  a  single  fundamental  region  for  T. 

Again,  we  will  require  a  correction  to  compensate  for  the  effect  of  the  Euclidean 
projection  pk.  Call  this  correction  Qk,  and  note  that  this  map  is  a  quasisymmetry 
on  the  intersection  with  the  lift  of  the  geodesic  D  and  is  (1  +  e^j-quasiconformal,  as 
in  Case  Ilia.  Again,  the  quasiconformality  constant  ek  is  independent  of  the  choice 
for  U. 

Corresponding  to  the  set  Qk  (Mkx (cpkvi(p[(Pk{U i)))))  C  C,  we  have  a  map  from 
the  conformal  structure  on  T,  y?"1,  taking  this  set  to  an  open  set  on  T  intersecting  the 
geodesic  L\  (Again,  the  coordinate  chart  really  comes  from  an  open  set  containing 
the  set  of  interest.)  We  now  apply  P  ,  to  obtain  a  set  Vj  C  T. 

Now,  consider  the  subset  of  U  whose  image  borders  the  discrete  ( secondary ) 
geodesic  l\\  call  this  subset  C/2.  We  map  this  set  to  T  in  essentially  the  same  manner 
as  we  mapped  U\  to  T.  The  differences  lie  in  distinct  choices  from  the  conformal 
structures  on  Tk  and  T,  a  different  choice  for  the  conformal  map  M,  and  in  not  ap¬ 
plying  a  correction  map  Q  (since  we  need  only  correct  on  one  side  of  the  geodesic) . 
Thus  we  will  have  a  map 

P~l  o  Pvt  °  Mk2  °  'fkV2  o  p\  oPko 

which  takes  (72  to  a  set  14  C  T 

The  way  we  have  constructed  these  maps,  in  particular  the  application  of  the 
map  Qk,  ensures  that  the  set  V  =  VjU  V2  C  T  is  an  element  of  the  collection  LA 
as  described  in  Section  5.2.  Corresponding  to  this  set  V  we  have  a  coordinate  chart 

from  the  conformal  structure  on  T  given  in  Section  5.2.  This  gives  the  following 
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piecewise  map  life  :  Tk  — >  T: 


I  o  P  1  o  (pyi  o  QU  O  Ml?1  o  <pkvi  op*fe  oPko  ipkl{u)  if  u  G  Ui  C  Tk 

nfc(«)  =  < 

O  P~l  O  fjk  oPko  ipk^{u)  if  u  G  U2  C  ffc. 

The  first  of  these  maps,  for  u  G  Ui,  is  (1  +  £fc)(l  +  £jt)(l  +  Sfc)-quasiconformal.  The 
second  map,  for  u  E  C/2,  is  1  -\-ek)(l  +  £fc)-quasiconformal.  As  we  have  shown,  each  of 
these  constants  is  independent  of  the  choice  of  U.  Now,  consider  the  dilatation  at  a 
point  u  €  U.  If  u  is  not  contained  in  the  set  Ek(Pk),  then  Lemma  3.2.4,  the  packing 
Lemma,  guarantees  that  the  dilatation  at  u  goes  to  1  since  each  of  £k,ek,Sk,  £jk  — >  0 
as  k  ->  00.  Thus,  this  piecewise  map  II ^  from  Tk  to  T  converges  to  a  conformal 
map  II,  except  perhaps  on  the  image  of  the  discrete  geodesic  tk  under  the  earthquake 
Ek.  But  this  set  has  measure  zero,  and  the  map  is  continuous.  Thus,  the  map  II  is 
conformal  on  all  of  U  [33]. 

Thus,  as  A;  — >  00  the  maps  from  Tk  to  T  are  becoming  conformal,  and,  by  con- 

- - 

struction,  the  markings  on  surfaces  T  and  Tk  are  consistent.  Therefore,  the  surfaces 
Tfc  converge  to  a  surface  conformally  equivalent  to  with  equivalent  markings;  but 
this  is  equivalent  to  sayirig  %  — >  T  (in  the  Teichmiiller  metric)  as  k  — >  00,  and  the 
result  is  shown.  □ 

We  now  give  several  examples  of  finite  earthquakes,  both  shearing  and  grafting, 
on  a  compact  torus  T,  and  we  approximate  these  earthquakes  using  the  circle  packing 
techniques  given. 

Example  5.4.1.  Let  T  be  the  point  in  the  normalized  Teichmiiller  space  of  tori 
defined  by  the  generators  z  +  1  and  z  +  u,  where  w  =  |  e  H.  We  place  a 

geodesic  L  on  this  surface  such  that  the  segment  between  and  is  an 

arc  of  the  lift  of  that  geodesic  through  the  canonical  fundamental  region  for  T.  This 
fundamental  region  for  T  and  the  geodesic  arc  described  are  shown  in  Figure  5.16a. 
Place  a  weight  of  a  ( L )  =  —  ^  on  the  geodesic.  The  action  of  the  (right)  earthquake 
induced  by  this  finite  geodesic  lamination  on  T  is  shown  in  Figure  5.16b.  Note 
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/*s 

that  though  the  generators  of  the  image  surface  T  are  not  the  canonical  generators 
associated  with  the  image  point  in  the  Teichmiiller  space  of  tori,  the  position  in 
Teichmiiller  space  is  determined  by  the  point  Cj  =  1  +  E  H  [29,  30]. 


Figure  5.16:  Explicit  Shearing  Action  on  a  Torus 


Figure  5.17:  Packing,  Discrete  Lamination,  and  Earthquake  on  a  Torus 


The  torus  T  shown  in  Figure  5.16a  is  itself  a  packable  surface,  so  we  may  place 
a  packing  on  the  surface  as  shown  in  Figure  5.17a.  To  remain  consistent  with  the 
notation  of  Theorem  5.4.1,  we  will  call  this  packed  surface  though  there  is  no 
difference  between  Ti  and  T  itself.  Also  in  Figure  5.17a  we  show  the  geodesic  L  (as 
a  dashed  segment)  and  the  discrete  geodesic  approximating  L  (as  a  solid  segment). 
Now,  we  execute  the  discrete  (combinatorial)  earthquake  E\  determined  by  the  dis- 
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Crete  geodesic  and  the  weight  on  L  to  obtain  the  surface  Ti  shown  in  Figure  5.17b. 
Note  that  the  discrete  earthquake,  in  this  case  a  right  shearing  action,  respects  the 
combinatorics  of  the  packing  on  Ti,  as  described  in  Figure  3.12.  The  discrete  earth¬ 
quake  results  in  a  normalized  point  in  the  Teichmiiller  space  of  tori  determined  by 
the  point  1  +  €  H.  Thus,  the  discrete  approximation  exactly  mimics  the  action 

of  the  true  earthquake.  This  will  always  occur  when  the  packing  on  our  surface  is  as 
regular  as  the  packing  shown  (a  regular  hex  packing)  and  the  shearing  respects  the 
combinatorics  of  the  triangulation. 

Example  5.4.2.  In  this  example,  we  repeat  the  process  given  in  Example  5.4.1, 
changing  only  the  weight;  here  we  use  a  ( L )  =  |.  All  other  parameters  for  the 
surface  and  (right)  earthquake  described  in  Example  5.4.1  are  unchanged.  When  we 
execute  the  earthquake  E  induced  by  these  parameters,  we  obtain  the  surface  T  in 
the  normalized  Teichmiiller  space  of  tori  determined  by  the  point  u  —  —  -f  -~-i  €  HI. 


(a)  Packing  and  Discrete  Lamination  on  T\  (b)  E\  (Tj)  =  Ti 


Figure  5.18:  Discrete  Shearing  Map  on  a  “Coarsely”  Packed  Torus 

Again,  the  surface  T  is  itself  packable,  and  we  place  a  packing  on  the  surface  as 
shown  in  Figure  5.18a  to  obtain  the  surface  Ti,  identical  to  T.  Now,  the  seemingly 
minor  change  in  the  weight  on  L  results  in  a  combinatorial  earthquake  which  does 
not  respect  the  combinatorics  of  the  triangulation  on  T\.  Thus,  we  must  refine  the 
packing,  as  shown  in  Figure  4.3.  This  results  in  new  circles  added  to  the  packing 
when  we  execute  the  combinatorial  earthquake  E\,  which  should  alter  the  result 
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obtained  in  Example  5.4.1.  In  Figure  5.18b  we  see  that,  after  normalization,  the 
action  of  the  E\  on  the  packable  torus  Tx  results  in  a  torus  determined  by  the  point 
wi  =  1.1  4-  0.615973*  G  EL 

Since  the  combinatorial  earthquake  E1  did  not  result  in  the  desired  image  T, 
we  refine  the  packing  to  obtain  a  new  packable  surface  T2.  While  we  continue  with 
notation  consistent  with  Theorem  5.4.1,  T2  is  again  identical  to  T.  The  refinement  we 
use  is  the  hex  refinement  described  in  Section  3.3.  This  new  packing,  the  geodesic, 
and  discrete  geodesic  are  shown  in  Figure  5.19a.  Note  that  the  refinement  results 
in  the  selection  of  a  new  edge  path  for  the  discrete  geodesic  which  is  closer  to  true 
geodesic.  Now,  we  execute  the  combinatorial  earthquake  E2  on  T2  to  obtain  a  new 
torus  determined,  by  the  point  Q 2  =  1.05  +  0.654396*  G  H,  as  shown  in  Figure  5.19b. 

ImCz) 

Im (z) 


Re(z) 

(a)  Packing  and  Discrete  Lamination  on  T2  (b)  E2  (T2)  =  T2 

Figure  5.19:  Discrete  Shearing  Map  on  a  “Finely”  Packed  Torus 


Note  that  we  are  now  “closer”  to  the  desired  destination  surface  T.  That  is,  if 
we  denote  the  Teichmiiller  distance  between  two  points  S  and  S  in  the  normalized 
Teichmiiller  space  of  tori  by  dr  S^j ,  and  we  recall  that  this  distance  corresponds 
to  the  hyperbolic  distance  between  the  points  cus,  ccg  G  H  [29],  we  have  the  following: 


dT  (t,Tx)  =  dw  (w.wx)  =  dm  ^  +  2~h  1  1  +  0.615973i  j  =  0.1641895 
dT  (T,f2)  =  dn  (u>,  a>2)  =dn(^  +  1.05  +  0.654396^  =  0.06806798. 
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Thus,  one  iteration  of  hex  refinement  improves  our  approximation  by  a  factor  of 
approximately  2.41214.  In  these  computations,  we  use  the  formula 


z,  w )  =  cosh  1  ( 1  + 


‘0 


z  —  w 


(5.4.1) 


2Im(^)Im(t(;)  J  ’ 

giving  the  hyperbolic  distance  in  the  half-plane  model  between  points  z,wGl  [2]. 


Example  5.4.3.  As  in  Example  5.4.1  and  Example  5.4.2,  we  let  T  be  the  point  in 
the  normalized  Teichmuller  space  of  tori  defined  by  the  generators  z  +  1  and  z  +  to, 
where  to  —  \  +i^-  €  H.  We  place  a  geodesic  L  on  this  surface  such  that  the  segment 
between  ||  +  i- ~  and  ||  +  is  an  arc  of  the  lift  of  that  geodesic  through  the 
canonical  fundamental  region  for  T.  This  fundamental  region  for  T  and  the  geodesic 
arc  described  are  shown  in  Figure  5.20a.  Place  a  weight  of  a  ( L )  =  \i  on  the  geodesic. 
The  action  of  the  (grafting)  earthquake  induced  by  this  finite  geodesic  lamination  on  T 
is  shown  in  Figure  5.20b.  Note  that  though  the  generators  of  the  image  surface  T  are 
not  the  canonical  generators  associated  with  the  image  point  in  the  Teichmuller  space 
of  tori,  the  position  in  Teichmuller  space  is  determined  by  the  point  to  —  |  + 

[29,  30]. 


Figure  5.20:  Explicit  Grafting  Action  on  a  Torus 


The  surface  T  is  itself  packable,  and  we  place  a  packing  on  the  surface  as  shown 
in  Figure  5.21a  to  obtain  the  surface  7\,  identical  to  T. 


117 


(a)  Packing  and  Discrete  Lamination  on  T\ 


(b)  Packing  on  Inserted  Cylinder 


Figure  5.21:  Packing  and  Discrete  Lamination  on  a  Torus  and  a  Packed  Cylinder 

Now,  using  the  technique  shown  in  Figure  4.4,  we  graft  together  the  combinatorics 

in  the  packings  in  Figure  5.20a  and  Figure  5.20b.  That  is,  we  execute  the  discrete 

_ _ 

earthquake  E\  to  obtain  the  packed  torus  Ti  shown  in  Figure  5.22a.  The  torus  T\ 
is  normalized  so  it  is  uniquely  described  by  the  point  —  .499124  +  1.00916*  G  H. 
Thus  we  may  calculate  the  distance  in  the  Teichm  idler  metric  between  T  and  Ti, 
dT  (f ,  fi)  =  0.122204. 

Im(z) 


ReU) 

(a)  Ei  (Ti)  =  Ti  (b)  E2  (T2)  =  f2 

Figure  5.22:  Discrete  Grafting  Maps  on  a  Torus  T 

As  in  Example  5.4.2,  we  may  improve  our  approximation  of  T  through  a  refinement 
of  the  original  surfaces.  We  hex  refine  the  circle  packing  on  T  to  obtain  a  new 
surface  T2;  as  before,  T2  is  identical  to  the  original  torus  T,  but  we  wish  to  remain 


ImU) 


Re(z) 


consistent  with  the  notation  of  Theorem  5.4.1.  Note  that  we  have  only  refined  the 
packing  on  the  torus;  the  packing  on  the  cylinder  remains  unchanged  in  the  grafting 
earthquake  E2  on  T2.  The  image,  E2  (T2)  =  T2,  of  the  earthquake  E2  is  shown  in 
Figure  5.22b.  The  torus  T2  is  normalized  so  that  it  is  uniquely  described  by  the 
point  u>2  =  0.404856  +  0.899675i  e  H.  The  Teichmiiller  distance  between  T  and  T2 
is  given  by  dr  (T,T2^j  —  0.0471377.  Thus,  one  iteration  of  hex  refinement  improved 

**'S. 

our  approximation  of  T  by  a  factor  of  2.5925. 

5.5  Earthquakes  on  Other  Surfaces 

We  are  now  in  a  position  to  give  a  complete  characterization  of  earthquakes  on 
compact  Riemann  surfaces.  We  recall  that,  as  a  consequence  of  the  Uniformization 
Theorem  [15],  compact  Riemann  surfaces  fall  into  a  finite  number  of  categories.  We 
have  described  the  action  of  earthquakes  on  the  Teichmiiller  spaces  of  compact  tori 
and  compact  hyperbolic  Riemann  surfaces  (n-holed  tori  with  n  >  2).  Among  the 
remaining  surfaces,  only  the  sphere  is  compact.  The  Teichmiiller  space  of  the  sphere 
is  the  trivial  (one-point)  space.  We  may  thus  easily  define  an  earthquake  on  the 
Riemann  sphere  and  obtain  trivial  transformations  in  the  Teichmiiller  space  of  the 
sphere. 
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APPENDIX 


APPENDIX  A:  DRAWING  THE  LAMINATION  ON  O 

A  well  known  property  of  cross-ratios  found  in  most  complex  analysis  textbooks  (see, 
for  example,  [13])  states  that  the  relation 

(w  -  Wi)(w2  -  w3)  _  {z  -  zl){z2  -  z3)  1 

(w  -w3)(w2  -  Wi)  {z  -  z3)(z2  -  Zi) 

implicitly  defines  the  unique  Mobius  transformation  (or  fractional  linear  transfor¬ 
mation)  which  maps  distinct  points  zi,z2,a nd  z3  in  the  finite  complex  plane  onto 
distinct  points  w i,  w2,  and  w3  in  the  finite  complex  plane.  Solving  this  relation  for  w, 
we  obtain  an  explicit  expression  for  this  unique  Mobius  transformation.  This  action 
is  realized  in  the  Mathematica  code  below.  In  the  notation  of  the  code  given,  the  map 
obtained  has  the  property  that  a  —  Z\,b  =  z2,  c  =  z3,d  =  W\,e  =  w2,  and /  =  w3. 

Mobius [{a_ , b_ , c_> , {d_ , e_ , f _>] : = 

Chop [N [ 

( (ade-dbe-adf +bef +dcf-ecf ) z+ (-adec+dbec+adbf -abef-dbcf +aecf ) ) / 
( (-db+ae+dc-ec-af +bf ) z+ (adb-abe-adc+bec+acf -bcf ) ) ] ] 

The  Mathematica  code  below  takes  two  points  m  and  n  on  the  unit  circle  and 
computes  a  parametric  description  for  the  unique  hyperbolic  geodesic  between  the 
two  points.  We  do  this  by  composing  a  parametrization  of  the  diameter  between 
z  =  —  1  and  z  =  1,  q(t)  =  —  1  +  2t,  0  <  t  <  1,  and  the  appropriate  Mobius  transfor¬ 
mation  sending  this  diameter  to  the  desired  hyperbolic  geodesic,  computed  using  the 
command  Mobius []  from  above. 

Eqn[m_,n_] := 

Block [{p,q,r}, 

p=Exp [I (Min [Arg [m] , Arg [n] ] ) ] ; 
q=Exp [I (Min [Arg [m] , Arg [n] ] ) ] ; 
r=Exp [I (Arg [m] +Arg [n] ) /2] ; 

FullSimplif y [{Re [Mobius [{ 1 , I , -1} , {p , r , q>] / . z-> (-l+2t ) ] , 

Im [Mobius [{1 , I , -1} , {p , Exp [Ir] , q>] / . z-> (-l+2t ) ] >] ] 

The  following  Mathematica  function  takes  a  list  of  n  >  1  pairs  of  distinct  points 
on  the  unit  circle,  {{zi,i,  Zij2}  ,  {z2)i,  z2,2}  ,  •  •  • ,  \zn,i:  zn,2}},  and  draws  the  lamination 
consisting  of  the  geodesics  connecting  zitl  and  zij2, 1  <  i  <  n. 

Lamination [pts_] := 

ParametricPlot [Evaluate [ 

Union  [{{Cos  [2  Pi  t]  ,  Sin  [2  Pi  t]». 

Table [Eqn [pt s [ [i] ] ] , {i , 1 , Length [pt s] }] ] ] , {t ,  0,  1} , 
AspectRatio->Automat ic ,  Axes->None ,PlotRange->All] 
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APPENDIX  B:  ACTION  OF  A  FINITE  LEFT  EARTHQUAKE 

The  following  Mathematica  code  explicitly  defines,  as  a  function  of  z,  the  compar¬ 
ison  map  associated  with  the  ith  geodesic  having  left  and  right  endpoints  x  and  y, 
respectively,  and  weight  u. 

Compar isonMap  [i_ ,  {{x_ ,  y_> ,  u_} ,  z_]  :  = 

Simplify [ ( (yExp [u] -x) z+yx (1-Exp [u] ) ) / ( (Exp [u] -1) z+ (y-xExp [u] ) ) ] 

The  following  code  defines  a  piecewise  function  for  the  application  of  the  shearing 
map  for  a  left  earthquake.  Note  that  this  map  is  the  comparison  map,  above,  on  one 
domain  and  the  identity  map  elsewhere. 

M[i_,z_]:=  Which [ 

Arg  [Q  C  Ci ,  1 , 2]  ]  ]  <Arg  [Q  [  [i ,  1 , 1]  ]  ]  , 

Which [ 

Or  [Arg  [z]  <Arg  [Q  [  [i ,  1 , 2]  ]  ]  ,  Arg  [z]  >Arg  [Q  [  [i ,  1 , 1]  ]  ]  ]  , 

ComparisonMap [i , {{Q [ [i , 1 , 1] ] ,Q[[i,l,2]]},Q[[i,2]]},z] , 
True , z] , 

Arg  [Q  [  [i ,  1 , 2]  ]  ]  >Arg  [Q  [  [i ,  1 , 1]  ]  ]  , 

Which [ 

And [Arg [z] <Arg [Q [ [i , 1 , 2] ] ] , Arg [z] >Arg [Q [ [i , 1 , 1] ] ] ] , 

ComparisonMap [i , {{Q [ [i , 1 , 1] ] ,Q[[i,l,2]]>,Q[[i,2]]},z] , 
True , z] ] 

The  Mathematica  function  below  takes  as  its  input  a  list  of  points  S  C  c©  and 
a  description,  L,  of  the  finite  measured  geodesic  lamination  associated  with  S.  S  = 
{si,  s2,  •  •  • ,  sn}  is  a  hst  of  complex  numbers  and  may  be  expressed  in  either  polar  or 
rectangular  form.  L  is  a  list  of  the  form 

L  =  {{{ih,i )  Pi, 2}  ,  Ah}  ,  {{P2, i,P2$}  >  M2}  ,  •  •  • ,  {{Pn-3,1 5^71-3,2}  >  Mn- 3}}  > 

where  (p*, 1,^,2}  describes  the  ith  geodesic  in  the  tree  of  triangles,  and  fii  is  the  weight 
associated  with  the  ith  geodesic.  (Recall  that  a  finite  lamination  on  n  vertices  is  com: 
prised  of  n  —  3  disjoint,  hyperbolic  geodesics.)  In  this  notation,  pi}  1  and  pj)2, 1  <  i  < 
n  —  3,  are  the  left  and  right  endpoints  of  the  fth  geodesic,  respectively.  Note  that  this 
code  calls  the  function  M[],  above,  which  then  calls  the  function  ComparisonMap[]. 

Earthquake [S_ , L_] : = 

Block [{P=S,Q=L}, 

Do[ 

Do [P [ [k] ]  =M [ j  , P [ [k] ] ]  , -[k,l, Length [P]>]  ; 

Do [Q [ [k , 1 , 1] ] =M [ j ,Q[[k,l,l]]] ,{k,l, Length [Q]>] ; 

Do [Q [ [k , 1 , 2] ] =M [ j , Q [ [k , 1 , 2] ] ] , {k,l, Length [Q]>] ; , 

{ j,l, Length [Q]>] ; 

{P,Q >] 
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APPENDIX  C:  COMPUTING  A  MEASURED  GEODESIC  LAMINATION 


The  Complement  []  command  in  Mathematica  delivers  the  complement  of  one 
set  relative  to  another.  The  result  returned  by  the  function  is  sorted  in  canonical 
order.  The  following  function  returns  the  list  x\y  without  altering  the  original  order 
of  elements  in  x.  [55] 

UnsortedComplement [x_List ,y_List] := 

Replace [x , Dispatch [ (# : ->Sequence [] ) &/@Union [y] ] , 1] 

The  following  Mathematica  function  is  a  cross-ratio  given  in  [24]  linked  to  the 
amount  points  of  S  are  moved  to  the  left  by  a  finite  left  earthquake.  This  command 
is  used  in  computing  the  lamination  and  associated  measure. 

cr[a_,b_,c_,d_] :=((d  -  c) (b  -  a))/((c  -  b) (a  -  d)) 

Now,  given  a  set  of  points  S  =  {sl5  s2, . . .  sn}  and  their  images  h(S)  the  algorithm 
below  generates  the  unique  finite  lamination  £  and  a  non-negative  measure  a  associ¬ 
ated  with  S  such  that,  up  to  post-composition  by  a  Mobius  transformation  h  is  the 
restriction  to  S  of  the  finite  left  earthquake  ha.  This  inductive  algorithm  is  suggested 
as  an  exercise  in  [24]  and  examined  in  more  detail  in  [23]. 

FindLam[pre_,post_] := 

Block [{StartPts, ImagePts ,T,a,b,c,al,bl,cl,L,A,q,test ,i, j ,k,s,K,TTemp}, 
StartPts=pre ; ImagePts=post ;L={} ; 

While [Length [StartPts]  >=4 , 

K=Table [ 

{Mod [ i , Length [St art P t  s] ] + 1 , 

Mod [i+1, Length [StartPt s]]+l, 

Mod [i+2 , Length [StartPt s] ] +1} , {i , 0 , Length [StartPts] -1}] ; 

T=Table [ 

{{StartPts [[K[k,l]]]] , StartPts [[K[[k,2]]]] , 

StartPts [ [K [ [k , 3] ] ] ] > , 

{ImagePts [ [K [k , 1] ] ] ] , ImagePts [ [K [ [k , 2] ] ] ] , 

ImagePts [ [K [ [k , 3] ] ] ] » , {k , 1 , Length [K] >] ; 
dist= 

Table [Table [ 

Log[ 

Chop [N [ 

cr[T[[i,2,l]] ,T[[i,2,2]] ,T[[i,2,3]] , 

UnsortedComplement [ImagePts ,T[[i, 2] ]][[]]]]/ 
cr [T [ [i , 1 , 1] ] ,T[[i,l,2]] ,T[[i,l,3]] , 

UnsortedComplement [StartPt s,T[[i,l]]] [[j]]]]]] , 

{j ,1, Length [StartPt s]-3>] , {i , 1 , Length [T] }] ; 
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q=Infimty; 

test=Table [True , {i , 1 , Length [dist [ [1] ] 1 >3 ; 

Do[Which[ 

And [Table [dist [[j , ill >=0 , {i , 1 , Length [dist [ [ j] ] ] >1 - test 

Min  [dist  [[j]]]<q]  , 

q=Min[dist[[j]]] ;]  ,{j ,1 .Length [dist]}] ; 

AppendTo [L , {{T [ [Posit ion [dist , q] [[1,1]] ,1,3]] , 

T [[Position [dist ,q] [[1,1]] ,1,1]]} >q}] > 

ImagePts=Drop [ImagePts , 

{Position[StartPts, 

T[ [Position [dist  ,q]  [[1]]  [[1]]  ,1,2]]]  [[1,1]]}]  , 
StartPts=Drop [StartPts , 

{Position [StartPts , 

T [[Position [dist ,q] [[!]] [[!]] ,1,2]]] [[1,1]]}] , 


APPENDIX  D:  COMPUTING  THE  PACKING  LABEL 


One  of  the  fundamental  results  in  the  theory  and  practice  of  circle  packing  is  that 
given  any  assignment  of  positive  numbers  to  the  boundary  vertices  of  a  triangulation 
/C,  there  exists  an  essentially  unique  locally  univalent  circle  packing  for  K.  whose 
boundary  circles  have  these  numbers  as  their  radii  [42].  This  essentially  unique  circle 
packing  is  determined  by  the  so-called  packing  label,  a  list  of  the  vertices  in  the 
triangulation  K  and  an  associated  radius  for  each.  The  first  difficulty  in  programming 
a  circle  packing  algorithm  is  the  calculation  of  these  radii. 

The  following  Mathematica  function  is  an  implementation  of  the  law  of  cosines 
giving  the  angle  determined  by  two  legs  of  a  triangle  created  by  connecting  the  centers 
of  three  mutually  tangent  circles.  The  values  x,y,and  z  represent  the  radii  of  the  three 
circles.  The  angle  calculated  is  the  angle  at  the  center  of  the  circle  with  radius  x. 

T[x_,y_,z_] :=N[ArcCos[((x+y)+(x+z)-(y+z))/(2(x+y) (x+z))]] 

The  Mathematica  function  below  takes  a  description  of  the  complex  K,  and  a 
list  of  boundary  vertices  along  with  an  assigned  radius  for  each  boundary  vertex, 
and  it  returns  the  essentially  unique  Euclidean  packing  label  associated  with  those 
inputs.  The  algorithm  which  this  function  implements  was  suggested  by  the  meta¬ 
code  described  in  Practicum  I  of  [44]. 

PackingLabel [K_ , B_] :=  Block[{e,n,m,R,A,o,G,c,t}, 
e=. 00001; 
n=Length [K] ; 
m=Length[B] ; 

R=Table[l. 1.  n}] ; 

Do[R[[B[[i,l]]]]=B[[i,2]]  ; 

A=Complement[Table[i,-[:Lj:L.nKl , Table [B[[i,l]]  ; 

o=n-m; 

G=Table [2*Pi , {i , 1 , o}] ; 

Do  [G[[i]]  = 

T[R[[K[[A[[i]]  ,1]]]]  , R [ [Last [K [ [A [ [i] ]  ,2]]]]]  , 
R[[First[K[[A[[i]],2]]]]]]  + 

Sum[T[R[[K[[AC[i]l  ,1]]]]  ,R[[K[[A[[i]]  ,2,  j]]]]  , 

R[[K[[A[[i]]  ,2,  j+1]]]]]  , 

{ 1, Length [K[ [A [[i]] ,2]]]-l>] , 

{i,l,o}] ; 
c=2; 

t=Max [Table [Abs [2*Pi-G [ [i] ] ] , {i , 1 , o>] ] ; 

While [t>e , 

Do  [If  [G  [  [i]  ]  <2*Pi  ,R[[A[[i]]]]=R[[A[[i]]]]-l/c, 
R[[A[[i]]]]=R[[A[[i]]]]+l/c]  , 

{i,l,0}]  ; 
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Do [If [E [ [A [ [i] ]  ]  ]  <0 ,  R  [  [A  [ [i] ] ] ] =N [1/ (2c) ] ] ,{i,l,o>] ; 
c=c+l; 

Do  [G[[i]]  = 

T[R[[K[[A[[i]]  ,1]]]]  , R [ [Last  [K [ [A [ [i] ]  ,2]]]]]  , 
R[[First[K[[A[[i]],2]]]]]]  + 
Sum[T[R[[K[[A[[i]]  ,1]]]]  ,R[[K[[A[[i]]  ,2,  j]]]]  , 
R[[K[[A[[i]]  ,2,  j+1]]]]]  , 

{ j,  1, Length [K[ [A C [i] ]  ,2]]]-l>]  , 

{i,l,o}]  ; 

t=Max [Table [Abs [2*Pi-G [ [i] ] ] , {i , 1 , o>] ] ; ]  ; 


The  key  to  this  algorithm  is  the  adjustment  of  the  radii  at  each  step.  In  the 
implementation  shown,  we  add  and  subtract  terms  in  a  harmonic  series  as  appropriate 
to  correct  the  angle  sum;  a  correction  step  is  included  to  prevent  negative  radii. 
This  method,  chosen  for  simplicity,  is  quite  slow.  There  are,  however,  many  other 
algorithms  for  the  computation  of  packing  labels  [17,  37,  44]. 


Figure  D.l:  A  Triangulation  K 


In  Figure  D.l,  we  see  a  triangulation  K  with  20  vertices.  If  we  give  each  of  the 
boundary  vertices  a  positive  radius,  in  order  to  achieve  a  locally  univalent  packing 
we  need  to  compute  the  remaining  radii.  We  assign  each  boundary  vertex  a  radius  of 
1,  and  we  use  the  algorithm  described  to  generate  a  complete  list  of  the  radii  in  the 
packing,  given  in  the  second  column  of  Table  D.l.  As  stated  above,  the  algorithm 
is  very  slow;  448397  iterations  were  required  to  achieve  results  to  within  10“5.  If  we 
assign  different  radii  to  the  boundary  vertices,  we  obtain  a  different  set  of  values  for 
the  radii  of  the  interior  vertices.  That  is,  the  packing  label  is,  as  we  have  already 
stated,  unique.  For  example,  if  we  assign  a  radius  of  1  to  the  vertices  1,11,13,  and 


129 


20  and  a  radius  of  |  to  the  vertices  2,4,7,10,12,14,17,  and  19  we  obtain  the  packing 
label  given  in  the  third  column  of  Table  D.l. 


Table  D.l:  Computed  Packing  Labels  for  /C 


Vertex 

Computed  Radius  for  B1 

Computed  Radius  for  jB2 

1 

1 

1 

1 

2 

0.663959 

0.377904 

1 

1 

2 

0.767825 

0.394196 

6 

0.502727 

0.253292 

7 

1 

1 

2 

8 

0.409361 

0.205468 

9 

0.493900 

0.247294 

10 

1 

1 

2 

11 

1 

l 

12 

1 

1 

2 

13 

1 

l 

14 

1 

1 

2 

15 

0.461162 

0.231158 

16 

0.645202 

0.324809 

17 

1 

1 

2 

18 

0.154701 

0.093836 

19 

1 

1 

2 

20 

1 

1 

As  we  have  already  stated,  there  are  other  algorithms  by  which  we  may  calcu¬ 
late  the  packing  label.  While  the  algorithm  given  above  is  certainly  convergent,  its 
convergence  rate  is  extremely  slow.  The  Mathematica  command  PackingLabel2[], 
given  below,  takes  as  its  input  the  same  arguments  as  the  PackingLabelfl,  but  the 
convergence  is  significantly  faster.  The  radii  computed  above  (requiring  448397  passes 
through  the  algorithm  coded  in  PackingLabelQ)  are  recalculated  using  the  algorithm 
coded  in  PackingLabel2[j,  requiring  only  96  passes  through  the  modified  algorithm. 
In  this  new  algorithm,  rather  than  changing  the  radius  associated  with  every  interior 
vertex  on  each  pass  through  the  algorithm,  here  we  make  only  one  change  with  every 
pass  through  the  algorithm.  Also,  rather  than  changing  the  radii  by  an  amount 
determined  by  a  convergent  alternating  series,  at  each  step  we  solve  the  equation 

n 

2^  =  5>M'  (D-1) 

where  r  is  the  radius  of  the  circle  at  the  center  of  a  flower  and  9i(r)  is  the  central 
angle  of  a  triangle  created  by  joining  the  legs  of  three  adjacent  circles  in  the  circle 
packing.  This  is  illustrated  in  Figure  D.2. 
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Figure  D.2:  Angle  Sums 


PackingLabel2[K_,B_3  :=  Block[{e,n,m,R,A,o,G,c,t ,d,r,i, j  ,a>, 
e= .  00001 ;  n=Length  [K]  ;  m=Length  [B]  ; 

R=Table[l. ,{i,  1,  n}3  ; 

Do[R[[B[[i,l3333=B[[i,233  ,-[i,l,m}3; 

A=Complement [Table [i,{i,l,n>]  , Table [B[[i,l33  ,{i,l,m>33  ; 
o=n~m; 

G=Table [2*Pi , {i , 1 , o>] ;  Do [G [ [i] ] = 

T[R[[K[[A[[i33  ,1]]]]  .RCCLast  [K[[A[[i]]  ,2]]]]]  , 

R[ [First  [K[[A[[i]]  ,2333333+ 

Sum [T [R [ [K [ [A [ [i] ]  ,11333  ,R[[K[[A[[i33  ,2,  J333D  , 
RCtKLtACLill ,2,j+l33333 , 

{1 , Length [K[ [A [[i]] ,2333-133 , 

{i , 1 , o>3 ; 

t=Max [Table [Abs [2*Pi— G [ [i3  33 , {i » 1 , °33  3  > 


While[t>e, 

d=Position[Table[Abs[2*Pi-G[[i333  ,{i,l,o}3 ,t3  L H » 1 J J > 

R[[K[[A[[d33  , 13333= 

FindRoot [2*Pi==T [r ,  R  [ [Last [K [ [A [ [d3  3  ,233333  , 
r[ [First  [K[[A[[d33  ,2333333  + 

Sum[T[r ,R[[K[[A[[d33  ,2,  j3333  ,R[[K[[A[[d33  ,2,  j+133333  , 
{j  ,l,Length[K[[A[[d33  , 2333-133 , 

•Cr,R[[K[[A[[d33 ,13333  ,e,Infinity>3  [[1,233  ; 

D°CGT[R[3[K[[A[[i]3  ,13333  , R [ [Last [K [ [A [ [i3 3  ,233333  , 

R[[First [K[[A[[i33  ,2333333  + 

Sum[T[R[[K[[A[[i33 ,13333  ,R[[K[[A[[i33 ,2, J3333  , 
R[[K[[A[[i33, 2,j+l33333  , 

{j ,1 ,Length[K[[A[[i33 ,2333-133 , 

{i,l,o>3; 

t=Max [Table [Abs [2*Pi-G [ [i3  3  3  ,-Ci,l,o>33  ;3  ;R3 


We  note  that  both  PackingLabelj]  and  PackingLabel2[]  can  accommodate 
and  compute  valid  packing  labels  for  triangulations  which  have  no  boundary  vertices. 
That  is,  we  may  compute  packing  labels  associated  with  triangulations  on  surfaces. 
For  example,  Figure  D.3  shows  a  simple  triangulation  on  a  torus. 


Figure  D.3:  A  Triangulation  on  a  Torus 

The  packing  label  for  this  triangulation  gives  radius  2.4142  to  the  vertices  labeled 
1, 2, . . . ,  12  and  a  radius  of  1  to  the  remaining  vertices  (from  the  initialization  values 
used  in  our  algorithms). 
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APPENDIX  E:  LAYING  OUT  THE  CIRCLE  PACKING 


Given  a  packing  label  R  associated  with  a  complex  1C,  we  wish  to  construct  a 
graphical  representation  of  the  essentially  unique  circle  packing  they  induce.  This 
process  is  not  a  numerical  process;  we  merely  lay  out  the  circles  in  their  respective 
positions  in  the  complex  plane  C.  This  is  achieved  through  an  accumulation  process 
implemented  in  the  Mathematica  function  Layout  [].  This  function  takes  as  its  argu¬ 
ments  the  variables  K,  R,  B,  a,  and  g.  K  is  a  description  of  the  complex  K,  as  a  list 
of  vertices  and  their  neighbors  in  the  form 


Note  that  each  vertex  is  simply  identified  as  an  integer  1,2 The  vertices 
{u^,  v\, . . . ,  v\. }  surrounding  the  ith  must  be  listed  in  counterclockwise  order  to  ensure 
that  the  orientation  remains  consistent.  R  is  a  row  vector  listing,  in  order,  the  radii 
of  each  vertex  in  the  packing.  B  is  a  subset  of  {1, 2, . . . ,  n}  naming  the  boundary 
vertices.  This  information  would  be  enough  to  generate  a  packing,  unique  up  to 
rotation  and  translation  (i.e.,  essentially  unique).  The  remaining  parameters,  a  and 
g,  are  elements  of  {1,2,...,  n}  that  place  the  center  of  a  circle  in  the  packing  at  the 
origin  (the  circle  associated  with  the  vertex  a)  and  place  the  center  of  an  adjacent 
circle  in  the  packing  on  the  real  axis  (the  circle  associated  with  the  vertex  g).  The 
function  determines  the  unique  set  of  centers  for  the  circles  in  the  packing  under  the 
normalization  imposed  by  a  and  g.  This  algorithm  was  suggested  by  meta-code  given 
in  Practicum  I  from  [44]. 


Layout [K_ , R_ , B_ , a_ , g_] : =Block [{n ,S,P,L,Q,q,u,i,j,b,m,l,p}, 
n=Length[R] ;  S=Table [Circle [{0,  0},R[[i]]] ,{i,  1,  n>] ; 

S [ [g] ] =Cir cle [{R [ [a] ] +R [ [g] ] , 0} , R [ [g] ] ] ;  P={a , g> ;  While [Length [P] <n , 
L=Complement [Table [i , {i , 1  ,n>] ,P] ; 

Do[ 

Q=Table[{L[[j]]  ,K[[L[[j]]  ,2,i]]  ,K[[L[[j]]  ,2,i+l]]}, 

{i, 1, Length [K[[L[[j]] ,2]]]-l>] ; 

If [Not [MemberQ [B ,  L  [  [  j  ]  ]  ]  ]  , AppendTo [Q , 

{L[[j]]  ,Last  [K[[L[[j]]  ,2]]]  , First [K[[L[[j]]  ,2]]]>]  ;]  ; 
u=Length[Q]  ; 
i=l; 
q=0; 

While [And [i<=u , q===0] , 

If [And [MemberQ [P,Q[[i, 2]]] , MemberQ [P,Q[[i, 3]]]] , 
b=T [R [ [Q [ [i , 2] ] ] ] ,R[[Q[[i,3]]]] ,  R[[Q[[i,l]]]]] ; 
m=N [S [ [Q [ [i , 2] ] ] ]  [[l,l]]+I*S[[Q[[i,2]]]]  [[1,2]]]  ; 

1=N  [S  [  [Q  [  [i ,  3]  ]  ]  ]  [  [1 , 1]  ]  +I*S  [  [Q  [  [i ,  3]  ]  ]  ]  [  [1 , 2]  ]  ]  ; 

If [m===0 . , 

p=  (R [ [Q  [ [i ,  1]  ]  ]  ]  +R [ [Q  [ [i , 2] ] ] ] )  *Exp  [I*  (Arg [1]  +b) ]  ; 
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S[[L[[j]]]]=Circle[{Re[p]  ,Im[p]}>R[[L[[j]]]]]  ; 

P=Sort [AppendTo  [P ,  L  [  [  j  3  ]  ]  ]  ; 
q— 1  >  > 

p=  (R  [  CQ  [  [i ,  1]  ]  3  3  +R  C  [Q  [  [i ,  2]  ]  ]  ]  )  *Exp  [I*  ( Arg  [1-m]  +b)  ]  +m ; 
S  [  [L  [  [ j ]  ]  ]  ]  =Cir cle  [{Re  [p]  ,  Im [p]  > ,  R  [  [L  [  [ j]  ]  ]  ]  ]  ; 

P=Sort [AppendT o[P,L[[j]]]]; 

q  =  l; 

3, 

i=i+l;3 

3, 

{j ,1, Length [L]>] ; 

];S] 


Figure  E.l:  Distinct  Circle  Packings  Induced  by  Different  Boundary  Radii 


As  an  example  of  this  packing,  consider  the  complex  K  described  in  Figure  D.l 
and  the  two  sets  of  associated  radii  given  in  Table  D.l.  If  we  let  a  =  1  and  g  =  2  in 
each  case,  the  function  Layout []  gives  the  arrangements  of  circles  shown  in  Figure 
E.l. 

We  may  also  lay  out  a  packing  on  a  surface.  Consider  the  complex  described  in 
Figure  D.3  and  the  associated  radii  computed  in  Appendix  E.  As  when  we  “cut  open” 
a  surface  along  generators  of  a  fundamental  group  in  order  to  lay  out  a  fundamental 
region  for  that  surface,  we  open  the  complex  by  removing  tangencies,  then  use  the 
Layout  []  (with  a  =  1  and  g  =  2)  to  arrange  the  circle  packing  in  the  plane.  The 
result  is  shown  in  Figure  E.2. 
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